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→GCA) crpDNA‑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.

GACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGG‑CGAACTGGGCCTGTTTGAAGAGGGCCA‑GGAACGTAGCG‑CATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGA  >  NC_000913/3486144‑3486634
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                        
gACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGG‑CGAACTGGGCCTGTTTGAAGAGGGCCA‑GGAACGTAGCG‑CATGGGTACGTGCGAAAACCGc                                                                                                                                                                                                                                                  <  2:226252/251‑1 (MQ=255)
         ctcGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGG‑CGAACTGGGCCTGTTTGAAGAGGGCCA‑GGAACGTAGCG‑CATGGGTACGTGCGAAAACCGCCTGTGAAGt                                                                                                                                                                                                                                         >  1:222673/1‑251 (MQ=255)
           cGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGG‑CGAACTGGGCCTGTTTGAAGAGGGCCA‑GGAACGTAGCG‑CATGGGTACGTGCGAAAACCGCCTGTGAAGTgg                                                                                                                                                                                                                                       >  2:509449/1‑251 (MQ=255)
               tggttatcgcctacctGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGG‑CGAACTGGGCCTGTTTGAAGAGGGCCA‑GGAACGTAGCG‑CATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTg                                                                                                                                                                                                                                    >  1:455548/16‑250 (MQ=255)
                 gTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGG‑CGAACTGGGCCTGTTTGAAGAGGGCCA‑GGAACGTAGCG‑CATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGaa                                                                                                                                                                                                                                  >  1:319078/1‑250 (MQ=255)
                      tttttaaTTGCCACATTCATAAGTCCCCACCCAAGAGCCCGTTTATCCTCCAGGGTAAAAAGGCGAAAACGCTGTAGTACACCTTTAAAGTCTCTGGGGCTGTGCTGATCAAAGACGAAGAGTTTAAAGAAAGTACCCTCTCCTATCTGAATCAGGGTGTTGTTTTGGG‑CAATCGGGGCCGGTTTGAAGAGGGCCA‑GGACGGGACCG‑CAGGGGAACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTc                                                                                                                                                                                                                             <  1:351325/244‑1 (MQ=255)
                                ccACATTCATAAGACCCCACCAAAGATCAGGCTTATCCACAAGGGTAAAAAAGCGAAACCTCGGTACAACATCTTAAAAGGCTTTGTGGCAGTGCTGATCAAAGCCAAAGAGGGAAAAGAAATGATCCTCTCCTATCGGAA‑CCGGGTTGTTTTTTTGGTCGATCGGGGCTTTTTTGAGGAGGGCCA‑GGAACGTAGCG‑CATGGGTACGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAAt                                                                                                                                                                                                                  <  1:275635/251‑1 (MQ=255)
                                             cccccATCAAAGACCACGTTTATTCACCGGGGTGAAAAAGCGGAAACGCTGTACAACATCTATAAAGCCTTTGTGGCTGTGCTGATAAAAGAGGAAGGGGGTAAAAAAATGATCCTCTCCTTTCTGAATCAGGGTGTTTTTATGGT‑CGACCTGGGCCTGTTTGAAGAGGGCCA‑GGACCGTACCG‑CATGGGTCCGTGAAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGAt                                                                                                                                                                                                     <  1:8227/249‑1 (MQ=255)
                                               cccATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGG‑CGAACTGGGCCTGTTTGACGAGGGCCA‑GGAACGTAGCG‑CATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAACTTTCGTACAACAAATTTCGCCAATTGATTc                                                                                                                                                                                                   >  1:365929/1‑251 (MQ=255)
                                                        gagCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGGGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGG‑CGAACTGGGCCTGTTTGAAGAGGGCCA‑GGAACGTAGCG‑CATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAACTTTCGTACAAAAACTTTCGACACTTGATTCAGGTAAAcc                                                                                                                                                                                          >  2:17907/1‑251 (MQ=255)
                                                        aagCACGCTTATTCCCCAGGGTGAAAAAGGGAAAACGGTGTCCGACATCGTTAAAGGCCCGGTGGCAGTGCGGACCAAAGACGAAGAGTGTAAAGAAATGACCCTCCCCAGTCTAAATCGGGGTGATTTTATTGG‑CGAACTGGCCTTGTTTGAAGAGGGCCA‑GGAAGGTAGCG‑CATGGGTTCTTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAAcc                                                                                                                                                                                          <  1:56926/250‑1 (MQ=255)
                                                                ttATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGG‑CGAACTGGGCCTGTTTGAAGAGGGCCA‑GGAACGTAGCG‑CATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACAtt                                                                                                                                                                                  >  1:260522/1‑251 (MQ=255)
                                                                        cAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGG‑CGAACTGGGCCTGTTTGAAGAGGGCCA‑GGAACGTAGCG‑CATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACACGGACATTCTGATGCg                                                                                                                                                                          >  1:459830/1‑251 (MQ=255)
                                                                              aaaaaaGCGGAAACGCTGTCCTACATCGTTAAAGGCTCTGTGGAAGTGCTGACAAAAGACGAAGAGGGAAAAGAAATGACCCCCTCCTATCTGAATCAGGGTGATTTTATTGG‑CGAACGGGGCCTGTTTGAAGAGGGCCA‑GGAACGTAGCG‑CATGGGTACGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTc                                                                                                                                                                    <  1:440494/250‑1 (MQ=255)
                                                                                  aaggggAAACGTTGTCCTACATCGTTAAAGGCGCTGGGGCAGTGCTGAACAAAGACGAAGAGGGTAAAGAAATGATCCTCCCCTATCTGAATCAGGGTGTTTTTATTGG‑CGAACTGGGCTTGTTTAAAGAGGGCCA‑GAACCGTAGCG‑CAGGTGTCCGGGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGca                                                                                                                                                                <  1:530297/247‑1 (MQ=255)
                                                                                      ggAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGAGATTTTATTGG‑CGAACTGGGCCTGTTTGAAGAGGGCAA‑GGAACGTAGCG‑CATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGcac                                                                                                                                                               >  2:79867/1‑248 (MQ=255)
                                                                                            gCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCCGGGTGATTTTATTGG‑CGAACTGGGCCTGTTTGAAGAGGGCCA‑GGAACGTAGCG‑CATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGgcgc                                                                                                                                                      >  1:34356/1‑251 (MQ=255)
                                                                                                    acATCGATAAAGGCTCTTTGGCCGCGCTCATCTACGACGAAGAGGGTAAAGAAATGATCCTATCCTATCTGAATCAGGGTGATTTTATTGG‑CGAACTGGGCCTGTTTGAAGAGGGCAA‑GGAACGTAGCG‑CATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGGAAACCCGGACATTCTGATGCGTTTGTCAGCACAGCTGGCGCGTCGTCTg                                                                                                                                              >  2:112588/1‑251 (MQ=255)
                                                                                                                       ggCAGTGCTGATCAAAGACGCAGAGGGTAAAGAAATGATCCTCTCCTATCTGCATCAGGGTGATTTTATTGG‑CGAACTGGGCCTGTTTGAAGAGGGCCA‑GGAACGTAGCG‑CATGGGTACGTGCGAAAACCGACTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGaa                                                                                                                             >  1:216552/1‑249 (MQ=255)
                                                                                                                             gCTGATAAAAGACGAAGAGGATAAAGAAATGATCCTCTCCTTTCTGAATCAGGGTGATTTTATTGG‑CGAACTGGGCCTGTTTGAAGAGGGCCA‑GGAACTTAGCG‑CTTGGTTACTTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCa                                                                                                                     <  1:87191/251‑1 (MQ=255)
                                                                                                                                  acAAAGACGAAGAGGAAAAAAAAATGATCCTCTCTTAACTGAATAAGGGTATTTTTATTGG‑ATAATGGGGCCTGTGTGAAGAGGGCCA‑GGAACGTAGCG‑AAGGGGTCCGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTg                                                                                                                <  2:335825/250‑1 (MQ=255)
                                                                                                                                    aaaGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGG‑CGAACTGGGCCTGTTTGAAGAGGGCCA‑GGAACGTAGCG‑CATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCt                                                                                                                                               >  1:297227/1‑218 (MQ=255)
                                                                                                                                     aaGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGG‑CGAACTGGGCCTGTTTGAAGAGGGCCA‑GGAACGTAGCG‑CATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTAAGAGAAAGTGGGCAACCTGGCg                                                                                                             >  1:481676/1‑251 (MQ=255)
                                                                                                                                         ggAAGAGGGTAAAGAAATGACCCTCTCCTATCTGAATCAGGGTGATTTTATTGG‑CGAACTGGGCCTGTTTGAAGAGGGCCA‑GGAACGTAGCG‑CATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTcc                                                                                                         <  2:181411/250‑1 (MQ=255)
                                                                                                                                          gAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGG‑CGAACTGGGCCTGTTTGAAGAGGGCCA‑GGAACGTAGCG‑CATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGCGAAAGTGGGCAACCTGGCGTTCCt                                                                                                        >  2:517624/1‑251 (MQ=255)
                                                                                                                                           aaGAGGTTAAAGAAATGCTCCTCCCCTATTTGAATCCGGGGGATTTTATGGG‑CGACCTGGGCCGTTTTGAAG‑GGGCCAAGGAACGTAGCG‑CATGGGTCCGCGCAAATACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCGGCACAGATGGCGCGTCGGCTGCAAGTCACTTCAGATAAAGTGGTCAACCTGGCTTTCCTc                                                                                                       <  2:401096/251‑1 (MQ=255)
                                                                                                                                              agGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGG‑CGAACTGGGCCTGTTTGAAGAGGGCCA‑GGAACGTAGCG‑CATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGAc                                                                                                    >  1:21575/1‑251 (MQ=255)
                                                                                                                                                 gTAAAGAAATGATACTCTCCTATCTGAATCAGGGTGATTTTATTGG‑CGAACTGGGCCTGTTTGAAGAGGGCCA‑GGCACGTCGCG‑CATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGCACAAAAAATTTAGCCAATTGATTCCGGTAAACCCGGCCTTTCTGATGCGTTTGTCTGCCCCGATGGCGCGTCGTCTGCAAGTCCCTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGAAGTg                                                                                                 >  2:574648/1‑251 (MQ=255)
                                                                                                                                                      gAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGG‑CGAACTGGGCCTGTTTGAAGAGGGCCA‑GGAACGTAGCG‑CATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACggg                                                                                            >  1:226051/1‑251 (MQ=255)
                                                                                                                                                        aaTGATCCTCTCCTATCTGAATCAGGGTGCTTTTATTGG‑CGAACTGGGCCTGTTTGAAGAGGGCCA‑GGAACGTAGCG‑CATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGcc                                                                                          >  1:494699/1‑251 (MQ=255)
                                                                                                                                                        aaTGATCCTCTCCTATCGGAATCAGGGTGTTTTTATTGG‑AGAACTGGGCCTGTTTGAAGAGGGCCA‑GGAACGTAGCG‑CATGGGTACGGGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGcc                                                                                          <  2:319078/251‑1 (MQ=255)
                                                                                                                                                                tctcCTATCTGAATCAGGGTGATTTTATTGG‑CGAACTGGGCCTGTTTGAAGAGGGCCA‑GGAACGTAGCG‑CATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCACTTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGc                                                                                   >  1:507603/1‑250 (MQ=255)
                                                                                                                                                                      aTCTGAATCAGGGTGATTTTATTGG‑CGAACTGGGCCTGTTTGAAGAGGGCCA‑GGAACGTAGCG‑CATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAgcct                                                                            >  2:256100/1‑248 (MQ=255)
                                                                                                                                                                       tCTGAATCAGGGTGATTTTATTGG‑CGAACTGGGCCTGTTTGAAGAGGGCCA‑GGAACGTAGCG‑CATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGActc                                                                           >  1:271739/1‑251 (MQ=255)
                                                                                                                                                                        cggaataaggggtttttgaTTGG‑CGAACGGGGCCTGTTTGAAGAGGGCAA‑GGAACGTAGCG‑CAGGGGTACGGGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGActct                                                                          <  1:509449/233‑1 (MQ=255)
                                                                                                                                                                              cAGGGTGATTTTATTGG‑CGACCTGGCCCTGTTGGAAGAGGGCCA‑GGAACGTGGCG‑CATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGa                                                                     <  1:435782/250‑1 (MQ=255)
                                                                                                                                                                                        ttgagtgggaaacggggCCTGTTTGAAGAGGGCCA‑GAAACTTAGCG‑AATGGTTACGGCCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAAGTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCaa                                                            <  1:256970/236‑1 (MQ=255)
                                                                                                                                                                                          aTTGG‑CGAACTGGGCCTGTTTGAAGAGGGCCA‑GGAACGTAGCG‑CATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACaaaaa                                                                                                                                                                                                                    >  1:37978/1‑95 (MQ=255)
                                                                                                                                                                                          aTTGG‑CGAACGGGGCTTGTTTGAAGAGGGCCA‑GGAACGTAGCG‑CATGGGTACTTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACaaaaa                                                                                                                                                                                                                    <  2:37978/95‑1 (MQ=255)
                                                                                                                                                                                                       ggCCTGTTTGAAGAGGGCCA‑GGAACGTAGCG‑CATGGGTACGTGCGCAAACCGCCTGTGACGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTAt                                            >  2:135510/1‑251 (MQ=255)
                                                                                                                                                                                                        ccgttttggAAGGGGGCGA‑GAAACGTAGCG‑CATGGGTATGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATg                                           <  2:222673/243‑1 (MQ=255)
                                                                                                                                                                                                               tGAAGGGGTCCA‑GGAACGT‑GCGGCAGGGGTCCGTGGAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCAc                                     <  1:126027/250‑1 (MQ=255)
                                                                                                                                                                                                               tGAAGAGGGCCA‑GGAACGTAGCG‑CATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCAcc                                    >  1:344819/1‑251 (MQ=255)
                                                                                                                                                                                                                gaagacgcccA‑GGAACGTAGCG‑CATGGGCACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAACCAACCAGCCGCTATGCCTCAccc                                   >  2:293918/9‑251 (MQ=255)
                                                                                                                                                                                                                 aaGAGGGCCA‑GGAACGTAGCG‑CATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCg                                  >  1:436412/1‑251 (MQ=255)
                                                                                                                                                                                                                    agGGCCA‑GGAACGTAGCG‑CATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGAc                               >  1:505187/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gatagatcagggggactggCGAAAACCGACTGTGAAGTGGCTGAAATTTCGTACAAAAAAGTTCGCCAAGTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGGCACTTCAGAGAAAGTGGGCACCCTGGCGTTCCTCGACGGGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCCCCCGGACGGTATGCAAAt                    <  1:546457/233‑1 (MQ=255)
                                                                                                                                                                                                                                   aGCG‑CATGGGTACGGGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAAAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACACCCAAACGCTATGACTAACCCGGACGGTTTGCAAATCa                  >  2:376789/1‑250 (MQ=255)
                                                                                                                                                                                                                                     cg‑cATGGGTACGCGCGAAAACCGCCTGTGAAGTGGCTGAAATTTAGTACAAAAAATTTCGCCAATTGATTAAGGGAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTAGGACTCACCCGGACGGGATGCAAATCaaa                >  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)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                        
GACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGG‑CGAACTGGGCCTGTTTGAAGAGGGCCA‑GGAACGTAGCG‑CATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGA  >  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: ‑