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

GGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTAT‑CTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATG‑ATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACA  >  NC_000913/4604096‑4604586
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
ggATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttct                                                                                                                                                                                                                                                       <  1:494303/250‑1 (MQ=255)
  aTAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttt                                                                                                                                                                                                                                                    <  1:310548/251‑1 (MQ=255)
    aaaTAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATTCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGTTCATtttttctttttc                                                                                                                                                                                                                                                  <  2:89078/251‑1 (MQ=255)
    aaaTAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATTCCATGTCGCGGATCTCGTTTGCGGCGGTCATttttacgttttc                                                                                                                                                                                                                                                  <  1:299646/251‑1 (MQ=255)
      aTAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTttttctttttctt                                                                                                                                                                                                                                                >  1:392024/1‑250 (MQ=255)
         gcgcgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCa                                                                                                                                                                                                                                             <  1:44576/251‑1 (MQ=255)
              cAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATgaga                                                                                                                                                                                                                                        <  1:497483/251‑1 (MQ=255)
                       agagAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgcgc                                                                                                                                                                                                                               <  2:486576/251‑1 (MQ=255)
                                  cttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGGGCGCAGAGAgggggt                                                                                                                                                                                                                    >  1:508982/1‑248 (MQ=255)
                                                     gaTGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCCTGAGAAGTGAGGGGAGAGAGGG‑GTTTT‑CTTGCctgactgggatg                                                                                                                                                                                               >  2:285930/1‑246 (MQ=255)
                                                     gaTGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTAT‑CTTGCctgactg                                                                                                                                                                                                    <  1:285930/246‑1 (MQ=255)
                                                             aGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGGTTATTCAGAAAGTTTGTGGGTAACCCTTATAAGCTCCAGACCTATACCAGTCAGGAGTGATTTCAGGATGCCATGTCGCGGATCTCTTTTGCGGCGGTCATTTATTCGTTTTCTGCCATGAGAAGGGAGCGCAGAGAGGT‑ATTAT‑CTTGCCTGACTGAACTGGCGATTaa                                                                                                                                                                                       <  1:479288/251‑1 (MQ=255)
                                                                       acaGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTAT‑CTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGt                                                                                                                                                                            <  2:508982/251‑1 (MQ=255)
                                                                         aGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTAT‑CTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTAc                                                                                                                                                                          <  2:240237/251‑1 (MQ=255)
                                                                           ccGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTAT‑CTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTAccc                                                                                                                                                                        <  2:464167/251‑1 (MQ=255)
                                                                                     aagaGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCGCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTAT‑CTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGt                                                                                                                                                              <  2:271869/251‑1 (MQ=255)
                                                                                      agagAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGCAGTGAGCGCAGAGAGGG‑AGTAT‑CTTTCCTGACTTAACTGGCGGCTAA‑GTTTCCGCGGACCCGGCGTTTAGtt                                                                                                                                                             >  2:423894/1‑251 (MQ=255)
                                                                                      agagAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTAT‑CTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGtt                                                                                                                                                             <  1:412536/251‑1 (MQ=255)
                                                                                                cGCACTCATAGAAAAATTCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTAT‑CTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGgat                                                                                                                                                    <  2:182439/250‑1 (MQ=255)
                                                                                                cGCAATCATAGAAAAATGCGTCATTAGTAATATCTGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGATCCATACCTATAACAGTCAGGATTCATTTCATGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTAT‑CTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCTTACCCGGCGTTTAGTTATTGCGgat                                                                                                                                                    <  2:365250/250‑1 (MQ=255)
                                                                                                         aGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTTCCCTGAGGAGTGGACGCAGAGAGGG‑ATTAT‑CTTTCCTGGCTGAAATGGGGATTAAAGTTTC‑GCGTAACCCGCGTTTTGTTATTTCGgatggataaattt                                                                                                                                          >  1:307719/1‑242 (MQ=255)
                                                                                                          gAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAGCGCAGAGAGGG‑GTTAT‑CTTGCCTTACTGAAATGGCGATTTA‑GGTTTCGCGTCCCCGGGGTTTTGTTATTGGGGATGATG‑ATATTg                                                                                                                                         >  2:421647/1‑250 (MQ=255)
                                                                                                           aaaaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGGGAGGG‑GTTAT‑CTTGCCTGACTGAACTGGCGGTTAA‑GGTTTCGCGTACCCGGCGTTTAGGTATTGCGGATGATG‑GTAttta                                                                                                                                        >  1:251317/1‑248 (MQ=255)
                                                                                                                  cGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCCGAGAGGG‑GTTAT‑ATTTCCTGACTGAACTGGCGGTTAA‑GGTTTCGCGTCCCCGGCGTTTTGTTTTTGCGGAGGGTG‑ATATTGAAGCTccgc                                                                                                                                >  2:55790/1‑248 (MQ=255)
                                                                                                                  cGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGGGAGCGCAGAGAGGG‑GTTTT‑CTTTCCCTCCTGGACTGGCGGTTAA‑GTTTCCGCGTACCCCGCGGTTTGTTTTTGCGGATGTTG‑ATTTTGAAGCTCGGc                                                                                                                                >  1:304826/1‑251 (MQ=255)
                                                                                                                  cGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTAT‑CTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATG‑ATATTGAAGCTCGGc                                                                                                                                <  2:390644/251‑1 (MQ=255)
                                                                                                                   gTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTGATTTTTTTTTTTTTTTCCCTGAGAAGTGAGGGCAGAGGGGG‑GTTAT‑TTTGCCTGGCTGAACTGGGGATTAA‑GGTTCCGCGTAACCGGCGGTTTGGTATTTCGGGGGATG‑ATATTGTAGCTCGGcg                                                                                                                               >  2:548884/1‑250 (MQ=255)
                                                                                                                          gtagtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTAT‑CTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATG‑ATATTGAAGCTCGGCTGATTGGt                                                                                                                        <  1:133604/251‑1 (MQ=255)
                                                                                                                                  ggTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGACATTTTTTCTTTTTTTGCCATTAGAAGGGGGCGCAGAGAGGG‑ATTAT‑CTTTCCTGACTGAACTGGCGATTAA‑GATTTCGCGTACCCCGCGTTTAGTTATTGCGGGTGATG‑ATATTGAAGCTCGGCTGGTTGGTTccttttt                                                                                                                >  1:483346/1‑245 (MQ=255)
                                                                                                                                     aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACACTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCATGAGGAGCGAGGGGAGAGGGGGGATGTT‑CT‑GCCTGACTGAAATGGCGGTTTA‑GCTTCCGCGTACCCGGAGGTTTGCTATTGCGGGTGATG‑ATATTGAAGCTCGGCTGATTGGGTCATTTTCGcc                                                                                                             >  1:504074/1‑251 (MQ=255)
                                                                                                                                                    aGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGGGAAGTGGGCGCAGAGAGGG‑ATTAT‑CTTGCCTGCCTGGACTGGCGATTAA‑GTTTCCGCGTATCCGGCGTTTTGTTATTGCGGATGATG‑ATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACgg                                                                                              >  1:346545/1‑251 (MQ=255)
                                                                                                                                                           ttGCGGTTAACCCATATAAGCTCCATACCTATTCCAGGCATCAGTCATGTCAGCATGCTATGTGTCGGATCGCGGTTGCGTCTGTCATTTTTTCGTTGTCTGCCATGAGAAGTGAGCTCAGAGAGGG‑ATTAT‑CTTGCCTGACTGAACTGGCGATTAA‑GTGTCCGCCTACCCGGCGTTTAGTGATTGCGGATGATG‑ATATTGATGCTCGGCTGAGTGGTTCATTGTCGCCTTCTCCGCTGGACGGTGTAtt                                                                                        <  2:524265/250‑1 (MQ=255)
                                                                                                                                                                     cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTAT‑CTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATG‑ATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGt                                                                             <  1:421647/251‑1 (MQ=255)
                                                                                                                                                                       cTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCATGAGAAGGGAGCGCAGAGAGGG‑ATTAT‑CTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATG‑ATATTTAAGCTCGGCTGATTGGTTTATTGTCGCCATCACCGCTGGGCGGGGTATTTAGTAAAGCGTc                                                                            >  1:66075/1‑250 (MQ=255)
                                                                                                                                                                           taAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTCTGCCATGAGAAGTGCGCGCCGCGAGGG‑GGTTT‑TTTGCCTGACCGAACTGGGGATTTA‑GTTTCCCCCTACCCCGCGGTTAGTTATTTCCGATGATGGAATTTGAA‑CGCGGCCGGATGGTTCATTTGCGCCATCACCGCTGGTCGgtgtgtttaggaaaacgaccacgc                                                                       >  2:399735/1‑231 (MQ=255)
                                                                                                                                                                            aaGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTCTTTTTTTGTCCTGAGAAGTGAGCGCAGAGGGGG‑GTTAT‑CTTGCCTGACTGAACTGGCGATTTA‑GTTTTCGCGGCCCCCGCGGTTTGTTTTTTCGGGTGGTG‑GTATTTAAGCCTGGCTGATTGGGTTATTTTTGCCATCACCCCTGGGCgggggttttaggaaaacggccaagct                                                                      >  2:473657/1‑227 (MQ=255)
                                                                                                                                                                                tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGCCATTTTTTTTTTTTTTTCCATGAGAAAGGGGGGCAGAGAGGG‑GTTAT‑TTTGCCTGACTGAACTGGCGATTTA‑GTTTTCGCGGACACGGCGTTTTGTTATTGTGGGTGGTT‑ATATTTAAGCTCGGCTGATTTGTTTTTTTTTGGCATCACCGCTGGCAGGGGGATTTAGTAAAGCGACAACGCTgg                                                                    >  2:286261/1‑249 (MQ=255)
                                                                                                                                                                                         ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTCTTTTTTTTCCACGAGAAATGAGGGACGAGAGGG‑ATTTT‑CTTTCCTGTCTGAACTTGCGCTTAA‑GTTTTTGCCTACCCGGCGTTTTGTTATTTCCGGTGATG‑ATATTTTAGCTCGGCTGATTTGTTCCTTTTCGCCATCACACCTTGTCGGGGTCTTTAGTAAAGCGTCCACGCTggcggattttcat                                                         >  1:11791/1‑240 (MQ=255)
                                                                                                                                                                                         ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTAT‑CTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATG‑ATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTtca                                                          <  1:423894/250‑1 (MQ=255)
                                                                                                                                                                                           aCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAGTTTTTTTTTTTTTGCCATGAGAAGGGAGGGCAGAGGGGG‑GTTAT‑CTTGCCTGACTGAACTGGCGGTTTA‑GTTTTCGCGTACCCGGGGTTTAGTTATTGCGGGTGATG‑ATATTTAAGCTCGGGTGGTTTGGTCATTGTAGCCCTCACCGCTGGGCGGGGTATTAAGTAAAGGGTCAACGCTGGGGATTTTtcatc                                                        >  1:629547/1‑250 (MQ=255)
                                                                                                                                                                                           aCCAGTCAGGAGTCATTTCAGGAGGACATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTAT‑CTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATG‑ATATTGAAGCTCGGCGGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTtcatc                                                        <  2:77892/250‑1 (MQ=255)
                                                                                                                                                                                             cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTCTTTTTTTGCCCTGCGGAGGGAGCGCGGAGAGGG‑ATTATTCTT‑CCTGACTTTACTGGCGACTTA‑GTTTCCGCGTACCCGGCGTTTTGTTTTTTCGGGAGATG‑AAATTGAAGGTCGGGTGGTTTGTTTATTTTCGCCATCACCGCCGGGCGGGGTATTTAGTAAAGAGTCAACGCTGGAGACTTTTCCTCAgg                                                     >  1:78567/1‑251 (MQ=255)
                                                                                                                                                                                                 cAGGATTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTAT‑CTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATG‑ATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAAct                                                 <  1:151092/251‑1 (MQ=255)
                                                                                                                                                                                                 cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTAT‑CTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATG‑ATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAAct                                                 <  2:392024/251‑1 (MQ=255)
                                                                                                                                                                                                      gTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTAT‑CTTGCCCGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATG‑ATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGc                                                                       <  1:524576/224‑1 (MQ=255)
                                                                                                                                                                                                             cAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTAT‑CTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATG‑ATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCAtt                                     <  1:609071/251‑1 (MQ=255)
                                                                                                                                                                                                                          cgcgGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTAT‑CTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATG‑ATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGt                        <  1:161358/251‑1 (MQ=255)
                                                                                                                                                                                                                             ggATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGGAGAGAGCGCAGAGAGGG‑GTTAT‑CTTGCCTGACTGAACTGGCGATTAA‑GGTTCCGCGTACCCGGCGTTTTGTTATTGCGGGTGATG‑ATATTGAAGCTCGGCTGATTTGTTCATTGTCGCCCTCACCGCTGGACGGTGGATTAAGTAAAGCGTCAACGCTGGAGGTTTTTTATCAGGAACTCTTTTTGTCATTAAATGGTGGACGTc                       >  1:366362/1‑249 (MQ=255)
                                                                                                                                                                                                                                          ggcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTAT‑CTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATG‑ATATTGAAGCTCGGCTGATTGGTTCAATGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGcgaccgaccgac        <  1:514389/251‑1 (MQ=255)
                                                                                                                                                                                                                                            cggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTAT‑CTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATTATG‑ATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGTCACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTg      >  1:301363/1‑251 (MQ=255)
                                                                                                                                                                                                                                            cggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTAT‑CTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATG‑ATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTg      >  1:242365/1‑251 (MQ=255)
                                                                                                                                                                                                                                            cggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTAT‑CTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATG‑ATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTg      >  2:64169/1‑251 (MQ=255)
                                                                                                                                                                                                                                                cATTTGTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTAT‑CTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATG‑ATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCATGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAAc   <  2:201753/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 aTTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTAT‑CTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATG‑ATATTGAAGCTCGGCTGCTTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACa  >  1:580678/1‑250 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
GGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTAT‑CTTGCCTGACTGAACTGGCGATTAA‑GTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATG‑ATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACA  >  NC_000913/4604096‑4604586

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