Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,486,3860GA28.3% 79.5 / 11.1 46A89A (GCG→GCAcrpDNA‑binding transcriptional dual regulator CRP
Reads supporting (aligned to +/- strand):  ref base G (26/7);  new base A (0/13);  total (26/20)
Fisher's exact test for biased strand distribution p-value = 7.62e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTC  >  NC_000913/3486137‑3486630
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    
gCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCAACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCg                                                                                                                                                                                                                                                      >  2:492786/1‑250 (MQ=255)
                 tcGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGAAAGATTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTg                                                                                                                                                                                                                                    >  1:624257/1‑251 (MQ=255)
                     aTGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTg                                                                                                                                                                                                                                >  2:571768/1‑251 (MQ=255)
                      tGGTTCTTGCCTCTTTGCCACATTCATAAGTACCCACCCAAGAGCAAGCTTATTCACAAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGTTTTTATTGGCGAACTGGCCCTGTTTGAAGAGGGCCAGGAACGTAGCGCAGGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTg                                                                                                                                                                                                                                <  1:571768/250‑1 (MQ=255)
                          tCTTGTCTCATTGCCACATTCATACCAACCCATCCAAGAGCAAGCTTATTCACAAGGGTGAAAACGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGACGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAAtt                                                                                                                                                                                                                           >  1:559929/1‑251 (MQ=255)
                                                 aaGTCCCCACCCAAGACCAAGCTTATTCCCCAGGGTAAAAAAGGGGAAACGCTTTATTACATCGTTAAAGGTTCGGGGGCAGTGCTAAACAAAGAGGAAGAGGGTAAAGAAAGGACCCCCTCTTATCTAAATAAGGGTTTTTTAATTGGCGAACTGGGCTTGTTTAAAGAGGGCAAGAAACGTAGCGAAGGGGTACGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAAtt                                                                                                                                                                                                    <  2:528696/251‑1 (MQ=255)
                                                     aCCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGAGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGAtt                                                                                                                                                                                                >  1:480604/1‑251 (MQ=255)
                                                         aTCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGCAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAgg                                                                                                                                                                                            >  1:462277/1‑251 (MQ=255)
                                                          tCCAAGAGCAAGTTTATTCACGAGGGAAAAAAAGGGGAAACTCTGTACTACATGTTTAAGGACTTTGTGGCAGTGCTGATAAAAGACGAAGAGGGAAAAGAAAAGTTCCTCCCCTATCGGAATCAGGTTGATTTTATTGGGAAACGGGGCTTGTTGTAGGAGGGCAAGAAACGTAGCGCAGGGTTACGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAg                                                                                                                                                                                             <  2:485754/249‑1 (MQ=255)
                                                              agagCAAGCTTATTAACCAGGGCGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACaa                                                                                                                                                                                                                   >  2:84437/1‑223 (MQ=255)
                                                                       ttATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACAtt                                                                                                                                                                              >  1:69847/1‑251 (MQ=255)
                                                                        tatccACCAGGGTGAAAAAGGGGAACCTCTGTACAACGTTTTTAAAGGCTCTGTGGAGGTGCGGATAAAAGACGAGGGGGGAAAAAAAAGGCCCCGCTCCTATCTAAAGCGGGGTGATTTTATGGGCGAAGTGGGCTTGTTGGAAGGGGGCCAGAAACGTAGCGAAGGGGTAGGTGCAAAACCCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACAt                                                                                                                                                                               <  2:669028/245‑1 (MQ=255)
                                                                             cccAGGGTGAAAAAGCGAAACCGCTGTACTAAATCGTAAAAGGCTCGGTGGCAGTGCTGAACAAAGACGAAGAGGGAAAAGAAATGCCCCTCTCCTATCTGATTCGGGTTGTTTTTATTGGCAACCGGGGCCTGTTTGAAGGGGCCCAGGAACGTAGCGCATGGGTACGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATg                                                                                                                                                                        <  1:368914/250‑1 (MQ=255)
                                                                                    ggAAAAAGCGGAACCGCTGTGCTAAAGCGTAAAAGGCTCTGGGGCAGTGCTGATAAAGAACGAAGGGGTAAAAGAAATGATCCTCTCCTATCAGAATCAGGGTGATTTTATTGGCGAACTGGCCTTGTTGGAAGAGGGCCAGAACCGTAGCGTTGGGTTACGTGTGAAATCCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCTCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGt                                                                                                                                                                 <  1:321388/250‑1 (MQ=255)
                                                                                        aaaGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGc                                                                                                                                                             >  2:203397/1‑251 (MQ=255)
                                                                                          aGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGcac                                                                                                                                                           >  2:74788/1‑251 (MQ=255)
                                                                                                    cTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCATCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGcgcg                                                                                                                                                 >  1:243179/1‑251 (MQ=255)
                                                                                                       tactacACCGTAAAAGGCGCTGTGGCAGTGTTGATAAAAGAGAAAGAGGTTAAAAAAATGACCCTCTCTATGCTGAATCAGGGTGTTTTTATTGGCGAACGGGGCCTGTTGGAAGAGGGCCAGGACGGTAGCGCGGGGGTACGGGAGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGcgtc                                                                                                                                               <  2:127601/250‑1 (MQ=255)
                                                                                                              ccttaaaaGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGTAAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGTTTTTATTGGCGAACGGGGCCTGTTGAAAGAGGGCCAGGAACGTAGCGCATGGGTAGGTGGAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCaa                                                                                                                                       <  1:74788/246‑1 (MQ=255)
                                                                                                                    aaGGCTCTGTGGCAGTGCGGATCAAAGACGAAGAGGGGAAAAAAATGTTCCTCTCCTATCTGAATCAGGTTGATTTTATGGGCAAACGGGGCCTTTTTGAAGAGGGCCAGGAACGTGGCGCATGGGTCCGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACt                                                                                                                                 <  1:224029/251‑1 (MQ=255)
                                                                                                                           atgtgGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCag                                                                                                                             >  2:520500/2‑249 (MQ=255)
                                                                                                                                gCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGaaatt                                                                                                                      >  1:164726/1‑249 (MQ=255)
                                                                                                                                gCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCACCTCTCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGt                                                                                                                      >  1:599144/1‑251 (MQ=255)
                                                                                                                                 cAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATCCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTg                                                                                                                     >  1:211658/1‑251 (MQ=255)
                                                                                                                                     gcggATCAAAGACGAAGAGGGAAAAAAAATGCCCCCCCCCTATCTGAATCAGGGTGTTTTTATGGGCGACGGGGGCCTGTTGGAAGGGGGCCAGGAACGTAGCGAAGGGGTACGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGc                                                                                                                  <  1:8043/247‑1 (MQ=255)
                                                                                                                                     gCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCTTGTTGGAAGAGGGCCAGGAACGTAGCGCAGGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCa                                                                                                                 <  1:50994/251‑1 (MQ=255)
                                                                                                                                        gATCAAAGACGCAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCCGAGAAAGTGGGCAAcc                                                                                                              >  1:241807/1‑251 (MQ=255)
                                                                                                                                         aTCAAAGACGAAGAGGGTAAAGAACTGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGACGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCt                                                                                                             >  2:575692/1‑251 (MQ=255)
                                                                                                                                           cAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTgg                                                                                                           >  1:554240/1‑251 (MQ=255)
                                                                                                                                                     agagGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGCAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGCTGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCg                                                                                                  >  2:577019/1‑251 (MQ=255)
                                                                                                                                                     agagGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCg                                                                                                  >  1:479092/1‑251 (MQ=255)
                                                                                                                                                            aaaGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGAc                                                                                           >  1:509790/1‑251 (MQ=255)
                                                                                                                                                              aGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACg                                                                                          >  1:232397/1‑250 (MQ=255)
                                                                                                                                                              aGAAATGATCCCCCCCTAACTGAATCGGGGTGTTTTTATTGGCGAACTGGCCCTGTTTGAAGGGGGCCAGGAACGAGGCGCAGGGTTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACg                                                                                          <  2:232397/250‑1 (MQ=255)
                                                                                                                                                                           ctcCTATCTGAATAGGGGTGTTTGTATTGGCGAACGGGGTCTTATTGAAGAGGCCAAGGAACGTAGGCTTTGGTTAGTTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATTCGTTTGTCTGCACAGATGGCGCGTCGTCGGCAAGTCTCTGCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGGATTGc                                                                               <  2:55196/249‑1 (MQ=255)
                                                                                                                                                                                              tgtttttATGGGCGAACGGGGCTTGTTTGAAGAGGGCCAGAAACGTCGCGCATGGGTACGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTg                                                            <  1:203397/247‑1 (MQ=255)
                                                                                                                                                                                                           gAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGCGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTAGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAgg                                              >  2:105253/1‑250 (MQ=255)
                                                                                                                                                                                                                      gttagAAGAGGGCAGGGACCGTAGTGCATGGGTGTGTGAAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGGTTGGCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTTCACAGACGCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTc                                   <  1:319877/247‑1 (MQ=255)
                                                                                                                                                                                                                      gTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGt                                                                                                                      >  1:484462/1‑168 (MQ=255)
                                                                                                                                                                                                                       tttGAAGAGGGCCAGGAACATAGCGCATGGGAACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTa                                         >  2:406426/1‑244 (MQ=255)
                                                                                                                                                                                                                       ttggAAGGGGGCCAGGACCTTAGCGAAGGGTTACGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTa                                         <  1:406426/241‑1 (MQ=255)
                                                                                                                                                                                                                         tGAAGAGGGCCAGGAACGTAGCGCAGGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCAc                                 <  1:496059/250‑1 (MQ=255)
                                                                                                                                                                                                                              ggggCCAGAAAGGTAGGGGATGGGTACGGGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGAc                           <  1:63216/250‑1 (MQ=255)
                                                                                                                                                                                                                                          gtggcgcATGGGGACGAGGAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAAGTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCATACGCTTTGACTCATCCGGACGGTATTCAAATc               <  2:559929/248‑1 (MQ=255)
                                                                                                                                                                                                                                           tAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATc               >  1:467089/1‑250 (MQ=255)
                                                                                                                                                                                                                                              cgcATGGGTACGTACAAAAGCCGCCTGTGAAGTGGCTGAAAGTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGATAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAAGCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCaaaa           <  1:490106/251‑1 (MQ=255)
                                                                                                                                                                                                                                                           gTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTc  >  1:481394/1‑249 (MQ=255)
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    
GCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTC  >  NC_000913/3486137‑3486630

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