Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,275,6110AC21.5% 187.6 / 10.8 79E110D (GAA→GACgarDgalactarate dehydratase
Reads supporting (aligned to +/- strand):  ref base A (47/15);  new base C (0/17);  total (47/32)
Fisher's exact test for biased strand distribution p-value = 7.08e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGAAATTATTCGTTATGGCGAAGTGATTGGTTACGCCGTGCGTGCAATCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTG  >  NC_000913/3275473‑3275750
                                                                                                                                           |                                                                                                                                            
tGAAATTATTCGTTATGGCGAAGTGATTGGTTACGCCGTGCGTGCAATCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAAc                                                                                                                                             >  1:2259903/1‑140 (MQ=255)
 gAAATTATTCGTTATGGCGAAGTGATTGGTTACGCCGTGCGTGCAATCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAAcc                                                                                                                                            >  1:244691/1‑140 (MQ=255)
        ttCGTTATGGCGAAGTGATTGGTTACGCCGTGCGTGCAATCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTAccg                                                                                                                                    >  1:2193858/1‑140 (MQ=255)
         tggTTAGGGCAAAGTGTTTGGTTCCGCCGTGCGTGAAACCCCACGGGGAAGCTGGATCGACAAATAAATGTTTGTACTCCCGGAGGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGCCCCGGACCCCTTAccgc                                                                                                                                   <  2:1243644/138‑1 (MQ=255)
            ttATGGCGAAGTGATTGGTTACGCCGTGCGTGCAATCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACcgccgc                                                                                                                                >  1:2270376/1‑140 (MQ=255)
               gggCGAGGTGTTTGTTTACGCGGGGCGTCCAACCCCCCGGGGAAGCTGGATCGACAAATAAATGTTTTTCCTCCGGGAGGCCCCGCCGTTACACCCGCTGCCACTGGCAACCAAAGCCCCGGACCCCTTACCGCCGCTgg                                                                                                                             <  2:2103688/139‑1 (MQ=255)
                ggCGAAGTGATTGGTTACGCCGTGCGTGCAATCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGTCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGa                                                                                                                            >  1:484800/1‑140 (MQ=255)
                 gggaggtgtTTGTTTCCCCCGGGCGTGCAACCCCCCGCGGAGGCGGGACCGACGAACAAATGGTTTTACACCGGGAGGCGCCGCCGTTACACCCGCTGCCACTGGCAACCAAAGTCCCGGACCCCTTACCGCCGCTGGaa                                                                                                                           <  2:1040898/135‑1 (MQ=255)
                 gagAAGTGATTGTTCCCGCCGTGCGTGCAACCCCCCGCGGAAGCGGGATCGACAAATCAATGTTTGTACTCCCGGAGGCGCCGCCGTTCCACCCGCTCCCACTGGCAACCAAAGTCCCGGACCCCTTACCGCCGCTGGaa                                                                                                                           <  2:568310/138‑1 (MQ=255)
                 gCGAAGTGATTGGTTACGCCGTGCGTGCAATCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGaa                                                                                                                           >  1:573470/1‑140 (MQ=255)
                     aGTGATTGGTTACGCCGTGCGTGCAATCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGaactgt                                                                                                                       >  1:1014764/1‑136 (MQ=255)
                      gTGATTGGTTACGCCGTGCGTGCAATCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATa                                                                                                                      >  2:1902483/1‑140 (MQ=255)
                         gttGGTTCCCCCGTGCGTCAAACCCCCCGCGGAAGCGGGACCGACAAATAAATGGTTGTCCTCCGGGAGGCCCCGCCGTTACCCACGCTGCCACTGGCAACCAAAGCCCCGGACCCCTTACCGCCGCTGGAAGGATacac                                                                                                                   <  2:1842592/139‑1 (MQ=255)
                         aTTGGTTACGCCGTGCGTGCAATCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCAGAACCCTTACCGCCACTGGAAGGATacac                                                                                                                   >  1:612545/1‑140 (MQ=255)
                                ccgcgGTGCGTGCAACCCCCCGGGGAGGCTGGATCGACAAATAAATGTTTTTACTCCCGGAGGCGCCGCCGTTACACACGCTGCCCCTGGCAACCAAAGCCCCGGACCCCTTACCGCCGCTGGAAGGATACACCTTTGAg                                                                                                            <  2:1418371/139‑1 (MQ=255)
                                 cGCCGTGCGTGCAATCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAgg                                                                                                           >  2:974235/1‑140 (MQ=255)
                                   ccGTGCGTGCAATCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGc                                                                                                         >  2:2016085/1‑140 (MQ=255)
                                   ccGTGCGTGCAATCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGc                                                                                                         >  2:1141203/1‑140 (MQ=255)
                                      tgcgtgcAATCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTAt                                                                                                      >  1:1370990/1‑140 (MQ=255)
                                       gcgtgcAATCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAATGATACACCTTTGAGGGCTATc                                                                                                     >  1:1281367/1‑140 (MQ=255)
                                       gcgtgcAATCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATc                                                                                                     >  1:1649227/1‑140 (MQ=255)
                                       gcgtgcAATCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATc                                                                                                     >  2:956359/1‑140 (MQ=255)
                                       gcgtgcAATCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATc                                                                                                     >  1:2039681/1‑140 (MQ=255)
                                        cgtgcAATCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCg                                                                                                    >  1:1016224/1‑140 (MQ=255)
                                        cgtgcAATCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCg                                                                                                    >  1:654456/1‑140 (MQ=255)
                                        cgtgcAATCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCg                                                                                                    >  2:1042662/1‑140 (MQ=255)
                                          ggcAATCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTCCTCCCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGCCCCGGACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCa                                                                                                  <  1:577349/139‑1 (MQ=255)
                                          ggaaaccccACGGGAAAGCTGGACCGACAAATAAATGTTTGTACTCCGGGAGGCCCCGCCGTTACACACGCTGCCACTGGCAACAAAAGTCCCGGACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCa                                                                                                  <  2:1260310/134‑1 (MQ=255)
                                          gccAACCCCCCGGGGAGGCGGGATCGACAAATAAAGGGTTGTCCTCCCGGAGGCCCCGCCGTTACACCCGCTGCTACTGGCAACCAAAGTCCCGGACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCa                                                                                                  <  2:128066/138‑1 (MQ=255)
                                                    cgcgGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATgg                                                                                        >  1:179203/1‑140 (MQ=255)
                                                    cgcgGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATgg                                                                                        >  1:1012384/1‑140 (MQ=255)
                                                    cgcgGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGAAATGCCGATgg                                                                                        >  2:1123329/1‑140 (MQ=255)
                                                    cgcgGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCACCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATgg                                                                                        >  1:1049482/1‑140 (MQ=255)
                                                         aaGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCg                                                                                   >  1:1647719/1‑140 (MQ=255)
                                                         aaGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCg                                                                                   >  1:1602487/1‑140 (MQ=255)
                                                               gATCGATGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCa                                                                             >  1:702735/1‑140 (MQ=255)
                                                                   gacgaATCTATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCaa                                                                         >  1:1767208/1‑140 (MQ=255)
                                                                   gacgaATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCaa                                                                         >  2:527200/1‑140 (MQ=255)
                                                                     cgaATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCaaaa                                                                       >  2:381722/1‑140 (MQ=255)
                                                                       aaTAAATGTTTGTACTCCGGGAGGCCCCGCCGTTACACACGCTGCCACTGGCAACCAAAGCCCCGGACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAAc                                                                     <  2:1286793/140‑1 (MQ=255)
                                                                           aaTGTTTGTACTCCCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACAAAAGTCCCGGACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGc                                                                 <  1:1223728/140‑1 (MQ=255)
                                                                           aaTGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGCCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGc                                                                 <  1:1226073/140‑1 (MQ=255)
                                                                             tgtttGTACTCCGGGAGGCCCCGCCTTTCCCCCCGCCGCCACGGGCAACCAAGGCCCCGGACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTc                                                               <  2:1647719/137‑1 (MQ=255)
                                                                              gtttGTCCTCCGGGACGCGCCGCCGTTACACCCGCTGCCACTGGCAACAAAAGCCCCGGACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCg                                                              <  2:1196235/138‑1 (MQ=255)
                                                                                  gTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTAt                                                          >  2:122605/1‑140 (MQ=255)
                                                                                   tactacCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATc                                                         >  2:2168362/1‑140 (MQ=255)
                                                                                   tactacCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATc                                                         >  2:1601127/1‑140 (MQ=255)
                                                                                     ctccgggAGGCCCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATcac                                                       <  2:968501/135‑1 (MQ=255)
                                                                                       cccGGAGGCCCCGCCGTTACACACGCTGCCACTGGCAACCAAAGCCCCGGACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATcacca                                                     <  2:697361/139‑1 (MQ=255)
                                                                                         cgataGCCCCGCCTTACCAAACGCATCCATTGGCACCAAAGCTCCCGGACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGTTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCaccacc                                                   <  2:1467520/136‑1 (MQ=255)
                                                                                         cGGAAGCCCCGCCTTTACACACGCTGCCACTGGCAACCAAAGTCCCGGACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCaccacc                                                   <  1:753711/140‑1 (MQ=255)
                                                                                          ggAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCAccacca                                                  >  2:1738999/1‑140 (MQ=255)
                                                                                            aggccccgccgTTACACCCGCTGCCACGGGCAACCAAAGCCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGc                                                <  1:899699/135‑1 (MQ=255)
                                                                                               cgccgcCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTc                                             <  1:398011/140‑1 (MQ=255)
                                                                                               ccccgccgTTACCCACGCTGCCACTGGCAACCAAAGCCCCGGACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTc                                             <  2:1649227/138‑1 (MQ=255)
                                                                                               ccccgccgTTACACCCGCTGCCACTGGCAACCAAAGCCCCGGACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTc                                             <  2:591924/138‑1 (MQ=255)
                                                                                                gccgccGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTcc                                            >  1:507222/1‑140 (MQ=255)
                                                                                                gccgccGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTcc                                            >  2:633406/1‑140 (MQ=255)
                                                                                                  cgccgTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCAc                                          >  1:964956/1‑140 (MQ=255)
                                                                                                   gccttaccacccGCGGCCACGGCAAACAAAAGCCCCAAACCCTTTACCGCCACTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACt                                         <  2:612545/133‑1 (MQ=255)
                                                                                                    cggTTACACGCGCTTTCACTGCCAACCAAAGTCCCGGACCCTTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACtg                                        <  2:484800/138‑1 (MQ=255)
                                                                                                          cacacGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCACCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGc                                  <  2:1049482/140‑1 (MQ=255)
                                                                                                            cacGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAg                                <  2:1012384/140‑1 (MQ=255)
                                                                                                            cacGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAg                                <  1:956359/140‑1 (MQ=255)
                                                                                                               gCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCg                             >  2:1037676/1‑140 (MQ=255)
                                                                                                                cTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCgt                            >  1:1169439/1‑140 (MQ=255)
                                                                                                                 tGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCgtg                           <  2:1602487/140‑1 (MQ=255)
                                                                                                                 tGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCgtg                           >  2:447022/1‑140 (MQ=255)
                                                                                                                 ggCCCCTGGCAACAAAAGCCCCGGACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCgtg                           <  2:573470/139‑1 (MQ=255)
                                                                                                                 gcccACTGGCAACCAAAGCCCCGGACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCgtg                           <  2:1747335/138‑1 (MQ=255)
                                                                                                                 cgCCAGTGCCAACCAAGCTCCGGGACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCgtg                           <  2:1767208/139‑1 (MQ=255)
                                                                                                                 cgCCACTGGCAACCAAAGTCCCGGACCCTTTACCGCCGCTGGAAGGATTCACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCgtg                           <  1:974235/139‑1 (MQ=255)
                                                                                                                      ttGGCAGCCAAAGTCCCGGACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGa                      <  2:86980/139‑1 (MQ=255)
                                                                                                                      cTGGCAACCAAAGTCCCGGACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGa                      <  1:1141203/140‑1 (MQ=255)
                                                                                                                      cTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGa                      <  1:1902483/140‑1 (MQ=255)
                                                                                                                      cTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGa                      <  1:1042662/140‑1 (MQ=255)
                                                                                                                      cTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGa                      <  1:381722/140‑1 (MQ=255)
                                                                                                                      cTGGCAACCAAAGCCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGa                      <  1:1654430/140‑1 (MQ=255)
                                                                                                                      cTGGCAACAAAAGTCCCGGACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCATCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGa                      <  1:44686/140‑1 (MQ=255)
                                                                                                                      cTGCCAACCAAAGCCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGa                      <  2:2044061/140‑1 (MQ=255)
                                                                                                                       tGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGAc                     >  2:1387157/1‑140 (MQ=255)
                                                                                                                        ggCAACCAAAGTCCCGGACCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACt                    <  1:51119/140‑1 (MQ=255)
                                                                                                                         gCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGa                      <  2:2168420/137‑1 (MQ=255)
                                                                                                                         gCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGa                      >  1:2168420/1‑137 (MQ=255)
                                                                                                                          cAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTAt                  <  1:527200/140‑1 (MQ=255)
                                                                                                                             ccAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATg                 >  1:149553/1‑138 (MQ=255)
                                                                                                                             ccAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATg                 <  2:149553/138‑1 (MQ=255)
                                                                                                                               aaaGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATgtagt             <  1:2016085/140‑1 (MQ=255)
                                                                                                                                    tCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTaaaa         >  1:255659/1‑140 (MQ=255)
                                                                                                                                    tCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTaaaaa        >  1:1587284/1‑140 (MQ=255)
                                                                                                                                     cccGGACCCCTAACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTaaaaa        <  1:1123329/140‑1 (MQ=255)
                                                                                                                                     cccGGAACCCTTACCGCCGCTGGAATGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTaaaaa        <  2:1281367/140‑1 (MQ=255)
                                                                                                                                        ggAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCa     <  1:122605/140‑1 (MQ=255)
                                                                                                                                        ggAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTAGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCa     >  2:261537/1‑140 (MQ=255)
                                                                                                                                         gAACCCTTACCGCCGCTGGAAGTATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCAt    >  1:67557/1‑140 (MQ=255)
                                                                                                                                         gAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCAt    <  1:1601127/140‑1 (MQ=255)
                                                                                                                                           aCCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTg  <  2:507222/140‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                            
TGAAATTATTCGTTATGGCGAAGTGATTGGTTACGCCGTGCGTGCAATCCCACGCGGAAGCTGGATCGACGAATCAATGGTTGTACTACCGGAAGCGCCGCCGTTACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTACCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTG  >  NC_000913/3275473‑3275750

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