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

AATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTATCTTGCCTGACTGAACTGGCG‑ATTAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAG‑ATTTTTCATC‑AGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTG  >  NC_000913/4604101‑4604593
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                              
aaTAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttcttttttt                                                                                                                                                                                                                                                         >  1:23798/1‑249 (MQ=255)
aaTAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttct                                                                                                                                                                                                                                                         <  1:336496/251‑1 (MQ=255)
 aTAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑ATctttttttttttttt                                                                                                                                                                                                                                                       >  2:23518/1‑247 (MQ=255)
  tAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttcttc                                                                                                                                                                                                                                                       >  1:311350/1‑249 (MQ=255)
   aGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGcc                                                                                                                                                                                                                                                      >  1:112233/1‑251 (MQ=255)
        gcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGTTCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATgag                                                                                                                                                                                                                                                 <  1:18573/251‑1 (MQ=255)
              tGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTCTTTTTTTTCCCTGAGAAGTgg                                                                                                                                                                                                                                           >  1:10115/1‑250 (MQ=255)
               gCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAtttttttttttttttccattggaagag                                                                                                                                                                                                                                            >  2:163486/1‑234 (MQ=255)
                    agaACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACAAAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCag                                                                                                                                                                                                                                     <  1:373300/251‑1 (MQ=255)
                      aacttaacttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttttc                                                                                                                                                                                                                                                          <  1:314620/228‑1 (MQ=255)
                      aacttaacttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCAtttttttttttttcaaacgggagaacggcgtgttgg                                                                                                                                                                                                                                   >  2:314620/1‑227 (MQ=255)
                                       gTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTATCTTGCctg                                                                                                                                                                                                                 <  1:334266/251‑1 (MQ=255)
                                              gagaTGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAGGGGAGAGGGGG‑GTTATCTTTCCTGGCTGaa                                                                                                                                                                                                           >  2:362192/1‑250 (MQ=255)
                                                     cGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGGGCGCCGGGAGGG‑GTTTTCttttctgtcgtaaattgggg                                                                                                                                                                                                    >  2:339067/1‑232 (MQ=255)
                                                                               gaagaGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGGGAAGGGAGAGCAGAGAGGG‑GTTATTTTTTCTGCGTGAACTGGCG‑ATTAAG‑ttttcccggcacccgcggtttg                                                                                                                                                                       >  1:362504/1‑232 (MQ=255)
                                                                                           cGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCCTGAGAAGTGAGCGCAGAGAGGG‑AGTATCTTGCCTGACTGAAATGGCG‑ATTAAG‑TTTCCGCGGAACCGGCGGTTTGTTATTTCGgatg                                                                                                                                                           >  1:405637/1‑251 (MQ=255)
                                                                                                        aaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCCGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTCTTTTTTTTCCATGAGAAGGGAGCGCCGCGGGGG‑GTTATCTTGCCTGTCTGCACTGGCC‑ATTTCG‑ATTCCTCCGAACCGGCGGTTTGTTAATGGGGATGATGATATTggaa                                                                                                                                               >  1:505525/1‑247 (MQ=255)
                                                                                                        aaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTATCTTGCCTGTCTGAACTGGCG‑ATTAAG‑TTTTCGCGGACCCGGCGTTTTGTTATTGCGGGGGATGTTATTTAAGc                                                                                                                                              >  2:150493/1‑251 (MQ=255)
                                                                                                        aaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTCTTTCATGGGAAGTGAGCGCCGAGAGGGGATT‑TTTTTCCTTGCTGAACCGGGG‑GTTAAG‑TTTTCGCGTACCCGGCGTTTTGGTTTTTCGGGTGATGATAtttaaac                                                                                                                                              >  2:126403/1‑246 (MQ=255)
                                                                                                            gCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTATCTTGCCTGACTGAACTGGCG‑ATTAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCg                                                                                                                                           <  2:10115/250‑1 (MQ=255)
                                                                                                              gTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTATCTTGCCTGACTGAACTGGCG‑ATTAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGTAGCTCGGCt                                                                                                                                        <  2:318107/251‑1 (MQ=255)
                                                                                                                                nttGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTAATTTTTTTTTTTTTTGCCATGAGAAGTGGGCGCCGAGAGGG‑ATTATCTTTCCTGACTCAACTGGCG‑CTTAAG‑GTTTCTCGGACCCGGCGGTTTGTTATTGTGGATGATGATATTGTAGCTCGGCTgctttggttattgttgcc                                                                                                                      >  2:242442/2‑234 (MQ=255)
                                                                                                                                          tATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTATCTTGCCTGACTGAACTGGCG‑ATTAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCt                                                                                                             <  2:112233/250‑1 (MQ=255)
                                                                                                                                               aGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTATCTTGCCTGACTGAACTGGCG‑ATTAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACgg                                                                                                       <  2:222607/251‑1 (MQ=255)
                                                                                                                                                    gTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑GTTATCTTGCCTGACTGGAATGGCG‑GTTAAG‑TTTTCGCGTACCCCGCGTTTAGTTATTTCGGGTGATGATATTGAAGCTCGGCGGATTGGTTCATTGTCGCCATCCCCCCCGGACgggggat                                                                                                  >  1:27202/1‑246 (MQ=255)
                                                                                                                                                              aaCCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTCTGCCCTGAGAAGTGAGCGCAGAGGGGG‑ATTTTTTTGCCTGTCCGAACTGGCGGTTTAAG‑TT‑CCGCGGACCCGGCGGTTTGGTATTTCGGGTGGTGGTATTTAAGCTCGGCTGGGTGGGTTATTTTCGCCCTTACCGCTGGGCGGTGTATTTAATAAAgg                                                                                        >  2:144920/1‑250 (MQ=255)
                                                                                                                                                                          cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGGGAGCGCCGAGAGGG‑GTTATCTTTCCTGACTTGACTGGCG‑ATTAAG‑TTTCCGCGTACCCCGCGTTTTGTTATTGCGGGTGATGGTATTGTAGCTCGGCTGGTTGGTTTCTTTTCGCCATCACCGCCGGGCGGGGTATTTAATAAAGCGGCAAAGCTggg                                                                            >  2:272753/1‑250 (MQ=255)
                                                                                                                                                                              aTACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTATCTTGCCTGACTGAACTGGCG‑ATTAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAG‑Att                                                                       <  1:150493/251‑1 (MQ=255)
                                                                                                                                                                                aCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTATCTTGCCTGACTGAACTTGCG‑ATTAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAG‑Atttt                                                                     <  2:382374/251‑1 (MQ=255)
                                                                                                                                                                                      aCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGGGAGCGCAGAGAGGG‑ATTATCTTGCCTGACTGAACTGGCG‑GTTTAGATTTCC‑CGTACCCCGCGGTTAGTTATTGCGGGTGTTGATATTTAAAATCGGCTGATTGGTTCATTGTTGCCATCACCGCTGGGCGGTGTATTTAGTAAAGCGGCAACGCTGGGG‑AGTTTtgatc                                                                >  2:269641/1‑250 (MQ=255)
                                                                                                                                                                                        cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGCCATTTTTTTTTTTTTTGCCATGAGAAGAGAGGCCAGGGGGGG‑GTTATTTTTCCTGGCTGAAATGGCG‑AGTAAG‑TTTCCGCGCACCCGGCGTTTAAGTATTGTGGATGATGATATTGCAGCTTGGGTGCATGGGTCATTTTCGCCATCACCGCTGGGAGGGGTATTAAGTAAAACGGCAAAGATGGAGGATTTTtcttagg                                                              >  1:401797/1‑246 (MQ=255)
                                                                                                                                                                                            cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTATCTTGCCTGACTGAACTGGCG‑ATTAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAG‑ATTTTTCATC‑AGGAAct                                                        <  2:401797/251‑1 (MQ=255)
                                                                                                                                                                                               gAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCCTGAGAAGTGGGCGCAGAGAGGG‑GTTATCTTTCCTGACTGAACTGGCG‑ATTAAG‑TTTTCGCGTACCCGGCGGTTTGTTTTTTCGGGTGATGATATTGAAGCTCGGCTGATTGGTTTATTGTCGCCATCACCGCTTGACGGTGTATTAAGTAAAGCGTCAACGCTGGAG‑ATTTTTCATC‑AGGAAATCtt                                                     >  1:350003/1‑251 (MQ=255)
                                                                                                                                                                                                aGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTATCTTGCCTGACTGAACTGGCG‑ATTAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAG‑ATTTTTCATC‑AGGAACTCttt                                                    <  2:405637/251‑1 (MQ=255)
                                                                                                                                                                                                   cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTATCTTGCCTGACTGAACTGGCG‑ATTAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAG‑ATTTTTCATC‑AGGAACTCttttt                                                  <  1:288583/250‑1 (MQ=255)
                                                                                                                                                                                                          ggATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTATCTTGCCTGACTGAACTGGCG‑ATTAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGGAAAGCGTCAACGCTGGAG‑ATTTTTCATC‑AGGAACTCTTTTTGTCATTaa                                          <  1:228626/251‑1 (MQ=255)
                                                                                                                                                                                                           gATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTATCTTGCCTGACTGAACTGGCG‑ATTAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAG‑ATTTTTCATC‑AGGAACTCTTTTTGTCATTaaa                                         <  2:311350/251‑1 (MQ=255)
                                                                                                                                                                                                           gATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTATCTTGCCTGACTGAACTGGCG‑ATTAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAG‑ATTTTTCATC‑AGGAACTCTTTTTGTCATTaaa                                         <  1:33235/251‑1 (MQ=255)
                                                                                                                                                                                                           gATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTATCTTGCCTGACTGAACTGGCG‑ATTAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAG‑ATTTTTCATC‑AGGAACTCTTTTTGTCATTaaa                                         <  1:362192/251‑1 (MQ=255)
                                                                                                                                                                                                               ccATGTCGCGGATCTCGTTTGCGGCAGTCATTTTTTCTTTTTCGGCCATGAGAAGTGAGCGCAGAGATGG‑ATTATCTTGCCTGACTGAACTGGCG‑ATTAAG‑TTTCCGCGTACCCGGCGTTGAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAG‑ATTTTTCATC‑AGGAACTCTTTTTTTCATTAAATgg                                      <  2:304049/250‑1 (MQ=255)
                                                                                                                                                                                                                       cgGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTATCTTGCCTGACTGAACTGGCG‑ATTAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAG‑ATTTTTCATC‑AGGAACTCTTTTTGTCATTAAATGGTGTACGTCa                             <  1:272753/251‑1 (MQ=255)
                                                                                                                                                                                                                        ggATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTATCTTGCCTGACTGAACTGGCG‑ATTAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAG‑ATTTTTCATC‑AGGAACTCTTTTTGTCATTAAATGGTGTACGTCAg                            <  1:339067/251‑1 (MQ=255)
                                                                                                                                                                                                                          aTCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTCCCCTCAGAAGGGAGCGCAGAGAGGG‑ATTTTCTTGCCTGACTGCACCGGCG‑ATTTAG‑TTTTCGCGTACCCGGCGTTTTGTTATTTCGGGTGATGATATTTTAGCTCGGCTGATTGGTTCATTGTTGCCCTCACCCCTGGACGGGGTATTTAGTAAAACGTCCACGCTGGTG‑ATTTTTTATC‑AGGAACTCTTTTTGTCATTTAATGGTGTACGTccggg                          >  1:124898/1‑247 (MQ=255)
                                                                                                                                                                                                                                   gcggcgGTCATTTTTTCTTTTTCTGCCATGAGAAGGGAGCGCAGAGGGGG‑ATTATCTTTGCTGACTGAACTGGCG‑ATTAAG‑TTTTCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTTTTGCCATCCCCGCTGGACGGGGGATTTAGAAAAAGGGCAACGCTTGGG‑GTTTTTCATCCGGGAAC‑CTTTTTTTCCTTAAAAGGGGTACGTCAGGcgaccgaccg                 >  2:561108/1‑251 (MQ=255)
                                                                                                                                                                                                                                      gcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTATCTTGCCTGACTGAACTGGCG‑ATTAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAG‑ATTTTTCATC‑AGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACt              >  1:231366/1‑251 (MQ=255)
                                                                                                                                                                                                                                         gTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTATCTTGCCTGACTGAACTGGCG‑ATTAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAG‑ATTTTTCATC‑AGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGa            >  2:421828/1‑250 (MQ=255)
                                                                                                                                                                                                                                           cATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTATCTTGCCTGACTGAACTGGCG‑ATTAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGATGGTGTATTAAGTAAAGCGTCAACGCTGGAG‑ATTTTTCATC‑AGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAAc          >  2:60079/1‑250 (MQ=255)
                                                                                                                                                                                                                                           cATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTATCTTGCCTGACTGAACTGGCG‑ATTAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAG‑ATTTTTCATC‑AGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAAc          >  1:450262/1‑250 (MQ=255)
                                                                                                                                                                                                                                            aTTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTATCTTGCCTGACTGAACTGGCG‑ATTAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAG‑ATTTTTCAGC‑AGGAACTCTTTTTGACATTAAATGGTTTACGTAAGGCGACCGACCGACTGAACa         >  1:160908/1‑250 (MQ=255)
                                                                                                                                                                                                                                                  tctttttctGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTATCTTGCCTGACTGAACTGGCG‑ATTAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAG‑ATTTTTCATC‑AGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTg  >  2:27085/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                              
AATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGG‑ATTATCTTGCCTGACTGAACTGGCG‑ATTAAG‑TTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAG‑ATTTTTCATC‑AGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTG  >  NC_000913/4604101‑4604593

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