Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,502,756 C→T 19.2% E226K (GAA→AAA)  insG ← IS4 transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,502,7560CT19.2% 101.7 / 4.2 50E226K (GAA→AAA) insGIS4 transposase
Reads supporting (aligned to +/- strand):  ref base C (2/38);  new base T (8/2);  total (10/40)
Fisher's exact test for biased strand distribution p-value = 3.46e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAGATCCCCCATTTCTCCTCCGGGGAAGCGCATGGCGTCCGTCATCGACGTCAGCAGATGGCAGACTTTTCC‑TTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGACCATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGCATGTGTCTGGCGGGGGA  >  NC_000913/4502513‑4502986
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                       
ccacaTCCCCCTTTTCCCCTCCGGGGAAGCGCAGGGGGCCCGTCACCGACGCCAGCAGTTGGAAGCTTTTTCC‑TTTGCCGGTCCCGGTCACCAGGCGGGCAGTCCCTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACC‑‑‑‑‑‑‑CTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTtctc                                                                                                                                                                                                                                       <  2:610719/235‑1 (MQ=255)
       ccccATTTCTCCTCCGGGGAAGCGCATGGCGTCCGTCATCGACGTCAGCAGATGGCAGACTTTTCC‑TTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCTTGCCTGCGGGCTGGTTTTCAGCTTACCCCGATGATCGCCTTTACCCAGTTTTCTGGTCTCTTCATATTGCCCTCCCTTCCTCGGGGGATTCATCCAGTGGCGGGTTtct                                                                                                                                                                                                                                        >  2:81351/1‑238 (MQ=255)
       ccccATTTCCCCACCGGGGAAGCGCAGGGCGTCCGTCATGGACGTCAGCAGATGGCAGACTTTTCC‑TTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTtctc                                                                                                                                                                                                                                       <  1:292430/239‑1 (MQ=255)
         ccATTTCTCCTCCGGGGAAGCGCATGGCGTCCGTCATCGACGTCAGCAGATGGCAGACTTTTCC‑TTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTTTcc                                                                                                                                                                                                                                      >  1:117236/1‑238 (MQ=255)
         ccATTTCTCCTCCGGGGAAGCGCATGGCGTCCGTCATCGACGTCAGCAGATGGCAGACTTTTCC‑TTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACC‑‑‑‑‑‑‑CTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGGGGTTTTCCGCCCgg                                                                                                                                                                                                                               >  1:610719/1‑238 (MQ=255)
             ttCCCCCCCGGGGAAGCCCATGGCGGCCGTCAGCGGCGTCAGTAGATGGCAGTCTTTTCC‑TTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCa                                                                                                                                                                                                                                 <  1:724208/239‑1 (MQ=255)
              tctcCTCCGGGGAAGCGCATGGCGTCCGTCATCGACGTCAGCAGATGGCAGACTTTTCC‑TTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGCTGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTTCGCTCCCTTTCTGAGAGGTATCATCCAGTGGGGGCGGTTTCTCGcccc                                                                                                                                                                                                                                >  2:750506/1‑237 (MQ=255)
                   tccGGGGAAGCGAGTGGCGTACGTCGACGACGTCAGCAGATGGCAGACTTTTCC‑TTTGCGCGTAACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCCCTTTTTGCGGGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTcc                                                                                                                                                                                                                           <  1:628689/239‑1 (MQ=255)
                   cccGGGGAAGCGCATGGCCGCCGTCATCGCCGGCATCAGATGGCAGACTTTTCC‑TTTGCGCGTCCCGGTCAGCGGGCGGGCAGGCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTcc                                                                                                                                                                                                                           <  2:824813/238‑1 (MQ=255)
                             cccctGGCGTCCGTTCTCGGCGCCAGCAGATGGCCGCTTTTTCC‑TTTGCGCGTCCCGGTCCGCCGGCGGGCAGTCATTTCATTTCCCAGTCCCGGCCCCTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTaa                                                                                                                                                                                                                 <  2:395462/235‑1 (MQ=255)
                              gcgtGGCGCCGGTCTTCGTCGTCAGCAGATGGCAGACTGTGCC‑TTTGTGCGTCACGGTCAGCTTGCGGGCAGTCACTTCATTCCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGTGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTaa                                                                                                                                                                                                                 <  1:351081/235‑1 (MQ=255)
                              ccATGGGGTCCGGCATCGACGTTCGCAGTTTGCAAACTTTTCC‑GTGGCGCGTCACGGTCAGGAGGCGGGGCGTCCCTTCATTTCCCAGTCCCGGCCACTTTTTTCGTCCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTTTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTaa                                                                                                                                                                                                                 <  2:30030/237‑1 (MQ=255)
                                      gccgtcAACGGCGTCAGGAGATGGCAGATTTTTCC‑TTCGCGCCGCACGGGCAGCAGGCGGGCAGTCACTTCATTTCCTAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAg                                                                                                                                                                                                        <  1:61481/238‑1 (MQ=255)
                                       ccgtcATCGACGTCAGCAGATGGCAGACTTTTCC‑TTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCTCTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGt                                                                                                                                                                                                       <  1:340211/239‑1 (MQ=255)
                                        cgtcATCGACGTCAGCAGATGGCAGACTTTTCC‑TTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGGATAATCCAGTGGCGGTGTTTTCCCCCCCGGTCCCCGGCCATTTACAGTGCCCGTg                                                                                                                                                                                                      >  1:855322/1‑239 (MQ=255)
                                        cgccATCGACGTCAGCAGTGTGCAGCCTTTTCC‑TTTGGGCGTTCCGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCCCCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGt                                                                                                                                                                                                       <  1:582100/235‑1 (MQ=255)
                                              cggcgtcgtcagATGGCATAGTTGTCT‑TTTCGGTCTCACTGTCCTCAGGCGGGCAGTCACGTCATTTCCCAGTCCCGGCCCCTTTTTTCGTGGCTGTGGGCTGGTTTTCAGCTTCACCAGGTGATCGCTTTGGCCAAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCATGCATTTAACAGTCCCAGTGAGtaac                                                                                                                                                                                                <  2:151120/230‑2 (MQ=255)
                                              cggcGTCAGCAGATGGCAGACTTTTCC‑TTTGCGCGTCACGGTCAGCCGGCGGGCAGTCCCTTAATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGtaat                                                                                                                                                                                                <  1:889627/236‑1 (MQ=255)
                                                        agATGGCAGACTTTTCC‑TTTGCGCGTCCCGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTAt                                                                                                                                                                                      <  1:16262/238‑1 (MQ=255)
                                                             ctcgttttttCCTTTTGCG‑GTCCCGGTCCACGGGCGGGCATTCCCTTCTTTTCCCAGCCCCGGCCACTTTTTTCGGGCCTGCGGGCTGGTTTTCAGCTTCCCCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCATTa                                                                                                                                                                                <  2:261719/233‑1 (MQ=255)
                                                                   ttttCC‑TGTGCGCGTCGCGGTCAGCAGGGGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTGCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCATTAACGTCa                                                                                                                                                                          <  1:243073/239‑1 (MQ=255)
                                                                      tCC‑TTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTTTCCCCCCCGGCTCCAGGCACTTTACAGTCCCCGTGAGTAATAAA‑CTTTATCCATTTACGTCataa                                                                                                                                                                       >  1:676470/1‑236 (MQ=255)
                                                                                   cACGGTCAGCGGGCGGGCAGTCACTCCA‑TCCCCAGTCCCGGCCACTTTTTCCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCATTAACGTCAGAGTGTTATCGCCGGt                                                                                                                                                          <  1:506043/239‑1 (MQ=255)
                                                                                      ggTCCGCAGGCGGGCAGTCCCTTCATTTCCCAGTCCCGGCCTCTTTTTTCGTGCCTGCGGGCGGGTTTTCAGTTTCCCCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCATTAACGTCAGAGTGTTATCGCCGGttt                                                                                                                                                        <  2:219754/239‑1 (MQ=255)
                                                                                        tCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTTTCCCGCCCAGCCCCACGCCTTTTACAATACCAGTGGGTAATAACACTTTTTCCCTT‑ACGTCAAAGGGTTTTCGCTGGtcctt                                                                                                                                                      >  1:18080/1‑235 (MQ=255)
                                                                                          agcagGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTaa                                                                                                                                                                                                                 <  1:914524/179‑1 (MQ=255)
                                                                                           gaagggggtaaGTCTCTTTATTTCCCAGTCCCGGCCACTTTTTGCGTGCCTGCGGGCTGGTTTGCAGCTTCACCAGATGATCGCCTTTACCCAGTTTGCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCt                                                                                                                                                   <  2:445415/229‑1 (MQ=255)
                                                                                           ccgggggggCATTCTCTTCCTTTCCCAGCCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCCCCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCt                                                                                                                                                   <  2:676470/233‑1 (MQ=255)
                                                                                            cagGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTTTCCCGCCACGCTCCCGGCATTTTAAAGTACCCGTGAGTGAAAAA‑CTTTTTCCCATAACCTCAGAAGGTTTTCTCCggggttttttt                                                                                                                                                  >  1:722299/1‑230 (MQ=255)
                                                                                                         aCTTCATTTCCCAGGCCCGGCAACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCGGAGCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGc                                                                                                                                     <  1:81351/239‑1 (MQ=255)
                                                                                                           ttCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCAttt                                                                                                                                                                                                                   <  1:919557/160‑1 (MQ=255)
                                                                                                                      ggcccGGCCCTTTTTGGCGTGCCTTGGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCg                                                                                                                        <  1:348699/237‑1 (MQ=255)
                                                                                                                                               cGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTTTCCCGCCACGCTCCCGGCCTTTTAAAATCCCCGGGAGTAATAAC‑CCTTATCCCTTAACGTCAGAGTGTTTTCGCCGGGTTGTTCTATAAATTTCTCAGCAAACTCATTTTCTCTGTTCTTCATCGTGCCGAAGGCt                                                                                                >  1:904830/1‑238 (MQ=255)
                                                                                                                                                  gCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCa                                                                                             <  1:873907/238‑1 (MQ=255)
                                                                                                                                                              cTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCg                                                                                                                        <  1:300333/199‑1 (MQ=255)
                                                                                                                                                                   ccAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTg                                                                           <  1:169062/239‑1 (MQ=255)
                                                                                                                                                                       atgatCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCATTAACGTCaga                                                                                                                                                                        <  1:384994/142‑1 (MQ=255)
                                                                                                                                                                                     ccAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAg                                                         <  1:750506/239‑1 (MQ=255)
                                                                                                                                                                                      ctgTTTTCTTTTCTTTTCATTTTGCCCTCCCTTTCTGGGAGGTACCATCCAGTGGCGGTGTTCCCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGa                                                        <  2:343748/237‑1 (MQ=255)
                                                                                                                                                                                         ttttCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTTACAGTCCCAGTGAGGAATAAA‑CTTTATCCCTTAACGTCAGAGTGGTATCTCCGGTTTGTTCTTTAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCTTTCCCAAGGCTCCAGCCGTCCCCAGATGGCTGGTCCATTCCCTCTTGcgggc                                                       >  1:670782/1‑233 (MQ=255)
                                                                                                                                                                                         ttttCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCg                                                                                                                        <  1:571614/172‑1 (MQ=255)
                                                                                                                                                                                         ttttCTGATCCCTTCATTTTGCGCTCCTTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCCCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGaa                                                                                                    <  2:776876/192‑1 (MQ=255)
                                                                                                                                                                                                  ctctTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCg                                                                                                                        <  1:774867/163‑1 (MQ=255)
                                                                                                                                                                                                    ctTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGGCCATTTTGACCTGc                                          <  1:75053/239‑1 (MQ=255)
                                                                                                                                                                                                         tatTGCGCTCCCTTTCTGGGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGACCATTTTGACCTGCGGGTa                                     <  1:328941/239‑1 (MQ=255)
                                                                                                                                                                                                                 tCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTTCCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCg                                                                                                      <  2:152541/166‑1 (MQ=255)
                                                                                                                                                                                                                 tCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGGTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCATTAACGTCAGAGTGGTATTGCCGGGTTGTTCTATAAATTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTTCCTCTTGCGGACCCCTTTTACCTTCGGGGGGAGCGccc                             >  1:443496/1‑238 (MQ=255)
                                                                                                                                                                                                                 tCCCTTGCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGACCATTTTGACCTGCGGGTAGAGCGCCg                             <  1:460279/239‑1 (MQ=255)
                                                                                                                                                                                                                           gagggTATCATCCGGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCATTAACGTCAGAGTTTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGACCATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGc                   <  1:415310/236‑1 (MQ=255)
                                                                                                                                                                                                                                     ccccggggCGGTTTTCTCCCGCCCGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGACCATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGCATGTGTCTgg         <  2:904830/233‑1 (MQ=255)
                                                                                                                                                                                                                                           ggCGTTGTTCTCCCGCGAGGCTCCAGGCATTTAACAGTCCAAGTGAGTAATAAC‑CTTTATCCAGTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTAAGTTCCATCTGGCAGACCATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGCATGTGTCTGGCgggg    <  1:139433/238‑1 (MQ=255)
                                                                                                                                                                                                                                           ggCGGTGTTCCCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGACCATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGCATGTGTCTGGCgggg    <  1:843129/238‑1 (MQ=255)
                                                                                                                                                                                                                                             ctgtgtTTCCCCGTCAGGCTCCAGGCTTGTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGACCATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGCATGTGTCTGGCGGGGGa  <  1:918345/236‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                       
ACAGATCCCCCATTTCTCCTCCGGGGAAGCGCATGGCGTCCGTCATCGACGTCAGCAGATGGCAGACTTTTCC‑TTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAAC‑CTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGACCATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGCATGTGTCTGGCGGGGGA  >  NC_000913/4502513‑4502986

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