Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,486,3860GA21.1% 98.9 / 10.8 52A89A (GCG→GCAcrpDNA‑binding transcriptional dual regulator CRP
Reads supporting (aligned to +/- strand):  ref base G (31/10);  new base A (0/11);  total (31/21)
Fisher's exact test for biased strand distribution p-value = 5.84e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGA  >  NC_000913/3486144‑3486634
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                        
gACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGc                                                                                                                                                                                                                                                  <  2:226252/251‑1 (MQ=255)
         ctcGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGt                                                                                                                                                                                                                                         >  1:222673/1‑251 (MQ=255)
           cGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTgg                                                                                                                                                                                                                                       >  2:509449/1‑251 (MQ=255)
               tggttatcgcctacctGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTg                                                                                                                                                                                                                                    >  1:455548/16‑250 (MQ=255)
                 gTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGaa                                                                                                                                                                                                                                  >  1:319078/1‑250 (MQ=255)
                      tttttaaTTGCCACATTCATAAGTCCCCACCCAAGAGCCCGTTTATCCTCCAGGGTAAAAAGGCGAAAACGCTGTAGTACACCTTTAAAGTCTCTGGGGCTGTGCTGATCAAAGACGAAGAGTTTAAAGAAAGTACCCTCTCCTATCTGAATCAGGGTGTTGTTTTGGGCAATCGGGGCCGGTTTGAAGAGGGCCAGGACGGGACCGCAGGGGAACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTc                                                                                                                                                                                                                             <  1:351325/244‑1 (MQ=255)
                                ccACATTCATAAGACCCCACCAAAGATCAGGCTTATCCACAAGGGTAAAAAAGCGAAACCTCGGTACAACATCTTAAAAGGCTTTGTGGCAGTGCTGATCAAAGCCAAAGAGGGAAAAGAAATGATCCTCTCCTATCGGAACCGGGTTGTTTTTTTGGTCGATCGGGGCTTTTTTGAGGAGGGCCAGGAACGTAGCGCATGGGTACGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAAt                                                                                                                                                                                                                  <  1:275635/251‑1 (MQ=255)
                                             cccccATCAAAGACCACGTTTATTCACCGGGGTGAAAAAGCGGAAACGCTGTACAACATCTATAAAGCCTTTGTGGCTGTGCTGATAAAAGAGGAAGGGGGTAAAAAAATGATCCTCTCCTTTCTGAATCAGGGTGTTTTTATGGTCGACCTGGGCCTGTTTGAAGAGGGCCAGGACCGTACCGCATGGGTCCGTGAAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGAt                                                                                                                                                                                                     <  1:8227/249‑1 (MQ=255)
                                               cccATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGACGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAACTTTCGTACAACAAATTTCGCCAATTGATTc                                                                                                                                                                                                   >  1:365929/1‑251 (MQ=255)
                                                        gagCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGGGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAACTTTCGTACAAAAACTTTCGACACTTGATTCAGGTAAAcc                                                                                                                                                                                          >  2:17907/1‑251 (MQ=255)
                                                        aagCACGCTTATTCCCCAGGGTGAAAAAGGGAAAACGGTGTCCGACATCGTTAAAGGCCCGGTGGCAGTGCGGACCAAAGACGAAGAGTGTAAAGAAATGACCCTCCCCAGTCTAAATCGGGGTGATTTTATTGGCGAACTGGCCTTGTTTGAAGAGGGCCAGGAAGGTAGCGCATGGGTTCTTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAAcc                                                                                                                                                                                          <  1:56926/250‑1 (MQ=255)
                                                                ttATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACAtt                                                                                                                                                                                  >  1:260522/1‑251 (MQ=255)
                                                                        cAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACACGGACATTCTGATGCg                                                                                                                                                                          >  1:459830/1‑251 (MQ=255)
                                                                              aaaaaaGCGGAAACGCTGTCCTACATCGTTAAAGGCTCTGTGGAAGTGCTGACAAAAGACGAAGAGGGAAAAGAAATGACCCCCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACGGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTc                                                                                                                                                                    <  1:440494/250‑1 (MQ=255)
                                                                                  aaggggAAACGTTGTCCTACATCGTTAAAGGCGCTGGGGCAGTGCTGAACAAAGACGAAGAGGGTAAAGAAATGATCCTCCCCTATCTGAATCAGGGTGTTTTTATTGGCGAACTGGGCTTGTTTAAAGAGGGCCAGAACCGTAGCGCAGGTGTCCGGGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGca                                                                                                                                                                <  1:530297/247‑1 (MQ=255)
                                                                                      ggAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGAGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCAAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGcac                                                                                                                                                               >  2:79867/1‑248 (MQ=255)
                                                                                            gCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCCGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGgcgc                                                                                                                                                      >  1:34356/1‑251 (MQ=255)
                                                                                                    acATCGATAAAGGCTCTTTGGCCGCGCTCATCTACGACGAAGAGGGTAAAGAAATGATCCTATCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCAAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGGAAACCCGGACATTCTGATGCGTTTGTCAGCACAGCTGGCGCGTCGTCTg                                                                                                                                              >  2:112588/1‑251 (MQ=255)
                                                                                                                       ggCAGTGCTGATCAAAGACGCAGAGGGTAAAGAAATGATCCTCTCCTATCTGCATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGACTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGaa                                                                                                                             >  1:216552/1‑249 (MQ=255)
                                                                                                                             gCTGATAAAAGACGAAGAGGATAAAGAAATGATCCTCTCCTTTCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACTTAGCGCTTGGTTACTTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCa                                                                                                                     <  1:87191/251‑1 (MQ=255)
                                                                                                                                  acAAAGACGAAGAGGAAAAAAAAATGATCCTCTCTTAACTGAATAAGGGTATTTTTATTGGATAATGGGGCCTGTGTGAAGAGGGCCAGGAACGTAGCGAAGGGGTCCGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTg                                                                                                                <  2:335825/250‑1 (MQ=255)
                                                                                                                                    aaaGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCt                                                                                                                                               >  1:297227/1‑218 (MQ=255)
                                                                                                                                     aaGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTAAGAGAAAGTGGGCAACCTGGCg                                                                                                             >  1:481676/1‑251 (MQ=255)
                                                                                                                                         ggAAGAGGGTAAAGAAATGACCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTcc                                                                                                         <  2:181411/250‑1 (MQ=255)
                                                                                                                                          gAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGCGAAAGTGGGCAACCTGGCGTTCCt                                                                                                        >  2:517624/1‑251 (MQ=255)
                                                                                                                                           aaGAGGTTAAAGAAATGCTCCTCCCCTATTTGAATCCGGGGGATTTTATGGGCGACCTGGGCCGTTTTGAAGGGGCCAAGGAACGTAGCGCATGGGTCCGCGCAAATACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCGGCACAGATGGCGCGTCGGCTGCAAGTCACTTCAGATAAAGTGGTCAACCTGGCTTTCCTc                                                                                                       <  2:401096/251‑1 (MQ=255)
                                                                                                                                              agGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGAc                                                                                                    >  1:21575/1‑251 (MQ=255)
                                                                                                                                                 gTAAAGAAATGATACTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGCACGTCGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGCACAAAAAATTTAGCCAATTGATTCCGGTAAACCCGGCCTTTCTGATGCGTTTGTCTGCCCCGATGGCGCGTCGTCTGCAAGTCCCTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGAAGTg                                                                                                 >  2:574648/1‑251 (MQ=255)
                                                                                                                                                      gAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACggg                                                                                            >  1:226051/1‑251 (MQ=255)
                                                                                                                                                        aaTGATCCTCTCCTATCTGAATCAGGGTGCTTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGcc                                                                                          >  1:494699/1‑251 (MQ=255)
                                                                                                                                                        aaTGATCCTCTCCTATCGGAATCAGGGTGTTTTTATTGGAGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGGGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGcc                                                                                          <  2:319078/251‑1 (MQ=255)
                                                                                                                                                                tctcCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCACTTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGc                                                                                   >  1:507603/1‑250 (MQ=255)
                                                                                                                                                                      aTCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAgcct                                                                            >  2:256100/1‑248 (MQ=255)
                                                                                                                                                                       tCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGActc                                                                           >  1:271739/1‑251 (MQ=255)
                                                                                                                                                                        cggaataaggggtttttgaTTGGCGAACGGGGCCTGTTTGAAGAGGGCAAGGAACGTAGCGCAGGGGTACGGGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGActct                                                                          <  1:509449/233‑1 (MQ=255)
                                                                                                                                                                              cAGGGTGATTTTATTGGCGACCTGGCCCTGTTGGAAGAGGGCCAGGAACGTGGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGa                                                                     <  1:435782/250‑1 (MQ=255)
                                                                                                                                                                                        ttgagtgggaaacggggCCTGTTTGAAGAGGGCCAGAAACTTAGCGAATGGTTACGGCCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAAGTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCaa                                                            <  1:256970/236‑1 (MQ=255)
                                                                                                                                                                                          aTTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACaaaaa                                                                                                                                                                                                                    >  1:37978/1‑95 (MQ=255)
                                                                                                                                                                                          aTTGGCGAACGGGGCTTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACTTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACaaaaa                                                                                                                                                                                                                    <  2:37978/95‑1 (MQ=255)
                                                                                                                                                                                                       ggCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGCAAACCGCCTGTGACGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTAt                                            >  2:135510/1‑251 (MQ=255)
                                                                                                                                                                                                        ccgttttggAAGGGGGCGAGAAACGTAGCGCATGGGTATGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATg                                           <  2:222673/243‑1 (MQ=255)
                                                                                                                                                                                                               tGAAGGGGTCCAGGAACGTGCGGCAGGGGTCCGTGGAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCAc                                     <  1:126027/250‑1 (MQ=255)
                                                                                                                                                                                                               tGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCAcc                                    >  1:344819/1‑251 (MQ=255)
                                                                                                                                                                                                                gaagacgcccAGGAACGTAGCGCATGGGCACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAACCAACCAGCCGCTATGCCTCAccc                                   >  2:293918/9‑251 (MQ=255)
                                                                                                                                                                                                                 aaGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCg                                  >  1:436412/1‑251 (MQ=255)
                                                                                                                                                                                                                    agGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGAc                               >  1:505187/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gatagatcagggggactggCGAAAACCGACTGTGAAGTGGCTGAAATTTCGTACAAAAAAGTTCGCCAAGTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGGCACTTCAGAGAAAGTGGGCACCCTGGCGTTCCTCGACGGGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCCCCCGGACGGTATGCAAAt                    <  1:546457/233‑1 (MQ=255)
                                                                                                                                                                                                                                   aGCGCATGGGTACGGGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAAAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACACCCAAACGCTATGACTAACCCGGACGGTTTGCAAATCa                  >  2:376789/1‑250 (MQ=255)
                                                                                                                                                                                                                                     cgcATGGGTACGCGCGAAAACCGCCTGTGAAGTGGCTGAAATTTAGTACAAAAAATTTCGCCAATTGATTAAGGGAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTAGGACTCACCCGGACGGGATGCAAATCaaa                >  2:543696/1‑250 (MQ=255)
                                                                                                                                                                                                                                            ggTACGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTAccc         <  1:321277/251‑1 (MQ=255)
                                                                                                                                                                                                                                                cGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCa     <  1:104025/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    cGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAACTTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGa  >  2:391418/1‑250 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                        
GACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGA  >  NC_000913/3486144‑3486634

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: