Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,726,045 C→A 5.6% intergenic (‑299/+24) kgtP ← / ← rrfG alpha‑ketoglutarate transporter/5S ribosomal RNA of rrnG operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,726,0450CA 5.6% 318.3 / 2.7 92intergenic (‑299/+24)kgtP/rrfGalpha‑ketoglutarate transporter/5S ribosomal RNA of rrnG operon
Reads supporting (aligned to +/- strand):  ref base C (43/42);  new base A (2/3);  total (45/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATAAATTTCTGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAA  >  NC_000913/2725898‑2726190
                                                                                                                                                   |                                                                                                                                                   
gATAAATTTCTGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTc                                                                                                                                                     <  1:865682/148‑1 (MQ=255)
gATAAATTTCTGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCa                                                                                                                                                    >  4:23167/1‑149 (MQ=255)
  tAAATTTCTGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑gg                                                                                                                                                 <  3:931540/149‑1 (MQ=255)
    aaTTTCTGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGa                                                                                                                                                <  2:126144/148‑1 (MQ=255)
     aTTTCTGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGAt                                                                                                                                               <  1:844896/148‑1 (MQ=255)
      tttCTGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATgg                                                                                                                                             >  4:868928/1‑149 (MQ=255)
       ttCTGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATgg                                                                                                                                             <  1:393698/148‑1 (MQ=255)
       ttCTGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATgg                                                                                                                                             >  2:550223/1‑148 (MQ=255)
        tCTGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGcc                                                                                                                                           <  3:516252/149‑1 (MQ=255)
         cTGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCTGAC‑GGATGGcc                                                                                                                                           <  2:14998/148‑1 (MQ=255)
          tGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCTGAC‑GGATGGCCtt                                                                                                                                         <  3:71729/149‑1 (MQ=255)
          tGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCtt                                                                                                                                         <  3:724478/149‑1 (MQ=255)
           ggTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCC‑TGACGGATGGCCttt                                                                                                                                        <  4:348647/149‑1 (MQ=38)
            gTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTc                                                                                                                                       >  4:399622/1‑149 (MQ=255)
             tAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCTGAC‑GGATGGCCtttt                                                                                                                                       <  1:809497/148‑2 (MQ=255)
             tAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCC‑TGACGGATGGCCtttt                                                                                                                                       <  2:67017/148‑2 (MQ=38)
               aggaggACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCa                                                                                                                                   >  3:905145/1‑149 (MQ=255)
                   ggaCACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATaa                                                                                                                                >  1:315653/1‑148 (MQ=255)
                       acGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTc                                                                                                                                       <  3:163906/138‑1 (MQ=255)
                       acGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTc                                                                                                                                       >  4:163906/1‑138 (MQ=255)
                        cGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGa                                                                                                                           >  1:448569/1‑148 (MQ=255)
                         gTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATGTAGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGTATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATg                                                                                                                         >  4:793626/1‑149 (MQ=255)
                           aTGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGc                                                                                                                        >  1:640766/1‑148 (MQ=255)
                               aaGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGc                                                                                                                   >  3:65346/1‑149 (MQ=255)
                                aGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCa                                                                                                                  <  4:883322/149‑1 (MQ=255)
                                   gggCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGt                                                                                                                >  2:122787/1‑148 (MQ=255)
                                      cAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTcc                                                                                                             <  2:448569/148‑1 (MQ=255)
                                       aaGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCt                                                                                                           <  4:477811/149‑1 (MQ=255)
                                        aGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTa                                                                                                          <  4:905145/149‑1 (MQ=255)
                                           tGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTAct                                                                                                        >  1:143476/1‑148 (MQ=255)
                                             gggAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTc                                                                                                              >  2:293210/1‑140 (MQ=255)
                                             gggAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTc                                                                                                              <  1:293210/140‑1 (MQ=255)
                                             gggAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTAct                                                                                                        >  1:784094/1‑146 (MQ=255)
                                             gggAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTAct                                                                                                        <  2:784094/146‑1 (MQ=255)
                                             gggAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACtctc                                                                                                     <  4:758292/149‑1 (MQ=255)
                                                 aGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCAt                                                                                                 <  3:18369/149‑1 (MQ=255)
                                                  gCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATg                                                                                                <  3:865425/149‑1 (MQ=255)
                                                   ccGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATg                                                                                                <  1:654985/148‑1 (MQ=255)
                                                   ccGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATg                                                                                                >  2:807211/1‑148 (MQ=255)
                                                    cGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATggg                                                                                              <  4:64995/149‑1 (MQ=255)
                                                           gTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGAccc                                                                                       <  4:398639/149‑1 (MQ=255)
                                                           gTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGAccc                                                                                       >  3:746065/1‑149 (MQ=255)
                                                              gCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCcac                                                                                    >  4:102692/1‑149 (MQ=255)
                                                               cTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCcaca                                                                                   >  4:868907/1‑149 (MQ=255)
                                                                tGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCcaca                                                                                   >  2:135260/1‑148 (MQ=255)
                                                                     cTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTAcc                                                                              >  1:466627/1‑148 (MQ=255)
                                                                     cTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTAcc                                                                              >  1:233611/1‑148 (MQ=255)
                                                                      tGGCATATGTGGGAGTATAAGACGCGCATCGTCGCATCAGGCATTTTCTTCTGCGCCACTGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTTTCGCATGGGGAGACCCCACACTACCAt                                                                            >  4:719709/1‑149 (MQ=255)
                                                                      tGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCAt                                                                            >  4:247031/1‑149 (MQ=255)
                                                                        gCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCt                                                                                                                      <  3:446783/105‑1 (MQ=255)
                                                                        gCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCt                                                                                                                      >  4:446783/1‑105 (MQ=255)
                                                                                gggAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAc                                                                   >  1:200701/1‑148 (MQ=255)
                                                                                gggAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACg                                                                  >  4:833769/1‑149 (MQ=255)
                                                                                 ggAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTc                                                                                                                                       >  3:936567/1‑80 (MQ=255)
                                                                                 ggAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTc                                                                                                                                       <  4:936567/80‑1 (MQ=255)
                                                                                       taAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGttt                                                            <  1:655599/148‑1 (MQ=255)
                                                                                        aaGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCa                                                          >  3:100757/1‑149 (MQ=255)
                                                                                          gACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCa                                                                                                                                   >  1:595912/1‑74 (MQ=255)
                                                                                          gACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCa                                                                                                                                   <  2:595912/74‑1 (MQ=255)
                                                                                           aCGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACt                                                        >  1:208397/1‑148 (MQ=255)
                                                                                           aCGCGCAGCGTCGCATCAGGCATGTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACt                                                        >  1:823203/1‑148 (MQ=255)
                                                                                             gcgcAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCt                                                     <  4:772860/149‑1 (MQ=255)
                                                                                              cgcAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATACAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCt                                                     >  1:793553/1‑148 (MQ=255)
                                                                                              cgcAGCGTCGCATCAGGCATTTTTTTCTGCACCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTg                                                    <  4:901080/149‑1 (MQ=255)
                                                                                                cAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGa                                                   <  2:793553/148‑1 (MQ=255)
                                                                                                   cGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTc                                               <  4:811325/149‑1 (MQ=255)
                                                                                                      cGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCgg                                             >  2:168023/1‑148 (MQ=255)
                                                                                                       gCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGc                                            <  2:200701/148‑1 (MQ=255)
                                                                                                          tCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATgg                                        <  4:100757/149‑1 (MQ=255)
                                                                                                              gCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTc                                               <  2:631810/138‑1 (MQ=255)
                                                                                                              gCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTc                                               >  1:631810/1‑138 (MQ=255)
                                                                                                              gCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTc                                    >  4:178477/1‑149 (MQ=255)
                                                                                                               cATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTTTGAGTTCGGCATGGGGTc                                    >  2:371415/1‑148 (MQ=255)
                                                                                                                aTTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAg                                  <  3:532851/149‑1 (MQ=255)
                                                                                                                  ttttttCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGt                                <  4:388581/149‑1 (MQ=255)
                                                                                                                  ttttttCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGt                                <  3:571874/149‑1 (MQ=255)
                                                                                                                   tttttCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGt                                <  2:937245/148‑1 (MQ=255)
                                                                                                                    ttttCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTg                               <  2:143476/148‑1 (MQ=255)
                                                                                                                    ttttCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTg                               <  1:135260/148‑1 (MQ=255)
                                                                                                                     tttCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTggg                             <  3:178477/149‑1 (MQ=255)
                                                                                                                             ccAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACcg                      <  2:709463/148‑1 (MQ=255)
                                                                                                                                aTGCAAAAAGGCCATCCGTCA‑AGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCgcgc                   >  1:570816/1‑148 (MQ=255)
                                                                                                                                    aaaaaGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACgg              >  4:104239/1‑149 (MQ=255)
                                                                                                                                     gaaaGGCCATCCTGAC‑GGATGGCCTTTTTGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACgg              >  3:511113/2‑149 (MQ=11)
                                                                                                                                       aaGGCCATCCTGAC‑GGATGGCCTTTC‑GCATAATTTGATGCTTGGGAATTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACTACCGCGCTACGgccg           >  3:808017/1‑149 (MQ=14)
                                                                                                                                         ggCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGgccgc          >  1:465812/1‑148 (MQ=255)
                                                                                                                                          gCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCa        <  4:341814/149‑1 (MQ=255)
                                                                                                                                           ccATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCa        >  2:268794/1‑148 (MQ=255)
                                                                                                                                           ccATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAg       <  4:509982/149‑1 (MQ=255)
                                                                                                                                             aTCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAgg      >  2:814500/1‑148 (MQ=255)
                                                                                                                                              tCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGgc             >  1:950213/1‑140 (MQ=255)
                                                                                                                                              tCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGgc             <  2:950213/140‑1 (MQ=255)
                                                                                                                                                cGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCaaa  <  3:102692/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GATAAATTTCTGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAA  >  NC_000913/2725898‑2726190

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑