Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913224,9890AT64.5% 77.2 / 112.1 93noncoding (1219/1542 nt)rrsH16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (10/23);  new base T (27/33);  total (37/56)
Fisher's exact test for biased strand distribution p-value = 1.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.02e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATACAAAGAGAA  >  NC_000913/224851‑225027
                                                                                                                                           |                                      
aaTGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAg                                       <  2:2462972/140‑1 (MQ=32)
aaTGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAg                                       <  2:1678432/140‑1 (MQ=32)
aaTGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAg                                       <  1:1653281/140‑1 (MQ=32)
 aTGTTGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTgg                                      >  1:901132/1‑140 (MQ=2)
    ttGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCt                                   >  1:507179/1‑140 (MQ=2)
     tGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTa                                  <  1:1250158/140‑1 (MQ=32)
     tGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTa                                  <  2:176353/140‑1 (MQ=2)
     tGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTa                                  <  2:412406/140‑1 (MQ=2)
     tGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTa                                  <  2:462778/140‑1 (MQ=2)
     tGGGTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTa                                  <  1:822668/140‑1 (MQ=2)
      gggTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTac                                 <  2:2176125/140‑1 (MQ=32)
      gggTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTac                                 <  2:311481/140‑1 (MQ=2)
      gggTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTac                                 >  2:368955/1‑140 (MQ=2)
      gggTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTac                                 <  1:870408/140‑1 (MQ=2)
      gggTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTac                                 <  2:824346/140‑1 (MQ=2)
       ggTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTaca                                >  2:1789018/1‑140 (MQ=33)
       ggTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAGGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTaca                                >  1:655803/1‑140 (MQ=11)
       ggTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTaca                                <  2:173230/140‑1 (MQ=2)
        gTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTacac                               <  2:2590446/140‑1 (MQ=32)
        gTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacac                               <  2:82380/140‑1 (MQ=2)
        gTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacac                               <  2:373725/140‑1 (MQ=2)
        gTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacac                               <  1:837310/140‑1 (MQ=2)
        gTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacac                               <  2:732687/140‑1 (MQ=2)
        gTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacac                               <  2:901132/140‑1 (MQ=2)
        gTTAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacac                               <  1:756122/140‑1 (MQ=2)
         ttAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTacaca                              <  1:1273315/140‑1 (MQ=32)
         ttAAGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGGCAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacaca                              <  1:364019/140‑1 (MQ=2)
           aaGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACg                            <  1:1644618/140‑1 (MQ=32)
           aaGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACg                            <  2:1974190/140‑1 (MQ=32)
           aaGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACg                            <  1:2430750/140‑1 (MQ=32)
           aaGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTCACGACCAGGGCTACACACg                            <  2:1657367/140‑1 (MQ=17)
           aaGTCCCGCAACGAGCGCAACCCTTATCTTTTGTCGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACg                            >  2:1917565/1‑140 (MQ=18)
           aaGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACg                            >  1:22817/1‑140 (MQ=2)
           aaGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACg                            >  1:979735/1‑140 (MQ=2)
            aGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGt                           <  2:1494758/140‑1 (MQ=32)
            aGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGt                           >  1:1283330/1‑140 (MQ=33)
            aGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGt                           >  1:1349088/1‑140 (MQ=33)
            aGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACTAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGt                           >  2:618998/1‑140 (MQ=2)
            aGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGt                           >  1:511647/1‑140 (MQ=2)
            aGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGt                           >  1:13538/1‑140 (MQ=2)
            aGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTac                                 >  1:41245/1‑134 (MQ=2)
            aGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTac                                 <  2:41245/134‑1 (MQ=2)
            aGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTac                                 <  2:749691/134‑1 (MQ=2)
            aGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTac                                 >  1:749691/1‑134 (MQ=2)
            aGTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGt                           >  2:60058/1‑140 (MQ=2)
             gTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGta                          <  1:199175/140‑2 (MQ=2)
             gTCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGta                          <  2:269586/140‑2 (MQ=2)
              tCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGc                         >  1:1231453/1‑140 (MQ=33)
              tCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCAAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGc                         >  2:1469315/1‑140 (MQ=17)
              tCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacac                               >  1:804148/1‑134 (MQ=2)
              tCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacac                               <  2:804148/134‑1 (MQ=2)
              tCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGta                          <  2:46362/139‑2 (MQ=2)
              tCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGta                          >  1:46362/1‑138 (MQ=2)
              tCCCGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGtac                         >  1:50373/1‑138 (MQ=2)
                 cGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacac                               <  1:852900/131‑1 (MQ=2)
                 cGCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTTGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacac                               >  2:852900/1‑131 (MQ=2)
                  gCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACa                     <  2:1603113/140‑1 (MQ=32)
                  gCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACa                     <  2:1626033/140‑1 (MQ=32)
                  gCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACa                     <  1:1375666/140‑1 (MQ=32)
                  gCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACa                     <  1:1436844/140‑1 (MQ=32)
                  gCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACa                     <  2:428408/140‑1 (MQ=2)
                  gCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACa                     >  2:646986/1‑140 (MQ=2)
                  gCAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACa                     <  1:859740/140‑1 (MQ=2)
                   cAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacac                               <  1:642751/130‑1 (MQ=2)
                   cAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACGTAAGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTacac                               >  2:642751/1‑130 (MQ=2)
                   cAACGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACaa                    >  2:276002/1‑140 (MQ=2)
                    aaCGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAAt                   <  1:1193249/140‑1 (MQ=32)
                    aaCGAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAAt                   <  2:474432/140‑1 (MQ=2)
                      cGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATgg                 <  1:1098351/140‑1 (MQ=32)
                      cGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATgg                 <  2:2138688/140‑1 (MQ=32)
                       gAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGgc                <  2:1073572/140‑1 (MQ=32)
                       gAGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAGGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgc                <  2:655803/140‑1 (MQ=2)
                        aGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGgcg               >  1:1620484/1‑140 (MQ=33)
                        aGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGgcg               <  2:2401635/140‑1 (MQ=32)
                        aGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGgcg               >  2:1014918/1‑140 (MQ=33)
                        aGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGta                          <  1:799215/129‑2 (MQ=2)
                        aGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGta                          >  2:799215/1‑128 (MQ=2)
                        aGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgcc               >  2:553398/1‑139 (MQ=2)
                        aGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgcc               >  1:504031/1‑139 (MQ=2)
                        aGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgcc               >  1:287132/1‑139 (MQ=2)
                        aGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgcc               >  2:395014/1‑139 (MQ=2)
                        aGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgcc               >  2:657882/1‑139 (MQ=2)
                        aGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgcc               >  2:287089/1‑139 (MQ=2)
                        aGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgc                >  1:792854/1‑139 (MQ=2)
                        aGCGCAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgc                <  2:792854/139‑1 (MQ=2)
                          cgcAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCa             >  1:1933861/1‑140 (MQ=33)
                          cgcAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgccag             <  2:287132/140‑4 (MQ=2)
                           gcAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgccagt            <  1:553398/140‑5 (MQ=2)
                            cAACCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAGAAAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgccagta           <  2:504031/140‑6 (MQ=2)
                              aCCCCTATTTTTAGTTGCCAGCACTTCGGGTGGGCACTCTAAAGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTACTACAATGgccagtaca         <  2:686442/140‑8 (MQ=2)
                               cccTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATACaa        <  1:1831309/140‑1 (MQ=32)
                                     tCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATACAAAGAGaa  <  1:1469829/140‑1 (MQ=32)
                                     tCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATACAAAGAGaa  >  2:2251126/1‑140 (MQ=32)
                                                                                                                                           |                                      
AATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATACAAAGAGAA  >  NC_000913/224851‑225027

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: