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

AAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATT  >  NC_000913/4604100‑4604600
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                              
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)
                                                                                                                                                                                                                                                      ttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTgg        <  2:346545/249‑1 (MQ=255)
                                                                                                                                                                                                                                                          ctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACAtt  >  2:253797/1‑251 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                              
AAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATT  >  NC_000913/4604100‑4604600

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: