Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,552,1910CA22.0% 194.2 / 19.5 109intergenic (+200/+207)fdnI/yddMformate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (29/56);  new base A (19/5);  total (48/61)
Fisher's exact test for biased strand distribution p-value = 1.27e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 110 total aligned reads displayed.

TATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCAT  >  NC_000913/1552043‑1552340
                                                                                                                                                    |                                                                                                                                                     
tatTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCa                                                                                                                                                      <  2:1829609/150‑1 (MQ=32)
tatTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCa                                                                                                                                                      <  1:27696/150‑1 (MQ=32)
tatTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCa                                                                                                                                                      <  2:118700/150‑1 (MQ=32)
   tGTAGACAACGTAGGCTTTTTTCATGCCGGCTGCGGCGTGTACGCCTTATCCGGCATGAATAGCCTTCATGTCCAATTGCTTGCTCTGAAGGTGTAGGCGTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAAt                                                                                                                                                   <  2:1107445/150‑1 (MQ=11)
   tGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGTCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGTGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAAt                                                                                                                                                   <  1:439828/150‑1 (MQ=21)
   tGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAAt                                                                                                                                                   >  1:808327/1‑150 (MQ=34)
      aGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCTTTTGTCATCAGTTTAAAATGgc                                                                                                                                                >  2:1986716/1‑150 (MQ=17)
      aGACAACGTAGGCTTTGTTCATGCCGGAGGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGgc                                                                                                                                                <  1:930045/150‑1 (MQ=21)
        aCAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGgcgc                                                                                                                                              <  1:1234018/150‑1 (MQ=34)
         catcctAGGCTTTGTTCATGCCGGATGCGGCGTGTACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCt                                                                                                                                             <  1:333755/145‑1 (MQ=19)
            cGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTa                                                                                                                                          <  1:2167474/150‑1 (MQ=37)
            cGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTa                                                                                                                                          <  1:919974/150‑1 (MQ=37)
              tAGGCTTTGTTCATGCCGGATGAGGCGTGAACGGCTTATACTGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTaaa                                                                                                                                        <  1:1407974/150‑1 (MQ=255)
                ggCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAg                                                                                                                                      <  2:1182302/150‑1 (MQ=40)
                ggCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTaaagg                                                                                                                                      >  2:36391/1‑148 (MQ=16)
                 gCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAgg                                                                                                                                     <  1:161772/150‑1 (MQ=40)
                   ttttttCATGCCGTATGCTGCGTTAACGCCTTATACGTCATGAAAACCCTTCAAATCCTATAGATTGCAGTGAACATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGAGTTTGTCTTCAGTTTCAAATGGCGCTGTAAAAGGCg                                                                                                                                   <  1:49950/150‑1 (MQ=14)
                   tttGTTCATGCCTGATGCGGCGTGAACGCCTTATACGTCATGAAAACCCTTCAAATCCAATAGATTGCATTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATATTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCg                                                                                                                                   <  1:46514/150‑1 (MQ=14)
                      gTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGAGTCa                                                                                                                                >  2:752688/1‑150 (MQ=255)
                          aTGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAACTCCAATAGATTGCAGTGAACGTGTAGGCCTGATACGCGTAGCGCATCAGGCAATGTTGCGTTTGTCAGCAGTTTCAAATGGCGCTGTAAAAGGCGTCAtttt                                                                                                                            <  1:288880/150‑1 (MQ=255)
                           tGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTc                                                                                                                           <  1:1566671/150‑1 (MQ=255)
                            gCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTATAAGGCGTCATTTTCa                                                                                                                          >  2:1712196/1‑150 (MQ=255)
                             ccGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCat                                                                                                                         <  1:697755/150‑1 (MQ=255)
                               ggATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTATAAAGAGACATTTTAatat                                                                                                                       >  2:2151617/1‑150 (MQ=11)
                                  tGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGt                                                                                                                    <  1:1680387/150‑1 (MQ=255)
                                    cggcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAg                                                                                                                  <  1:2135803/150‑1 (MQ=255)
                                     ggcgTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGa                                                                                                                 <  1:2133316/150‑1 (MQ=255)
                                      gcgTGAACGCCTTATCCGGCATGAAACCCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGACATTTTCTTATTGTTGAc                                                                                                                >  1:1454140/1‑150 (MQ=255)
                                       cgTGAACGCCTTATCCCGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACa                                                                                                               <  2:2165297/150‑1 (MQ=255)
                                          gAACGCCTTATCCGTCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACg                                                                                                            <  2:937731/150‑1 (MQ=255)
                                             cGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAg                                                                                                         >  1:857015/1‑150 (MQ=255)
                                              gCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAgg                                                                                                        <  2:2039767/150‑1 (MQ=255)
                                              gCCTTATCCGGCATGAAAACCCTTCAAACCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAATGCGTCATTTTGATATTGTAGAACACTTAgg                                                                                                        >  2:1677363/1‑150 (MQ=14)
                                               ccTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTAATATTGTAGACAACGTAGGc                                                                                                       >  2:956057/1‑150 (MQ=255)
                                               ccTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAAGAGGCATTTTATTATTGTATAAAACGTAGGc                                                                                                       >  1:1226235/1‑150 (MQ=11)
                                                cTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCt                                                                                                      <  2:2020933/150‑1 (MQ=255)
                                                cTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTTTAAAAGGCGTCATTTTAATATTGTAGACAAAGTAGGCt                                                                                                      >  1:1930660/1‑150 (MQ=255)
                                                    tCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTATAACGTGACATTTTTATATTGTATACAACGTCGGCttttt                                                                                                  >  2:1241487/1‑148 (MQ=11)
                                                      cGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTAATATTGTAGACAACGTAGTCTTTGTTc                                                                                                >  1:825924/1‑150 (MQ=255)
                                                       ttcATGAAAACCCTTCAAATACAATATATTGCAGTGAAAGTGTAGTCATGATACGCGTAGCGCATCAGTCAATGTTGCGTTTGTCATCAGATTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCa                                                                                               <  1:73945/148‑1 (MQ=21)
                                                       ggCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTTTTCa                                                                                               >  1:1629753/1‑150 (MQ=255)
                                                       ggCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCa                                                                                               >  1:673432/1‑150 (MQ=255)
                                                          aTGAAAAGGCTTGATATCCATTAGATTGCAGTGTACGTGTAGGCCTGAGAAGCGTAGCGATTGAGGCAATGTTTCTTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGc                                                                                            <  1:149301/150‑1 (MQ=25)
                                                           tGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGcc                                                                                           <  1:598850/150‑1 (MQ=255)
                                                              aaaCCCTTCAATTCCTATAGATTGCAGTGCACGTGGAGTCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGa                                                                                        <  2:65143/150‑1 (MQ=21)
                                                              aaaCCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGAATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGa                                                                                        <  2:1346004/150‑1 (MQ=255)
                                                                 cccTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCtt                                                                                                     >  1:44049/1‑134 (MQ=255)
                                                                 cccTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCtt                                                                                                     <  2:44049/134‑1 (MQ=255)
                                                                 cccTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATgt                                                                                     >  1:1852897/1‑149 (MQ=255)
                                                                   cTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTATAAGGGGTCATTTTAATATTGTATACAACGTAGGCTTAGTTCATTCCGGAAgccg                                                                                   >  2:1189934/1‑148 (MQ=17)
                                                                   aatCAAATCCAATAGATTGCAGTGAACGTTTAGGCCTTATACGCGTAGCGCATCATGCAATGTTGCTTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATgcgg                                                                                   <  2:1010905/148‑1 (MQ=25)
                                                                    ttCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCTCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATgcggc                                                                                  <  2:755457/150‑1 (MQ=255)
                                                                     tCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATgcggcg                                                                                 <  1:2241524/150‑1 (MQ=255)
                                                                        aaTCCTTTAGATTGCTGTGTACGTGTAGTCCTGATACGCGGAGCTAATCATGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGGCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGa                                                                              <  1:377606/150‑1 (MQ=21)
                                                                        aaTCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTAAAAGGCGACATTTTCATATTGTAGACACCGAAGGCTTAGTTCATGCCGGATGCGGCGTGa                                                                              >  2:1306641/1‑150 (MQ=21)
                                                                           ccAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTATTATTGTAGACAAAGTAGGCTTTGTTCATGCAGGATGCGGCGTGaaag                                                                           >  1:2205860/1‑148 (MQ=21)
                                                                            cAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAATGGGTCATTTTGATATTGTAGACAAAGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACgg                                                                          >  2:44887/1‑149 (MQ=21)
                                                                                    tGCATTGAAAGTGTAGGCATGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCg                                                                  <  2:526392/150‑1 (MQ=255)
                                                                                     gCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTATAAGGCGACATGTTCATATTGTAGACACCGTAGGCTTAGATCATTCCGAAAGCGTCGTGAAAGACTAATGCgg                                                                 >  2:1129424/1‑150 (MQ=17)
                                                                                      cAGTGAACGTGTAGTCCTGATAAGCGTAGCGCATCATGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                <  1:1969/150‑1 (MQ=255)
                                                                                      cAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                <  2:857015/150‑1 (MQ=255)
                                                                                       aGTGAACGTGTAGGCCTGATAAGCGTATCTCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCa                                                               <  2:938418/150‑1 (MQ=255)
                                                                                       aGTGAACCTTTAGGCCTGATAAGAGTAGCTAATCAGGCAATGTTTCGTTTCTTATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCa                                                               <  1:1382150/150‑1 (MQ=18)
                                                                                          gAACCTCTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGa                                                            <  1:1153377/150‑1 (MQ=255)
                                                                                               tgtTGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTTATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAACAcc                                                       >  2:455473/1‑150 (MQ=255)
                                                                                                gtAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTAAAACGAGGCATTTTAATATTGAAGACAACGTAGTCTTGTTTCATTACAGAGGAGACATGAAAGCCTGATCCGGCATGaaaaaac                                                      >  2:1568909/1‑147 (MQ=12)
                                                                                                 tcgtccTGATAAGCGTAGCGCATCAGGAAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCt                                                     <  2:1949332/146‑1 (MQ=255)
                                                                                                   ggCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTc                                                   <  2:625638/150‑1 (MQ=255)
                                                                                                     ccTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGCGACATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGAGGCGTGAAACACGAATCCGGAAAGAAAACCCATCaa                                                 >  2:242653/1‑150 (MQ=18)
                                                                                                          tAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGAGTCATTTTCATATTGTAGACCACGTAGGCTTTGTTCATGCCGGATGCGTCGTGAAAGGCTTTTACGGAAAGAAAACCCTTAAAATACa                                            >  1:1383602/1‑150 (MQ=18)
                                                                                                            aGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAgg                                                                                                        <  2:1172208/88‑1 (MQ=255)
                                                                                                            aGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAgg                                                                                                        >  1:1172208/1‑88 (MQ=255)
                                                                                                            aGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGAGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAAAGCCTTATCCGGCATGAAAACCCTTAAAATCCAAt                                          >  2:683170/1‑150 (MQ=255)
                                                                                                              ggTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGTCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAg                                        <  2:808327/149‑1 (MQ=255)
                                                                                                              cGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCTTGAACGCCTTTTCCGGCATGAAAACCCTTCAAATCCAATat                                        >  1:1671590/1‑149 (MQ=255)
                                                                                                                   ggcttCAGGTAATTTTTCGTTTTTATTCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGc                                   <  2:2074454/146‑1 (MQ=11)
                                                                                                                   cgaaTCAGGCAATGTTGAGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGc                                   <  2:2080569/147‑1 (MQ=255)
                                                                                                                      aTCAGGCAATTTTTCGTTTTTCATAAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGt                                <  2:1859413/150‑1 (MQ=255)
                                                                                                                        cAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATACGGCATGAAAACCCTTCAAATCCAATAGATTGCAGtta                              >  2:1863253/1‑148 (MQ=255)
                                                                                                                         agtaaaTGTTTAGTTTATTATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGaa                             <  2:1488640/146‑1 (MQ=11)
                                                                                                                         aGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAATGCGTCATTTTCTTATTGTAGAAAACGTAGGATTTGTTCATGCAGTATGCGGCGTGAAAGACTTATACGCCATGAAAACACTTAAAAACCAATAGAATTCAGTGaa                             >  1:352680/1‑150 (MQ=25)
                                                                                                                         aGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGAGTCATTTTAATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATACGGAATGAAAACCCTTCAAATCCAATAGATTGCAGTGaa                             >  2:800368/1‑150 (MQ=17)
                                                                                                                          tttaaTGTTGTTTTTGTGATCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAAc                            <  1:235241/147‑1 (MQ=14)
                                                                                                                          ggCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGACTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGAAGTGAAc                            >  1:1286524/1‑150 (MQ=255)
                                                                                                                          ggCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCAGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATACAATAGATTGCAGTGAAc                            >  1:1892356/1‑150 (MQ=255)
                                                                                                                               tGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGc                                   >  1:1611591/1‑138 (MQ=255)
                                                                                                                               tGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGc                                   <  2:1611591/138‑1 (MQ=255)
                                                                                                                                gTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAg                      >  1:1226128/1‑150 (MQ=255)
                                                                                                                                gTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGAAGCGGCGTGAACGCCTTATCCGGCATGAACACCCTTCAAATCCAATAGAGTGCAGTGAAGGTGTAg                      >  2:1125826/1‑150 (MQ=255)
                                                                                                                                    cGTTTGTCATCAGTTTCAAATGTCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGTGTGAACGCCTTATCCGGCATGAAAACCGTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCt                  <  1:321859/150‑1 (MQ=255)
                                                                                                                                     gTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTg                 >  1:813807/1‑150 (MQ=255)
                                                                                                                                         gTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATaa             >  1:1703920/1‑150 (MQ=255)
                                                                                                                                            atcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCg          >  2:582058/1‑150 (MQ=39)
                                                                                                                                             taaGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCGTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGt         <  2:442137/148‑1 (MQ=12)
                                                                                                                                                gTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCTTGCACGCCTTATCCGGCATGAAACCCCTTCAAATACAATAGATTGCAGTGAACGTGTAGGCGTCATAAGCGTAgc      >  1:1664641/1‑150 (MQ=11)
                                                                                                                                                gTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAgg                     >  1:1331645/1‑135 (MQ=37)
                                                                                                                                                gTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAgg                     <  2:1331645/135‑1 (MQ=37)
                                                                                                                                                gTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgc      <  2:730040/150‑1 (MQ=35)
                                                                                                                                                 tttCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgcg     <  2:2154814/150‑1 (MQ=35)
                                                                                                                                                    cAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCAt  >  2:563690/1‑150 (MQ=32)
                                                                                                                                                    |                                                                                                                                                     
TATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCAT  >  NC_000913/1552043‑1552340

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑