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.

GGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACA  >  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)
                                                     gaTGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCCTGAGAAGTGAGGGGAGAGAGGGGTTTTCTTGCctgactgggatg                                                                                                                                                                                               >  2:285930/1‑246 (MQ=255)
                                                     gaTGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctgactg                                                                                                                                                                                                    <  1:285930/246‑1 (MQ=255)
                                                             aGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGGTTATTCAGAAAGTTTGTGGGTAACCCTTATAAGCTCCAGACCTATACCAGTCAGGAGTGATTTCAGGATGCCATGTCGCGGATCTCTTTTGCGGCGGTCATTTATTCGTTTTCTGCCATGAGAAGGGAGCGCAGAGAGGTATTATCTTGCCTGACTGAACTGGCGATTaa                                                                                                                                                                                       <  1:479288/251‑1 (MQ=255)
                                                                       acaGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGt                                                                                                                                                                            <  2:508982/251‑1 (MQ=255)
                                                                         aGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAc                                                                                                                                                                          <  2:240237/251‑1 (MQ=255)
                                                                           ccGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAccc                                                                                                                                                                        <  2:464167/251‑1 (MQ=255)
                                                                                     aagaGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCGCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGt                                                                                                                                                              <  2:271869/251‑1 (MQ=255)
                                                                                      agagAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGCAGTGAGCGCAGAGAGGGAGTATCTTTCCTGACTTAACTGGCGGCTAAGTTTCCGCGGACCCGGCGTTTAGtt                                                                                                                                                             >  2:423894/1‑251 (MQ=255)
                                                                                      agagAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGtt                                                                                                                                                             <  1:412536/251‑1 (MQ=255)
                                                                                                cGCACTCATAGAAAAATTCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgat                                                                                                                                                    <  2:182439/250‑1 (MQ=255)
                                                                                                cGCAATCATAGAAAAATGCGTCATTAGTAATATCTGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGATCCATACCTATAACAGTCAGGATTCATTTCATGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCTTACCCGGCGTTTAGTTATTGCGgat                                                                                                                                                    <  2:365250/250‑1 (MQ=255)
                                                                                                         aGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTTCCCTGAGGAGTGGACGCAGAGAGGGATTATCTTTCCTGGCTGAAATGGGGATTAAAGTTTCGCGTAACCCGCGTTTTGTTATTTCGgatggataaattt                                                                                                                                          >  1:307719/1‑242 (MQ=255)
                                                                                                          gAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTTACTGAAATGGCGATTTAGGTTTCGCGTCCCCGGGGTTTTGTTATTGGGGATGATGATATTg                                                                                                                                         >  2:421647/1‑250 (MQ=255)
                                                                                                           aaaaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGGGAGGGGTTATCTTGCCTGACTGAACTGGCGGTTAAGGTTTCGCGTACCCGGCGTTTAGGTATTGCGGATGATGGTAttta                                                                                                                                        >  1:251317/1‑248 (MQ=255)
                                                                                                                  cGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCCGAGAGGGGTTATATTTCCTGACTGAACTGGCGGTTAAGGTTTCGCGTCCCCGGCGTTTTGTTTTTGCGGAGGGTGATATTGAAGCTccgc                                                                                                                                >  2:55790/1‑248 (MQ=255)
                                                                                                                  cGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGGGAGCGCAGAGAGGGGTTTTCTTTCCCTCCTGGACTGGCGGTTAAGTTTCCGCGTACCCCGCGGTTTGTTTTTGCGGATGTTGATTTTGAAGCTCGGc                                                                                                                                >  1:304826/1‑251 (MQ=255)
                                                                                                                  cGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGc                                                                                                                                <  2:390644/251‑1 (MQ=255)
                                                                                                                   gTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTGATTTTTTTTTTTTTTTCCCTGAGAAGTGAGGGCAGAGGGGGGTTATTTTGCCTGGCTGAACTGGGGATTAAGGTTCCGCGTAACCGGCGGTTTGGTATTTCGGGGGATGATATTGTAGCTCGGcg                                                                                                                               >  2:548884/1‑250 (MQ=255)
                                                                                                                          gtagtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGt                                                                                                                        <  1:133604/251‑1 (MQ=255)
                                                                                                                                  ggTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGACATTTTTTCTTTTTTTGCCATTAGAAGGGGGCGCAGAGAGGGATTATCTTTCCTGACTGAACTGGCGATTAAGATTTCGCGTACCCCGCGTTTAGTTATTGCGGGTGATGATATTGAAGCTCGGCTGGTTGGTTccttttt                                                                                                                >  1:483346/1‑245 (MQ=255)
                                                                                                                                     aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACACTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCATGAGGAGCGAGGGGAGAGGGGGGATGTTCTGCCTGACTGAAATGGCGGTTTAGCTTCCGCGTACCCGGAGGTTTGCTATTGCGGGTGATGATATTGAAGCTCGGCTGATTGGGTCATTTTCGcc                                                                                                             >  1:504074/1‑251 (MQ=255)
                                                                                                                                                    aGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGGGAAGTGGGCGCAGAGAGGGATTATCTTGCCTGCCTGGACTGGCGATTAAGTTTCCGCGTATCCGGCGTTTTGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACgg                                                                                              >  1:346545/1‑251 (MQ=255)
                                                                                                                                                           ttGCGGTTAACCCATATAAGCTCCATACCTATTCCAGGCATCAGTCATGTCAGCATGCTATGTGTCGGATCGCGGTTGCGTCTGTCATTTTTTCGTTGTCTGCCATGAGAAGTGAGCTCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTGTCCGCCTACCCGGCGTTTAGTGATTGCGGATGATGATATTGATGCTCGGCTGAGTGGTTCATTGTCGCCTTCTCCGCTGGACGGTGTAtt                                                                                        <  2:524265/250‑1 (MQ=255)
                                                                                                                                                                     cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGt                                                                             <  1:421647/251‑1 (MQ=255)
                                                                                                                                                                       cTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCATGAGAAGGGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTTAAGCTCGGCTGATTGGTTTATTGTCGCCATCACCGCTGGGCGGGGTATTTAGTAAAGCGTc                                                                            >  1:66075/1‑250 (MQ=255)
                                                                                                                                                                           taAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTCTGCCATGAGAAGTGCGCGCCGCGAGGGGGTTTTTTGCCTGACCGAACTGGGGATTTAGTTTCCCCCTACCCCGCGGTTAGTTATTTCCGATGATGGAATTTGAACGCGGCCGGATGGTTCATTTGCGCCATCACCGCTGGTCGgtgtgtttaggaaaacgaccacgc                                                                       >  2:399735/1‑231 (MQ=255)
                                                                                                                                                                            aaGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTCTTTTTTTGTCCTGAGAAGTGAGCGCAGAGGGGGGTTATCTTGCCTGACTGAACTGGCGATTTAGTTTTCGCGGCCCCCGCGGTTTGTTTTTTCGGGTGGTGGTATTTAAGCCTGGCTGATTGGGTTATTTTTGCCATCACCCCTGGGCgggggttttaggaaaacggccaagct                                                                      >  2:473657/1‑227 (MQ=255)
                                                                                                                                                                                tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGCCATTTTTTTTTTTTTTTCCATGAGAAAGGGGGGCAGAGAGGGGTTATTTTGCCTGACTGAACTGGCGATTTAGTTTTCGCGGACACGGCGTTTTGTTATTGTGGGTGGTTATATTTAAGCTCGGCTGATTTGTTTTTTTTTGGCATCACCGCTGGCAGGGGGATTTAGTAAAGCGACAACGCTgg                                                                    >  2:286261/1‑249 (MQ=255)
                                                                                                                                                                                         ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTCTTTTTTTTCCACGAGAAATGAGGGACGAGAGGGATTTTCTTTCCTGTCTGAACTTGCGCTTAAGTTTTTGCCTACCCGGCGTTTTGTTATTTCCGGTGATGATATTTTAGCTCGGCTGATTTGTTCCTTTTCGCCATCACACCTTGTCGGGGTCTTTAGTAAAGCGTCCACGCTggcggattttcat                                                         >  1:11791/1‑240 (MQ=255)
                                                                                                                                                                                         ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTtca                                                          <  1:423894/250‑1 (MQ=255)
                                                                                                                                                                                           aCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAGTTTTTTTTTTTTTGCCATGAGAAGGGAGGGCAGAGGGGGGTTATCTTGCCTGACTGAACTGGCGGTTTAGTTTTCGCGTACCCGGGGTTTAGTTATTGCGGGTGATGATATTTAAGCTCGGGTGGTTTGGTCATTGTAGCCCTCACCGCTGGGCGGGGTATTAAGTAAAGGGTCAACGCTGGGGATTTTtcatc                                                        >  1:629547/1‑250 (MQ=255)
                                                                                                                                                                                           aCCAGTCAGGAGTCATTTCAGGAGGACATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCGGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTtcatc                                                        <  2:77892/250‑1 (MQ=255)
                                                                                                                                                                                             cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTCTTTTTTTGCCCTGCGGAGGGAGCGCGGAGAGGGATTATTCTTCCTGACTTTACTGGCGACTTAGTTTCCGCGTACCCGGCGTTTTGTTTTTTCGGGAGATGAAATTGAAGGTCGGGTGGTTTGTTTATTTTCGCCATCACCGCCGGGCGGGGTATTTAGTAAAGAGTCAACGCTGGAGACTTTTCCTCAgg                                                     >  1:78567/1‑251 (MQ=255)
                                                                                                                                                                                                 cAGGATTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAAct                                                 <  1:151092/251‑1 (MQ=255)
                                                                                                                                                                                                 cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAAct                                                 <  2:392024/251‑1 (MQ=255)
                                                                                                                                                                                                      gTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCCGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGc                                                                       <  1:524576/224‑1 (MQ=255)
                                                                                                                                                                                                             cAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCAtt                                     <  1:609071/251‑1 (MQ=255)
                                                                                                                                                                                                                          cgcgGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGt                        <  1:161358/251‑1 (MQ=255)
                                                                                                                                                                                                                             ggATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGGAGAGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGATTAAGGTTCCGCGTACCCGGCGTTTTGTTATTGCGGGTGATGATATTGAAGCTCGGCTGATTTGTTCATTGTCGCCCTCACCGCTGGACGGTGGATTAAGTAAAGCGTCAACGCTGGAGGTTTTTTATCAGGAACTCTTTTTGTCATTAAATGGTGGACGTc                       >  1:366362/1‑249 (MQ=255)
                                                                                                                                                                                                                                          ggcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCAATGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGcgaccgaccgac        <  1:514389/251‑1 (MQ=255)
                                                                                                                                                                                                                                            cggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATTATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGTCACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTg      >  1:301363/1‑251 (MQ=255)
                                                                                                                                                                                                                                            cggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTg      >  1:242365/1‑251 (MQ=255)
                                                                                                                                                                                                                                            cggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTg      >  2:64169/1‑251 (MQ=255)
                                                                                                                                                                                                                                                cATTTGTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCATGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAAc   <  2:201753/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 aTTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGCTTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACa  >  1:580678/1‑250 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
GGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACA  >  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: