Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,888,212 Δ1 bp 7.6% intergenic (+20/‑223) mnmE → / → tnaC tRNA U34 5‑methylaminomethyl‑2‑thiouridine modification GTPase/tryptophanase leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,888,2120T. 7.6% 319.0 / 18.9 92intergenic (+20/‑223)mnmE/tnaCtRNA U34 5‑methylaminomethyl‑2‑thiouridine modification GTPase/tryptophanase leader peptide
Reads supporting (aligned to +/- strand):  ref base T (32/53);  new base . (3/4);  total (35/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.97e-03

CAAACAGAGCATGGGCTTTGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATTTAAACAATTTCAGAATAGACAAAAACTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCATTATGAATATCTTACATATATGTGTGACCTCAAAATGGTTCAATATTGAC  >  NC_000913/3887943‑3888482
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                               
cAAACAGTGCATGGGCTTTGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTg                                                                                                                                                                                                                                                                             <  2:414499/274‑1 (MQ=255)
  aaCAGGGCAATGGCTTGGACACAAACATGGGAGGCGGTTTCAGGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAAATTTGGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACt                                                                                                                                                                                                                                                                          <  4:478260/275‑1 (MQ=255)
   ccaagtattGGTCTTTGACCCCAACATGGAAGGCGGCGTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAGCAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGGCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACt                                                                                                                                                                                                                                                                          <  4:400297/266‑1 (MQ=255)
               ctttagaacCCACACGGTAGGGGGCTTCCTAGCGGGGCGTGCCGCCCTACAGGGGCTGGAACAGGCAGCGGAACATCNAAAGCGGGGCAAAGCGCAACTTTTGGGAGCCTGGGCAGGTGAACTGCTTGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTCAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGgct                                                                                                                                                                                                                                                            <  3:176409/268‑1 (MQ=255)
               cgtgaacacAAACAGGGAGGGGGTCTCCCNAGGGCGTCGTCGAAACCTAGAGGGGNGGGAACAGGCAGCGGAAAATGTACAACAGGGGAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGNTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTGCTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGGGCTTACGAAGCCGCATTCTGACTGTCAGATGCGgct                                                                                                                                                                                                                                                            <  3:137580/270‑1 (MQ=255)
                   gACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGc                                                                                                                                                                                                                                                                >  2:410128/1‑267 (MQ=255)
                       ccAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGgct                                                                                                                                                                                                                                                            <  4:410126/267‑1 (MQ=255)
                        cAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGCAGAGTTGCGTCTGGCACAGCCGAACATAAGCGAAATCTCCGGGGAATTTACTTCAGACGAACTGCTGGGGCGGATTTTCTCCTGCTTCTGTATTTGTAGGTAACCGCACTTACGAAGCCGCATTCTGACTGTAACATACGGCTTCCCTTCAt                                                                                                                                                                                                                                                   >  4:226061/1‑275 (MQ=255)
                              ggagcgcgggctacTGGCGCGTCGTCGAAACCGACAGGACACGGAACAGGCTGCGGGACATCTACAACAAGGCAAGGAGCTACTGTGTGAAGCCAGGGTAGGTGAAAGGCTAGCGGAGCAGTTGAGTCTGGAAAAGCAGAACTTAAGCGAAATCACAGGGGTATTGACTGCGGATGACCTGCAGGGGCGGATTTTATCCAGCTTCGGTAGTGGTAAGTAACCGCGCTTGCGAAGCCGTATTCTGACTGTCAGATGCGGCTTCGCTTCATGGTTAc                                                                                                                                                                                                                                             <  4:140655/262‑1 (MQ=255)
                                 aGGCGGGTTCCTGGCGGGGCGTGGCCACCTACCGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCAc                                                                                                                                                                                                                                          <  4:187116/275‑1 (MQ=255)
                                 aGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCAc                                                                                                                                                                                                                                          <  2:473112/275‑1 (MQ=255)
                                  ggcggcTTCCTGGGGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGGAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACt                                                                                                                                                                                                                                         <  2:385466/275‑1 (MQ=255)
                                       cTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTa                                                                                                                                                                                                                                              >  2:659786/1‑265 (MQ=255)
                                         ccccggCGGGTCGTCGCCCCCGCCAGGCGCTTGAACAGGCAGCGGAACATCTACACCAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGGATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTa                                                                                                                                                                                                                                  <  4:565978/271‑1 (MQ=255)
                                           cTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCAc                                                                                                                                                                                                                                          <  4:659784/265‑1 (MQ=255)
                                            tGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTAGCTTCATTTTTACCACTCCTGTTATTc                                                                                                                                                                                                                               >  1:139084/1‑275 (MQ=255)
                                                                       tGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTg                                                                                                                                                                                                                                                                             >  1:106159/1‑203 (MQ=255)
                                                                           aCAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTg                                                                                                                                                                                                                                                                         <  3:106158/203‑1 (MQ=255)
                                                                                 aGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGCTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGCGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTACCCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGGTTCATTGTTACCACTCCTGTTGTTCCTCAACCtttttt                                                                                                                                                                                                                  >  2:189314/1‑251 (MQ=255)
                                                                                            taaaacaGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTACAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGAGGCTTCGCTTCATTGTTACCAc                                                                                                                                                                                                                                          <  2:244074/213‑1 (MQ=255)
                                                                                                    ggCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCtt                                                                                                                                                                                                >  2:87925/1‑250 (MQ=255)
                                                                                                        aaGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGt                                                                                                                                                                                            <  4:87924/250‑1 (MQ=255)
                                                                                                              aaCTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCAtt                                                                                                                                                                                                                                                  >  1:104185/1‑190 (MQ=255)
                                                                                                                  gTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTa                                                                                                                                                                                                                                              <  3:104184/190‑1 (MQ=255)
                                                                                                                             gggCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGNGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTTTTATTCCTCANCCCTTTTTTTAAACATTAAAATTTTTTATTAATTTATAATTTTTTNANACAAAATTTAATATTTCTTCCCaaaaaa                                                                                                                                              >  3:82591/1‑272 (MQ=255)
                                                                                                                             gggCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTg                                                                                                                                                                                                                                                 >  2:515669/1‑176 (MQ=255)
                                                                                                                              ggCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTg                                                                                                                                                                                                                                                                             <  6:786798/148‑1 (MQ=255)
                                                                                                                              ggCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGa                                                                                                                                                                                                                                                                            <  7:654781/149‑1 (MQ=255)
                                                                                                                              ggCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATGTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTAc                                                                                                                                                                                                                                             <  2:539106/179‑1 (MQ=255)
                                                                                                                                 aGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTAc                                                                                                                                                                                                                                             <  4:515667/176‑1 (MQ=255)
                                                                                                                                   gTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTc                                                                                                                                                                                                                                                                       <  7:787038/149‑1 (MQ=255)
                                                                                                                                   gTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGgct                                                                                                                                                                                                                                                            >  4:773549/1‑159 (MQ=255)
                                                                                                                                     gAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAg                                                                                                                                                                                                                                                                     <  7:586213/149‑1 (MQ=255)
                                                                                                                                       aCTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGgcttcgc                                                                                                                                                                                                                                                        <  2:773551/159‑1 (MQ=255)
                                                                                                                                            tgGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATCTGACAGTCAGAATGCGgc                                                                                                                                                                                                                                                             <  7:107136/149‑1 (MQ=255)
                                                                                                                                              gCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGgct                                                                                                                                                                                                                                                            >  6:301240/1‑148 (MQ=255)
                                                                                                                                                ggAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGCATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATCTGACAGTCAGAATGCGgctcc                                                                                                                                                                                                                                                          <  7:747457/148‑3 (MQ=255)
                                                                                                                                                ggAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGgcttcg                                                                                                                                                                                                                                                         >  7:446869/1‑149 (MQ=255)
                                                                                                                                                ggAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATCTGACAGTCAGAATGCGgctcc                                                                                                                                                                                                                                                          >  8:747457/1‑146 (MQ=255)
                                                                                                                                                 gAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGAc                                                                                                                                                                                                                                                                           <  7:494920/131‑1 (MQ=255)
                                                                                                                                                 gAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGAc                                                                                                                                                                                                                                                                           >  8:494920/1‑131 (MQ=255)
                                                                                                                                                   agagTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATTTa                                                                                                                         <  1:82592/274‑1 (MQ=255)
                                                                                                                                                    gagTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGgcttcgcttc                                                                                                                                                                                                                                                     >  7:403391/1‑149 (MQ=255)
                                                                                                                                                    gagTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGgcttcgcttc                                                                                                                                                                                                                                                     <  7:699577/149‑1 (MQ=255)
                                                                                                                                                           gTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTAGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTa                                                                                                                                                                                                                                              <  7:289525/149‑1 (MQ=255)
                                                                                                                                                             cTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTAc                                                                                                                                                                                                                                             >  5:770807/1‑148 (MQ=255)
                                                                                                                                                               ggCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGNTGCGGCTTCGCTTCATTGTTACCNCTccn                                                                                                                                                                                                                                      >  4:722972/1‑152 (MQ=255)
                                                                                                                                                               ggCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCtttt                                                                                                                                                                                                                   >  4:73813/1‑172 (MQ=255)
                                                                                                                                                                  acaGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTc                                                                                                                                                                                                                                        <  6:556704/148‑1 (MQ=255)
                                                                                                                                                                   cagcagAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTa                                                                                                                                                                                                                                  <  2:722974/153‑1 (MQ=255)
                                                                                                                                                                   cagcagAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTa                                                                                                                                                                                                               <  2:73814/172‑1 (MQ=255)
                                                                                                                                                                       agAACTTAAGCGAAATCACCGGGGAATTTACTTCAGCCGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTa                                                                                                                                                                                                                                              <  8:81264/137‑1 (MQ=255)
                                                                                                                                                                       agAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTa                                                                                                                                                                                                                                              >  7:81264/1‑137 (MQ=255)
                                                                                                                                                                          aCTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTAtt                                                                                                                                                                                                                                <  5:289474/148‑1 (MQ=255)
                                                                                                                                                                             tAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATCTGACAGTCAGaatg                                                                                                                                                                                                                                                                <  6:872652/112‑4 (MQ=255)
                                                                                                                                                                             tAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATCTGACAGTCAGaatg                                                                                                                                                                                                                                                                >  5:872652/1‑109 (MQ=255)
                                                                                                                                                                                gCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAc                                                                                                                                                                                                                         <  7:428262/149‑1 (MQ=255)
                                                                                                                                                                                 cGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAc                                                                                                                                                                                                                         >  5:898433/1‑148 (MQ=255)
                                                                                                                                                                                      tCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCttt                                                                                                                                                                                                                    <  5:749890/148‑1 (MQ=255)
                                                                                                                                                                                      tCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCttt                                                                                                                                                                                                                    <  6:770807/148‑1 (MQ=255)
                                                                                                                                                                                         ccGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCttttttt                                                                                                                                                                                                                >  7:700488/1‑149 (MQ=255)
                                                                                                                                                                                            gggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTaa                                                                                                                                                                                                              <  5:301240/148‑1 (MQ=255)
                                                                                                                                                                                                  ttACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTa                                                                                                                                                                                                        <  6:779999/148‑1 (MQ=255)
                                                                                                                                                                                                      ttCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAAtt                                                                                                                                                                                                   >  7:245155/1‑149 (MQ=255)
                                                                                                                                                                                                      ttCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCtt                                                                                                                                                                                                >  3:74473/1‑152 (MQ=255)
                                                                                                                                                                                                       ccAGACGACCTGCTGGGGCGGATTTTCTTCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACCTTAAAATTc                                                                                                                                                                                                  <  8:779670/148‑1 (MQ=255)
                                                                                                                                                                                                          gacgacCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGt                                                                                                                                                                                            <  1:74474/152‑1 (MQ=255)
                                                                                                                                                                                                                    tgGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTtat                                                                                                                                                                                     <  8:446869/149‑1 (MQ=255)
                                                                                                                                                                                                                         cGGTTTTTTTCCAGGTTTTGTATTTTGAANTATCCTCATTTACAAAACCGAATTANGACTGTCAGATGCTGCTTCGCTTCATTGTTACCTCTCCTGTTATTTCTCAACCCTTTTTTTATACATTATAAATCTTACTNTATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATGTAAACAATTTCAGAATAGACAAAAACTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCAtt                                                  <  3:139083/275‑1 (MQ=255)
                                                                                                                                                                                                                          ggATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCtt                                                                                                                                                                               <  8:403391/149‑1 (MQ=255)
                                                                                                                                                                                                                              tttCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTaa                                                                                                                                                                            >  6:117227/1‑148 (MQ=255)
                                                                                                                                                                                                                                  tcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTaaaaaaa                                                                                                                                                                       <  8:178024/149‑1 (MQ=255)
                                                                                                                                                                                                                                   ccAGCTTCTGTAGTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAGTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATTTAAACAATTTCAGAATAGACAAAAACTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCATTATGAATATc                                         <  2:313123/274‑1 (MQ=255)
                                                                                                                                                                                                                                           tGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTaa                                                                                                                                                               <  6:898433/148‑1 (MQ=255)
                                                                                                                                                                                                                                            gTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAata                                                                                                                                                             <  8:847879/149‑1 (MQ=255)
                                                                                                                                                                                                                                            gTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAata                                                                                                                                                             <  8:700488/149‑1 (MQ=255)
                                                                                                                                                                                                                                               ttGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTg                                                                                                                                                          <  7:648139/149‑1 (MQ=255)
                                                                                                                                                                                                                                                 ggTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGAGTCACATTTAAACAATTTCAGAATAGACAAAAACTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGCGGCTAAGTGCATTATGAATATCTTACATATAtgtgtg                          >  1:154770/1‑275 (MQ=255)
                                                                                                                                                                                                                                                   taagtaaCCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCTTTTTTTTAAACATTaaa                                                                                                                                                                                                      <  6:669369/101‑1 (MQ=255)
                                                                                                                                                                                                                                                   taagtaaCCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCTTTTTTTTAAACATTaaa                                                                                                                                                                                                      >  5:669369/1‑101 (MQ=255)
                                                                                                                                                                                                                                                     agtaaCCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAAc                                                                                                                                                                                                            >  2:762982/1‑93 (MQ=255)
                                                                                                                                                                                                                                                       taaCCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTg                                                                                                                                                                                                                                                 >  2:59565/1‑54 (MQ=255)
                                                                                                                                                                                                                                                       taaCCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTAt                                                                                                                                                                                                                                 >  4:371982/1‑70 (MQ=255)
                                                                                                                                                                                                                                                         aCCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTa                                                                                                                                                                                                        <  4:762980/93‑1 (MQ=255)
                                                                                                                                                                                                                                                           cgcgCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTAc                                                                                                                                                                                                                                             <  4:59564/54‑1 (MQ=255)
                                                                                                                                                                                                                                                           cgcgCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCt                                                                                                                                                                                                                             <  2:371984/70‑1 (MQ=255)
                                                                                                                                                                                                                                                            gcgcTTACGAAGCCGCATCCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAACATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATTTAAACAATTTCAGAATAGATAAAAACTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCAtt                                                  >  1:11052/1‑240 (MQ=255)
                                                                                                                                                                                                                                                                  aCGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAAc                                                                                                                                                <  8:879317/140‑1 (MQ=255)
                                                                                                                                                                                                                                                                  aCGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAAc                                                                                                                                                >  7:879317/1‑140 (MQ=255)
                                                                                                                                                                                                                                                                       gCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTgat                                                                                                                                       <  5:765557/144‑1 (MQ=255)
                                                                                                                                                                                                                                                                       gCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTgat                                                                                                                                       >  6:765557/1‑144 (MQ=255)
                                                                                                                                                                                                                                                                          gCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAGACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATTTAAACAATTTCAGAATAGACAAAAACTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCATTATGAATATCTTACATATATGTTTGACCTCAAAATGGTTCAATATTGac  >  2:156635/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                               
CAAACAGAGCATGGGCTTTGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATTTAAACAATTTCAGAATAGACAAAAACTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCATTATGAATATCTTACATATATGTGTGACCTCAAAATGGTTCAATATTGAC  >  NC_000913/3887943‑3888482

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: