Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,604,3440CT36.1% 45.8 / 13.4 36S62F (TCT→TTT) bglJDNA‑binding transcriptional regulator BglJ
Reads supporting (aligned to +/- strand):  ref base C (6/17);  new base T (13/0);  total (19/17)
Fisher's exact test for biased strand distribution p-value = 1.28e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTG‑CCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCC‑GGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAAC‑GCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACA  >  NC_000913/4604108‑4604586
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                     
cgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCTCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTG‑CCATga                                                                                                                                                                                                                                        <  2:201552/251‑1 (MQ=255)
      gTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTG‑CCATGAGAAGTg                                                                                                                                                                                                                                  <  1:291040/251‑1 (MQ=255)
          aaGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTG‑CCATGAGAAGTGAgcg                                                                                                                                                                                                                              <  1:297199/251‑1 (MQ=255)
                 cttaacttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCtttttttttttttttccctgtgaaggggaggcggag                                                                                                                                                                                                                          >  1:317390/1‑225 (MQ=255)
                   taacttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTG‑CCATGTGAAGTGAGCGCCGAGAggg                                                                                                                                                                                                                     >  2:336571/1‑251 (MQ=255)
                    aacctgcattatGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATTCCTATACCAGTCATGAGTCATTTCAGGATGCCATGTCGCGGAGCTCGTTTGCGGGGGTCATTTTTTCTTTTTCTG‑CCATGAGAAGTGAGCGCAGAGAggg                                                                                                                                                                                                                     <  1:267453/239‑1 (MQ=255)
                             tCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTAGCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTG‑CCATGAGAAGTGAGCGCAGAGAGGGATTATCtt                                                                                                                                                                                                             <  2:309276/249‑1 (MQ=255)
                                gTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAAGGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTG‑CCATGAGAAGTGATCGCAGAGAGGGATTATCTTGc                                                                                                                                                                                                           <  1:203380/248‑1 (MQ=255)
                                    aTGGAGATGACGCCGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGACAGGAGTCATTTCAGGCTGCCATGTCGCGGATCTAGTTTGCGGCGGTCCTGTTTTTTTTTTCTT‑CCATGAGAAGTGAGAGGCGAGAGGGAGTATcattgctgactt                                                                                                                                                                                                    >  2:353531/1‑240 (MQ=255)
                                                    aTGGAACACAACCGCATTAAGAAGAGAAATGTCGAACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGCAACCCTTATAAGCTCCATACCTATCCCAGTCAGGAGTCATTTCAGGATGCCAGGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTCTTTTTTTG‑CCATGGGAAGAGAGAGCCGTGAGGGGTTATTTTGGTTGTCTGGACTGGCggtttaggt                                                                                                                                                                                    >  2:248523/1‑243 (MQ=255)
                                                            caGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATACGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTTCTTTTTTTTTTTTTTT‑CCATGTGAAGTGAGCGCCGAGAGGGATTTTTTTGGCTTACTTAACTGGCGGTTAAGTTTCCGCGt                                                                                                                                                                             >  2:164363/1‑250 (MQ=255)
                                                                      aagaagAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTG‑CCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCC‑GGCGttt                                                                                                                                                                 <  1:454972/251‑1 (MQ=255)
                                                                      aagaagAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTCTT‑CCATGAGGAGAGAGAGCGGAGGGGGATTATCTTGCCTGACTGAACAGGCGATTTAGTTTTCGGGGCCCC‑GGCGttt                                                                                                                                                                 >  2:397002/1‑251 (MQ=255)
                                                                        gaagaGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTG‑CCATGCGAAGTGAGCGCAGAGGGGGATTATCTTGCCTGACTGAACTGTCGGTTAAGTTTCCGCGTACCC‑GGCgatttt                                                                                                                                                               >  1:319536/1‑246 (MQ=255)
                                                                                                                        tATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTT‑CCTTGAGGAGTGGGCCGAGCGAGGGGTTTTTTTGCCCGGATGGACTGGCGGTTAAGTTTCCGCGGACCCCGGGGTTT‑GTTTTTTCGGATGGTGATAATGaaacgcagagcgggagggcatacgcct                                                                                                               >  1:452516/1‑226 (MQ=255)
                                                                                                                        tATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTG‑CCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCC‑GGCGTTTAGTTATTGCGGATGATGATATTGAAGCTc                                                                                                                                    <  2:452516/230‑1 (MQ=255)
                                                                                                                            gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTG‑CCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCC‑GGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATc                                                                                                           <  1:147297/251‑1 (MQ=255)
                                                                                                                               agTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCCCCATTCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCTGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTG‑CCATGAGAAGTGAGCGCAGAGAGGGATTAGCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCC‑GGCGTTTAGTTATTGCGGTTGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCAcc                                                                                                        <  1:566511/251‑1 (MQ=255)
                                                                                                                                gTTTATTCAGAGAGTTTGCGGGTCACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGCTGCCATGTCGCGGATCTCGTTTGCGGCGGTTACTTTTTCTTTTTTTGTCCTTG‑GAAGTGAGCGACGAGGGGGCTTATTTTTCCTTTATGCACTGGCGATTAAGTTTTCCCGGACCC‑GGCGTGTTGTTTTTGTGGATGATGGTAATGAAACGCGGCTGATTGGTTcttttttggcatccacc                                                                                                       >  2:87970/1‑235 (MQ=255)
                                                                                                                                    aTTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTG‑CCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCC‑GGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTgg                                                                                                   <  1:336571/251‑1 (MQ=255)
                                                                                                                                               ttGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTCTG‑CCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTTCCTGACTGTAATGGCGATTAAATTTCCGCGTAACC‑GGCGTTTAGTTATTGCGGGTGATGATATTGAAGCTCGGCTGTTTGGTTTATTTTCTCCATCACCGCTTGGCGGTGTAttt                                                                                        >  1:440705/1‑250 (MQ=255)
                                                                                                                                                    ggTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTG‑CCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCC‑GGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTaa                                                                                   <  2:160863/251‑1 (MQ=255)
                                                                                                                                                       aaCCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTTTTTTTTTTTTTTTT‑CCCTGCGAAGGGAGAGCCGAGAGGGGGTATTTTTCATGGCTGACCTGGCGCATAAGTTTTCGCGTCACC‑GCGGGTTAGTTTTTGCGGAGGGTGTTATTGAAGATCGGCTGATTGGGTCGTTGTCGCCATCCCCGCTGGGGGGGGTATTTAGTaaa                                                                                  >  1:294146/1‑249 (MQ=255)
                                                                                                                                                       aaCCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTG‑CCCTGAGAAGTGAGCGCCGAGGGGGGTTTTCTTGCCCGGCTGGAATGGCGATTAAGTTTTCGCGTACCC‑CGCGTTTTGTTTTTGCGGATGATGATATTTAAGCTCGGCTGGTTGGTTTCTTTTCGCCCTCACCCCTGGGCGGGGTATTTAGTaaaac                                                                                >  1:512899/1‑249 (MQ=255)
                                                                                                                                                         cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTG‑CCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCC‑GGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGt                                                                              <  1:409792/251‑1 (MQ=255)
                                                                                                                                                          ccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTG‑CCAAGAGAAAGGAGCGCAGAGAGGGGTTATATTTCCTGACTGAACTGGGGCTTAAGTTTCCTCGTACCC‑CGCGTTTAGTTTTTGAGGAAGATGATATTTAAGCTCGGATGATTGGTTTCTTGTAGCCATCACCGCTGGACGGGGTATGTAGTAAAACGtt                                                                             >  1:368093/1‑250 (MQ=255)
                                                                                                                                                                     ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTG‑CCCTGAGCAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTTCGCGTACCC‑GGCGTTTTGTTATTGCGGGTGGTTATATTGAAGCTCGGCCGATTGGTTCATTGTCGCCATCACCGCTGGGCGGGGGATTAAGGAAAGCGTCAAC‑GTTgg                                                                    >  1:338317/1‑248 (MQ=255)
                                                                                                                                                                     ccATACCTATACCAGGCATGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTG‑CCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCC‑GGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATGAAGTAAAGCTTCAAC‑GTTGgag                                                                  <  1:443435/250‑1 (MQ=255)
                                                                                                                                                                      cATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTT‑CCCTGAGAAGTGAGCGCAGAGAGGGGTTATCTTTCCTGACTTAACTGGCGATTAAATTTCCGCGTACCC‑CGCGTTTTGTTATTTCGGGTGATGATTTTGAAGCTCGGCTGGTTGGTTTATTTTCGCCCTCCCCGCTGGGCGGGGGATTAAGTAAAGCGGCAAC‑GTTGGGGAt                                                                >  1:318271/1‑251 (MQ=255)
                                                                                                                                                                                        gAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTCTTTTTCTG‑CCATGAGAAGTGAGCGCCGGGAGGGGTTATCTTGCCTGGCTGAACTGGGGCTTTAGTTTCCGCGTACCC‑GGCGTTTTGTTTTTGCGGATGGTGATATTTAAGCTCGGGTGATTGGTTCATTTTCGCCCTCCCCGCTGGACGGTGTATTAAGTAAAACGTCAAC‑GTTGGAGGTTTTTTATCCGGAAATCtt                                              >  1:129571/1‑251 (MQ=255)
                                                                                                                                                                                            cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTT‑GCATGGGAAGGGAGCGCAGAGAGGGGTTATTTTGCCTTACTGAACTGGCGGTTAAGGTTCCGCGTAACC‑GGCGGTTTGTTTTTGCGGGTGATGATATTTAAGCTCGGCTGAGTGGTTCATTTTCGCCATCCCCGCTGGGCGGGGGAATAAATAAAGCGTCAAC‑GGTGGAGATTTTTTATCCGGAAATCtttttt                                          >  2:94062/1‑250 (MQ=255)
                                                                                                                                                                                             aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTCTT‑CCATGGGAAGTGAGCGCCGAGAGGGGTTATCTTTCCTGACCGAACTGGCGATTAAATTTCCGCGGACCC‑GGCGTTTTGTTATTTCGGGTGATGATATTGAAGCTCGGCTGGTTGGTTTCTTGTCGCCATCACCGCTGGGAGGGGTATTTAGTAAAACGTCAACAGTTGG‑GGTTTTTCATAAGGGACTCttttttt                                         >  1:342140/1‑249 (MQ=255)
                                                                                                                                                                                               ttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTG‑CCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCC‑GGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTc                                        <  2:512899/250‑1 (MQ=255)
                                                                                                                                                                                                       gCCATGGCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTG‑CCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTGCCGCGTACCC‑GGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCTCCATCACCGCTGTACGGTGTATTAATTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTGGTCATTAAATg                                <  1:87970/250‑1 (MQ=255)
                                                                                                                                                                                                                            gcggcgGTCATTTTTTCTTTTTCTG‑CCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCC‑GGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGcgaccgacc           <  1:27812/250‑1 (MQ=255)
                                                                                                                                                                                                                              ggcggTCATTTTTTCTTTTTCTG‑CCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCC‑GGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGcgaccgaccgac        >  1:335086/1‑251 (MQ=255)
                                                                                                                                                                                                                                     aTTTTTTCTTTTTCTG‑CCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCC‑GGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGTAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACa  >  1:411592/1‑250 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                     
CGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTG‑CCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCC‑GGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAAC‑GCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACA  >  NC_000913/4604108‑4604586

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