Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,039,4770CT74.0% 30.8 / 52.1 77noncoding (1959/2905 nt)rrlA23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (5/15);  new base T (37/20);  total (42/35)
Fisher's exact test for biased strand distribution p-value = 3.51e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGGTGTGACGCCTGCCCGGTGCCGGAAGGTTAATTGATGGGGTTAGCCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCACGAATGGCGTAATGATGGCCAGGCTGTCTCCACCCGAGACTCAGT  >  NC_000913/4039339‑4039530
                                                                                                                                          |                                                     
cGGTGTGACGCCTGCCCGGTGCCGGAAGGTTAATTGATGGGGTTAGCCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCg                                                      <  2:1799575/140‑1 (MQ=40)
 gggctgACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGa                                                     <  2:644731/136‑1 (MQ=2)
   gctgACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGAcc                                                   >  2:128990/3‑140 (MQ=2)
   gctgACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGAcc                                                   >  1:1407004/3‑140 (MQ=2)
   gctgACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGAcc                                                   <  1:89270/138‑1 (MQ=2)
   gctgACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGAcc                                                   <  1:337601/138‑1 (MQ=2)
    ctgACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGa                                                     >  1:1123798/2‑137 (MQ=2)
    ctgACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGAccc                                                  >  2:1495522/2‑139 (MQ=2)
    ctgACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGAccc                                                  >  2:754249/2‑139 (MQ=2)
     tgACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGAcccg                                                 <  2:1801021/140‑3 (MQ=2)
     tgACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGAcccg                                                 <  2:1091069/140‑3 (MQ=2)
     tgACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGAcccg                                                 <  1:1947803/140‑3 (MQ=2)
      gACGCCTGCCCGGTGTTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGc                                                >  2:128160/1‑140 (MQ=2)
      gACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTTAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGc                                                >  2:1087763/1‑140 (MQ=2)
      gACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGAcc                                                   >  1:1006883/1‑137 (MQ=2)
      gACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGAcc                                                   <  2:1006883/137‑1 (MQ=2)
      gACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGc                                                >  2:737343/1‑140 (MQ=2)
      gACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGc                                                >  2:1043810/1‑140 (MQ=2)
      gACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGc                                                >  2:680091/1‑140 (MQ=2)
      gACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGc                                                >  2:1003164/1‑140 (MQ=2)
      gACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGc                                                >  2:1053332/1‑140 (MQ=2)
      gACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGc                                                >  1:894635/1‑140 (MQ=2)
      gACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGc                                                >  1:1548572/1‑140 (MQ=2)
      gACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGc                                                >  1:123924/1‑140 (MQ=2)
      gACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGc                                                >  1:810663/1‑140 (MQ=2)
      gACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGc                                                >  2:1067762/1‑140 (MQ=2)
      gACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGc                                                >  1:358576/1‑140 (MQ=2)
      gACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGc                                                >  1:1483088/1‑140 (MQ=2)
      gACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGc                                                >  2:211882/1‑140 (MQ=2)
      gACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGc                                                >  2:1088515/1‑140 (MQ=2)
      gACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGc                                                >  2:727257/1‑140 (MQ=2)
      gACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGc                                                >  2:1208510/1‑140 (MQ=2)
      gACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGc                                                >  1:1940612/1‑140 (MQ=2)
      gACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGc                                                >  2:48769/1‑140 (MQ=2)
      gACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGc                                                >  2:1274062/1‑140 (MQ=2)
      gACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGc                                                >  2:494232/1‑140 (MQ=2)
      gACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGc                                                >  1:1335936/1‑140 (MQ=2)
      gACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCAGc                                                >  2:1438683/1‑140 (MQ=2)
      gACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACACGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGc                                                >  1:365763/1‑140 (MQ=2)
       aCGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGCa                                               >  1:218764/1‑140 (MQ=2)
       aCGCCTGCCCGGTGCCGGAAGTCTAATTGATGGGGTTAGCCGCAAGGCGAAGCCCTTGATCGAAGCCCCGGTAAGCGGCGGCCGTAACTATAACGATCCTTAGGTAGCGAAATTCCTTGTCGGTATAGTTCCGatccgct                                               >  2:2310370/1‑134 (MQ=14)
       aCGCCTGCCCGGTGCCGGAAGGTTAATTGATGGGGTTAGCCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCa                                               >  1:1920853/1‑140 (MQ=40)
       aCGCCTGCCCGGTGCCGGAAGGTTAATTGATGGGGTTAGCCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCa                                               >  1:1793538/1‑140 (MQ=40)
        ctccTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGCAc                                              >  1:938435/3‑140 (MQ=2)
        cGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTGAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGCAc                                              >  1:1476284/1‑140 (MQ=2)
        cGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGCAc                                              >  2:2240882/1‑140 (MQ=2)
        cGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGCAc                                              >  2:1118924/1‑140 (MQ=2)
        cGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGCAc                                              >  2:1506967/1‑140 (MQ=2)
        cGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGCAc                                              <  1:1246400/140‑1 (MQ=2)
         gCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGTCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGCACg                                             <  2:1999310/140‑1 (MQ=2)
         gCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGCACg                                             <  1:1495522/140‑1 (MQ=2)
          cttGCCGGGTGCCTGAAGGTTAATTGATGGGGTTAGCCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGGAGCGTAATTCCTTGTCGGGTAAGTTCCGACCTGCACGa                                            <  2:2110298/138‑1 (MQ=14)
          ccTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGCACGa                                            <  2:1630136/140‑1 (MQ=2)
          ccTGCCCGGTGCCGGAAGGTTAATTGATGGGGTTAGCCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAATCGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCACGa                                            <  2:1584780/140‑1 (MQ=255)
           cTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGTGGTAACCCGAAGCTGTTAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGCACGaa                                           <  2:1217990/140‑1 (MQ=2)
           cTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGCACGaa                                           <  2:810663/140‑1 (MQ=2)
           cTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGCACGaa                                           <  1:1092688/140‑1 (MQ=2)
           cTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGCACGaa                                           <  1:184440/140‑1 (MQ=2)
           cTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGCACGaa                                           <  2:1194955/140‑1 (MQ=2)
           cTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGCACGaa                                           >  1:862917/1‑140 (MQ=2)
           cTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTATGACCCGCACGaa                                           <  2:824440/140‑1 (MQ=2)
           cTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTACTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGCACGaa                                           <  1:1352337/140‑1 (MQ=2)
           cTGCCCGGTGCCGGAAGGTTAATTGATGGGGTTAGCCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCACGaa                                           <  2:2267229/140‑1 (MQ=40)
            tGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGCACGAAt                                          <  2:1742508/140‑1 (MQ=2)
            tGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCAGCACGAAt                                          <  1:1438683/140‑1 (MQ=2)
            tGCCCGGTGCCGGAAGGTTAATTGATGGGGTTAGCCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCACGAAt                                          <  2:2318596/140‑1 (MQ=40)
                cGGTGCCGGAAGGTTAATTGATGGGGTTAGCCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCACGAATGGCg                                      <  1:1845371/140‑1 (MQ=40)
                cGGTGCCGGAAGGTTAATTGATGGGGTTAGCCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCACGAATGGCg                                      <  1:1925484/140‑1 (MQ=40)
                   tGCCGGAAGGTTAATTGATGGGGTTAGCCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCACGAATGGCGTaa                                   >  2:2144260/1‑140 (MQ=40)
                   tGCCGGAAGGTTAATTGATGGGGTTAGCCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCACGAATGGCGTaa                                   >  1:1757589/1‑140 (MQ=40)
                    gCCGGAAGGTTAATTGATGGGGTTAGCCGCAAGGCGAAGCTCTTGATCGAGGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCACGAATGGCGTAat                                  <  2:2270141/140‑1 (MQ=255)
                     cgggAAGGTTAATTGATGGGGTTAGCGGCAAGGGGAAGATCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTGGCGAAATTCCTTGTGGGGTAAGTTCCGACCTGCGCGAATGGCGTAatg                                 <  2:2277876/138‑1 (MQ=11)
                         aaGGTTAATTGATGGGGTTAGCCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCACGAATGGCGTAATGATgg                             <  1:2129266/140‑1 (MQ=40)
                              tAATTGATGGGGTTAGCCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCACGAATGGCGTAATGATGGCCAgg                        <  1:1892846/140‑1 (MQ=40)
                                    aTGGGGTTAGCCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCACGAATGGCGTAATGATGGCCAGGCTGtct                  <  1:1947330/140‑1 (MQ=40)
                                    aTGGGGTTAGCCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCACGAATGGCGTAATGATGGCCAGGCTGtct                  <  2:1608026/140‑1 (MQ=40)
                                    aTGGGGTTAGCCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCACGAATGGCGTAATGATGGCCAGGCTGtct                  <  2:1793538/140‑1 (MQ=40)
                                              ccGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCACGAATGGCGTAATGATGGCCAGGCTGTCTCCCCCCgaga        <  2:1819538/140‑1 (MQ=25)
                                                gCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCACGAATGGCGTAATGATGGCCAGGCTGTCTCCACCCGAGACt      >  1:1736038/1‑140 (MQ=33)
                                                   aGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCACGAATGGCGTAATGATGGCCAGGCTGTCTCCACCCGAGACTCAg   >  2:2087191/1‑140 (MQ=31)
                                                    ggCGAAGCTCTTGATCGAAGCCCCGGTAGACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCACGAATGGCGTAATGATGGCCAGGCTGTCTCCACCCGAGACTCAGt  <  1:2118960/140‑1 (MQ=11)
                                                                                                                                          |                                                     
CGGTGTGACGCCTGCCCGGTGCCGGAAGGTTAATTGATGGGGTTAGCCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCACGAATGGCGTAATGATGGCCAGGCTGTCTCCACCCGAGACTCAGT  >  NC_000913/4039339‑4039530

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