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

ACAGATCCCCCATTTCTCCTCCGGGGAAGCGCATGGCGTCCGTCATCGACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGACCATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGCATGTGTCTGGCGGGGGA  >  NC_000913/4502513‑4502986
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                       
ccacaTCCCCCTTTTCCCCTCCGGGGAAGCGCAGGGGGCCCGTCACCGACGCCAGCAGTTGGAAGCTTTTTCCTTTGCCGGTCCCGGTCACCAGGCGGGCAGTCCCTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACC‑‑‑‑‑‑‑CTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTtctc                                                                                                                                                                                                                                       <  2:610719/235‑1 (MQ=255)
       ccccATTTCTCCTCCGGGGAAGCGCATGGCGTCCGTCATCGACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCTTGCCTGCGGGCTGGTTTTCAGCTTACCCCGATGATCGCCTTTACCCAGTTTTCTGGTCTCTTCATATTGCCCTCCCTTCCTCGGGGGATTCATCCAGTGGCGGGTTtct                                                                                                                                                                                                                                        >  2:81351/1‑238 (MQ=255)
       ccccATTTCCCCACCGGGGAAGCGCAGGGCGTCCGTCATGGACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTtctc                                                                                                                                                                                                                                       <  1:292430/239‑1 (MQ=255)
         ccATTTCTCCTCCGGGGAAGCGCATGGCGTCCGTCATCGACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTTTcc                                                                                                                                                                                                                                      >  1:117236/1‑238 (MQ=255)
         ccATTTCTCCTCCGGGGAAGCGCATGGCGTCCGTCATCGACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACC‑‑‑‑‑‑‑CTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGGGGTTTTCCGCCCgg                                                                                                                                                                                                                               >  1:610719/1‑238 (MQ=255)
             ttCCCCCCCGGGGAAGCCCATGGCGGCCGTCAGCGGCGTCAGTAGATGGCAGTCTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCa                                                                                                                                                                                                                                 <  1:724208/239‑1 (MQ=255)
              tctcCTCCGGGGAAGCGCATGGCGTCCGTCATCGACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGCTGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTTCGCTCCCTTTCTGAGAGGTATCATCCAGTGGGGGCGGTTTCTCGcccc                                                                                                                                                                                                                                >  2:750506/1‑237 (MQ=255)
                   tccGGGGAAGCGAGTGGCGTACGTCGACGACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTAACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCCCTTTTTGCGGGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTcc                                                                                                                                                                                                                           <  1:628689/239‑1 (MQ=255)
                   cccGGGGAAGCGCATGGCCGCCGTCATCGCCGGCATCAGATGGCAGACTTTTCCTTTGCGCGTCCCGGTCAGCGGGCGGGCAGGCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTcc                                                                                                                                                                                                                           <  2:824813/238‑1 (MQ=255)
                             cccctGGCGTCCGTTCTCGGCGCCAGCAGATGGCCGCTTTTTCCTTTGCGCGTCCCGGTCCGCCGGCGGGCAGTCATTTCATTTCCCAGTCCCGGCCCCTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTaa                                                                                                                                                                                                                 <  2:395462/235‑1 (MQ=255)
                              gcgtGGCGCCGGTCTTCGTCGTCAGCAGATGGCAGACTGTGCCTTTGTGCGTCACGGTCAGCTTGCGGGCAGTCACTTCATTCCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGTGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTaa                                                                                                                                                                                                                 <  1:351081/235‑1 (MQ=255)
                              ccATGGGGTCCGGCATCGACGTTCGCAGTTTGCAAACTTTTCCGTGGCGCGTCACGGTCAGGAGGCGGGGCGTCCCTTCATTTCCCAGTCCCGGCCACTTTTTTCGTCCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTTTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTaa                                                                                                                                                                                                                 <  2:30030/237‑1 (MQ=255)
                                      gccgtcAACGGCGTCAGGAGATGGCAGATTTTTCCTTCGCGCCGCACGGGCAGCAGGCGGGCAGTCACTTCATTTCCTAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAg                                                                                                                                                                                                        <  1:61481/238‑1 (MQ=255)
                                       ccgtcATCGACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCTCTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGt                                                                                                                                                                                                       <  1:340211/239‑1 (MQ=255)
                                        cgtcATCGACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGGATAATCCAGTGGCGGTGTTTTCCCCCCCGGTCCCCGGCCATTTACAGTGCCCGTg                                                                                                                                                                                                      >  1:855322/1‑239 (MQ=255)
                                        cgccATCGACGTCAGCAGTGTGCAGCCTTTTCCTTTGGGCGTTCCGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCCCCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGt                                                                                                                                                                                                       <  1:582100/235‑1 (MQ=255)
                                              cggcgtcgtcagATGGCATAGTTGTCTTTTCGGTCTCACTGTCCTCAGGCGGGCAGTCACGTCATTTCCCAGTCCCGGCCCCTTTTTTCGTGGCTGTGGGCTGGTTTTCAGCTTCACCAGGTGATCGCTTTGGCCAAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCATGCATTTAACAGTCCCAGTGAGtaac                                                                                                                                                                                                <  2:151120/230‑2 (MQ=255)
                                              cggcGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCCGGCGGGCAGTCCCTTAATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGtaat                                                                                                                                                                                                <  1:889627/236‑1 (MQ=255)
                                                        agATGGCAGACTTTTCCTTTGCGCGTCCCGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTAt                                                                                                                                                                                      <  1:16262/238‑1 (MQ=255)
                                                             ctcgttttttCCTTTTGCGGTCCCGGTCCACGGGCGGGCATTCCCTTCTTTTCCCAGCCCCGGCCACTTTTTTCGGGCCTGCGGGCTGGTTTTCAGCTTCCCCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTa                                                                                                                                                                                <  2:261719/233‑1 (MQ=255)
                                                                   ttttCCTGTGCGCGTCGCGGTCAGCAGGGGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTGCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCa                                                                                                                                                                          <  1:243073/239‑1 (MQ=255)
                                                                      tCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTTTCCCCCCCGGCTCCAGGCACTTTACAGTCCCCGTGAGTAATAAACTTTATCCATTTACGTCataa                                                                                                                                                                       >  1:676470/1‑236 (MQ=255)
                                                                                   cACGGTCAGCGGGCGGGCAGTCACTCCATCCCCAGTCCCGGCCACTTTTTCCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGt                                                                                                                                                          <  1:506043/239‑1 (MQ=255)
                                                                                      ggTCCGCAGGCGGGCAGTCCCTTCATTTCCCAGTCCCGGCCTCTTTTTTCGTGCCTGCGGGCGGGTTTTCAGTTTCCCCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGttt                                                                                                                                                        <  2:219754/239‑1 (MQ=255)
                                                                                        tCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTTTCCCGCCCAGCCCCACGCCTTTTACAATACCAGTGGGTAATAACACTTTTTCCCTTACGTCAAAGGGTTTTCGCTGGtcctt                                                                                                                                                      >  1:18080/1‑235 (MQ=255)
                                                                                          agcagGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTaa                                                                                                                                                                                                                 <  1:914524/179‑1 (MQ=255)
                                                                                           gaagggggtaaGTCTCTTTATTTCCCAGTCCCGGCCACTTTTTGCGTGCCTGCGGGCTGGTTTGCAGCTTCACCAGATGATCGCCTTTACCCAGTTTGCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCt                                                                                                                                                   <  2:445415/229‑1 (MQ=255)
                                                                                           ccgggggggCATTCTCTTCCTTTCCCAGCCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCCCCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCt                                                                                                                                                   <  2:676470/233‑1 (MQ=255)
                                                                                            cagGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTTTCCCGCCACGCTCCCGGCATTTTAAAGTACCCGTGAGTGAAAAACTTTTTCCCATAACCTCAGAAGGTTTTCTCCggggttttttt                                                                                                                                                  >  1:722299/1‑230 (MQ=255)
                                                                                                         aCTTCATTTCCCAGGCCCGGCAACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCGGAGCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGc                                                                                                                                     <  1:81351/239‑1 (MQ=255)
                                                                                                           ttCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCAttt                                                                                                                                                                                                                   <  1:919557/160‑1 (MQ=255)
                                                                                                                      ggcccGGCCCTTTTTGGCGTGCCTTGGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCg                                                                                                                        <  1:348699/237‑1 (MQ=255)
                                                                                                                                               cGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTTTCCCGCCACGCTCCCGGCCTTTTAAAATCCCCGGGAGTAATAACCCTTATCCCTTAACGTCAGAGTGTTTTCGCCGGGTTGTTCTATAAATTTCTCAGCAAACTCATTTTCTCTGTTCTTCATCGTGCCGAAGGCt                                                                                                >  1:904830/1‑238 (MQ=255)
                                                                                                                                                  gCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCa                                                                                             <  1:873907/238‑1 (MQ=255)
                                                                                                                                                              cTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCg                                                                                                                        <  1:300333/199‑1 (MQ=255)
                                                                                                                                                                   ccAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTg                                                                           <  1:169062/239‑1 (MQ=255)
                                                                                                                                                                       atgatCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCaga                                                                                                                                                                        <  1:384994/142‑1 (MQ=255)
                                                                                                                                                                                     ccAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAg                                                         <  1:750506/239‑1 (MQ=255)
                                                                                                                                                                                      ctgTTTTCTTTTCTTTTCATTTTGCCCTCCCTTTCTGGGAGGTACCATCCAGTGGCGGTGTTCCCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGa                                                        <  2:343748/237‑1 (MQ=255)
                                                                                                                                                                                         ttttCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTTACAGTCCCAGTGAGGAATAAACTTTATCCCTTAACGTCAGAGTGGTATCTCCGGTTTGTTCTTTAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCTTTCCCAAGGCTCCAGCCGTCCCCAGATGGCTGGTCCATTCCCTCTTGcgggc                                                       >  1:670782/1‑233 (MQ=255)
                                                                                                                                                                                         ttttCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCg                                                                                                                        <  1:571614/172‑1 (MQ=255)
                                                                                                                                                                                         ttttCTGATCCCTTCATTTTGCGCTCCTTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCCCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGaa                                                                                                    <  2:776876/192‑1 (MQ=255)
                                                                                                                                                                                                  ctctTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCg                                                                                                                        <  1:774867/163‑1 (MQ=255)
                                                                                                                                                                                                    ctTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGGCCATTTTGACCTGc                                          <  1:75053/239‑1 (MQ=255)
                                                                                                                                                                                                         tatTGCGCTCCCTTTCTGGGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGACCATTTTGACCTGCGGGTa                                     <  1:328941/239‑1 (MQ=255)
                                                                                                                                                                                                                 tCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTTCCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCg                                                                                                      <  2:152541/166‑1 (MQ=255)
                                                                                                                                                                                                                 tCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGGTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGGTATTGCCGGGTTGTTCTATAAATTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTTCCTCTTGCGGACCCCTTTTACCTTCGGGGGGAGCGccc                             >  1:443496/1‑238 (MQ=255)
                                                                                                                                                                                                                 tCCCTTGCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGACCATTTTGACCTGCGGGTAGAGCGCCg                             <  1:460279/239‑1 (MQ=255)
                                                                                                                                                                                                                           gagggTATCATCCGGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTTTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGACCATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGc                   <  1:415310/236‑1 (MQ=255)
                                                                                                                                                                                                                                     ccccggggCGGTTTTCTCCCGCCCGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGACCATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGCATGTGTCTgg         <  2:904830/233‑1 (MQ=255)
                                                                                                                                                                                                                                           ggCGTTGTTCTCCCGCGAGGCTCCAGGCATTTAACAGTCCAAGTGAGTAATAACCTTTATCCAGTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTAAGTTCCATCTGGCAGACCATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGCATGTGTCTGGCgggg    <  1:139433/238‑1 (MQ=255)
                                                                                                                                                                                                                                           ggCGGTGTTCCCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGACCATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGCATGTGTCTGGCgggg    <  1:843129/238‑1 (MQ=255)
                                                                                                                                                                                                                                             ctgtgtTTCCCCGTCAGGCTCCAGGCTTGTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGACCATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGCATGTGTCTGGCGGGGGa  <  1:918345/236‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                       
ACAGATCCCCCATTTCTCCTCCGGGGAAGCGCATGGCGTCCGTCATCGACGTCAGCAGATGGCAGACTTTTCCTTTGCGCGTCACGGTCAGCAGGCGGGCAGTCACTTCATTTCCCAGTCCCGGCCACTTTTTTCGTGCCTGCGGGCTGGTTTTCAGCTTCACCAGATGATCGCCTTTACCCAGTTTTCTGATCTCTTCATATTGCGCTCCCTTTCTGAGAGGTATCATCCAGTGGCGGTGTTCTCCCGCCAGGCTCCAGGCATTTAACAGTCCCAGTGAGTAATAACCTTTATCCATTAACGTCAGAGTGTTATCGCCGGTTTGTTCTATAAGTTGCTCAGCAAGCTCATTTTCGCTGTTCTTCATCGTGCCGAAGGCTGCAGCCGTCAGCAGATGGCTGGTCAGTTCCATCTGGCAGACCATTTTGACCTGCGGGTAGAGCGCCGGGTTCCCGGCATGTGTCTGGCGGGGGA  >  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: