Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,039,4670GA73.0% 36.7 / 53.2 74noncoding (1949/2905 nt)rrlA23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (9/11);  new base A (32/22);  total (41/33)
Fisher's exact test for biased strand distribution p-value = 3.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.06e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGACGTATACGGTGTGACGCCTGCCCGGTGCCGGAAGGTTAATTGATGGGGTTAGCCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCACGAATGGCGTAATGATGGCCAGGCTGTCTCCACCCGAGACTCA  >  NC_000913/4039330‑4039528
                                                                                                                                          |                                                             
cgaagTATACGGGCTGACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCagg                                                             >  1:1179964/5‑140 (MQ=2)
 gaagTATACGGGCTGACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGt                                                            <  1:615947/137‑1 (MQ=2)
 gaagTATACGGGCTGACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGt                                                            <  1:115871/137‑1 (MQ=2)
 gaagTATACGGGCTGACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGt                                                            <  2:883764/137‑1 (MQ=2)
 gaagTATACGGGCTGACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGt                                                            <  1:791521/137‑1 (MQ=2)
 gaagTATACGGGCTGACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGt                                                            <  1:508100/137‑1 (MQ=2)
 gaagTATACGGGCTGACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGt                                                            <  2:226459/137‑1 (MQ=2)
 gaagTATACGGGCTGACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGt                                                            <  2:403318/137‑1 (MQ=2)
 gaagTATACGGGCTGACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGt                                                            <  1:180945/137‑1 (MQ=2)
 gACGTATTCGGTGTGACGCCTGCCCGGTGCCGGAAGGTTAATTGATGGGGTTAGCCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGt                                                            <  2:1579715/140‑1 (MQ=255)
 gACGTATACGGTGTGACGCCTGCCCGGTGCCGGAAGGTTAATTGATGGGGTTAGCCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGt                                                            <  1:1612933/140‑1 (MQ=40)
 gACGTATACGGTGTGACGCCTGCCCGGTGCCGGAAGGTTAATTGATGGGGTTAGCCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGt                                                            <  1:1646036/140‑1 (MQ=40)
  aagTATACGGGCTGACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTa                                                           >  1:258535/3‑140 (MQ=2)
  aagTATACGGGCTGACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTa                                                           >  1:428800/3‑140 (MQ=2)
       taCGGTGTGACGCCTGCCCGGTGCCGGAAGGTTAATTGATGGGGTTAGCCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTc                                                      >  2:1608690/1‑140 (MQ=40)
         cgggctgACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTctg                                                    >  2:523368/6‑138 (MQ=2)
         cgggctgACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTctg                                                    >  2:170352/6‑138 (MQ=2)
          gggctgACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTctga                                                   <  2:111608/136‑4 (MQ=2)
           gtgtGACGCCTGCCCGGTGCCGGAAGGTTAATTGATGGGGTTAGCCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGAc                                                  >  1:1591123/1‑140 (MQ=40)
            tgtgACGCCTGCCCGGTGCCGGAAGGTTAATTGATGGGGTTAGCCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGAcc                                                 >  1:1307996/1‑140 (MQ=40)
            tgtgACGCCTGCCCGGTGCCGGAAGGTTAATTGATGGGGTTAGCCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGAcc                                                 >  2:1528371/1‑140 (MQ=40)
            gctgACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTGAATTCTGAcc                                                 >  2:852736/3‑140 (MQ=2)
            gctgACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGAcc                                                 >  2:1114392/3‑140 (MQ=2)
            gctgACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGAcc                                                 <  1:170352/138‑1 (MQ=2)
            gctgACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGAcc                                                 >  2:236326/3‑140 (MQ=2)
             ctgACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGt                                                            >  1:1025832/2‑128 (MQ=2)
             ctgACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGt                                                            <  2:1025832/127‑1 (MQ=2)
             ctgACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGAccc                                                >  1:833887/2‑139 (MQ=2)
             ctgACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGAccc                                                >  2:251651/2‑139 (MQ=2)
             ctgACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGAccc                                                >  1:927598/2‑139 (MQ=2)
             ctgACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGAccc                                                >  2:387754/2‑139 (MQ=2)
             ctgACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGAccc                                                >  2:10872/2‑139 (MQ=2)
             ctgACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGAccc                                                >  2:66036/2‑139 (MQ=2)
             ctgACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCAGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTCTCAGGTAAATTCTGAccc                                                >  2:1153186/2‑139 (MQ=2)
              tgACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGAcccg                                               <  2:492363/140‑3 (MQ=2)
              tgACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGAcccg                                               >  2:529986/1‑138 (MQ=2)
              tgACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGAcccg                                               <  2:643249/140‑3 (MQ=2)
              tgACGCCTGCCCGGTGCCGGAAGGTTAATTGATGGGGTTAGCCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTg                                               <  2:1451309/140‑1 (MQ=40)
               gACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGc                                              >  2:243950/1‑140 (MQ=2)
               gACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGc                                              >  2:228430/1‑140 (MQ=2)
               gACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGc                                              >  2:546776/1‑140 (MQ=2)
               gACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGc                                              >  2:732345/1‑140 (MQ=2)
               gACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGc                                              >  2:996598/1‑140 (MQ=2)
               gACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGc                                              >  1:163935/1‑140 (MQ=2)
               gACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGc                                              >  1:475585/1‑140 (MQ=2)
               gACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGc                                              >  1:434180/1‑140 (MQ=2)
               gACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGc                                              >  1:558351/1‑140 (MQ=2)
               gACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGc                                              >  1:927876/1‑140 (MQ=2)
               gACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGc                                              >  2:1015781/1‑140 (MQ=2)
               gACGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGc                                              >  1:150360/1‑140 (MQ=2)
                aCGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGt                                                            >  2:556013/1‑125 (MQ=2)
                aCGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGt                                                            >  1:924935/1‑125 (MQ=2)
                aCGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGt                                                            <  2:924935/125‑1 (MQ=2)
                aCGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGt                                                            <  1:556013/125‑1 (MQ=2)
                 cGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGCAc                                            >  2:904253/1‑140 (MQ=2)
                 cGCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGCAc                                            >  2:931371/1‑140 (MQ=2)
                  gCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGCACg                                           <  1:237676/140‑1 (MQ=2)
                  gCCTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAG‑‑‑GCGAAATTCCTTGTCAGGTAAATTCTGACCCGCACGAAt                                        <  1:129428/140‑1 (MQ=2)
                   ccTGCCCGGTGCCGGAAGGTTAATTGATGGGGTTAGCCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCACGa                                          >  2:1314418/1‑140 (MQ=40)
                    cTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGCACGaa                                         <  1:1153186/140‑1 (MQ=2)
                    cTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGCACGaa                                         <  2:299612/140‑1 (MQ=2)
                    cTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGCACGaa                                         <  2:258535/140‑1 (MQ=2)
                    cTGCCCGGTGCTGGAAGGTTAAGGGGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGCACGaa                                         <  2:163935/140‑1 (MQ=2)
                     tGCCCGGTGCTGGAAGGTTAAGGTGAGCGGTTAGGGGTAACCCGAAGCTGTGAACCGAAGCCCCAGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCAGGTAAATTCTGACCCGCACGAAt                                        <  2:1201807/140‑1 (MQ=2)
                             gCCGGAAGGTTAATTGATGGGGTTAGCCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCACGAATGGCGTAat                                <  1:1608690/140‑1 (MQ=40)
                                   aGGTTAATTGATGGGGTTAGCCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCACGAATGGCGTAATGATGGc                          <  2:1764142/140‑1 (MQ=40)
                                   aGGTTAATTGATGGGGTTAGCCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCACGAATGGCGTAATGATGGc                          <  2:1293364/140‑1 (MQ=40)
                                    ggTTAATTGATGGGGTTAGCCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCACGAATGGCGTAATGATGGcc                         <  2:1725775/140‑1 (MQ=40)
                                      ttAATTGATGGGGTTAGCCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCACGAATGGCGTAATGATGGCCAg                       >  1:1508219/1‑140 (MQ=40)
                                      ttAATTGATGGGGTTAGCCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCACGAATGGCGTAATGATGGCCAg                       >  2:1643283/1‑140 (MQ=40)
                                       tAATTGATGGGGTTAGCCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCACGAATGGCGTAATGATGGCCAgg                      <  1:1405031/140‑1 (MQ=40)
                                              aTGGGGTTAGCCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCACGAATGGCGTAATGATGGCCAGGCTGtct                <  1:1528371/140‑1 (MQ=40)
                                                ggggTTAGCCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCACGAATGGCGTAATGATGGCCAGGCTGTCTcc              <  2:1338194/140‑1 (MQ=40)
                                                            aaGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCACGAATGGCGTAATGATGGCCAGGCTGTCTCCACCCGAGACTCa  >  1:1284364/1‑140 (MQ=32)
                                                            aaGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCACGAATGGCGTAATGATGGCCAGGCTGTCTCCACCCGAGACTCa  >  2:1505166/1‑140 (MQ=32)
                                                                                                                                          |                                                             
GGACGTATACGGTGTGACGCCTGCCCGGTGCCGGAAGGTTAATTGATGGGGTTAGCCGCAAGGCGAAGCTCTTGATCGAAGCCCCGGTAAACGGCGGCCGTAACTATAACGGTCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCACGAATGGCGTAATGATGGCCAGGCTGTCTCCACCCGAGACTCA  >  NC_000913/4039330‑4039528

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