Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,888,219 T→A 10.1% intergenic (+27/‑216) 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,2190TA10.1% 285.6 / 19.3 89intergenic (+27/‑216)mnmE/tnaCtRNA U34 5‑methylaminomethyl‑2‑thiouridine modification GTPase/tryptophanase leader peptide
Reads supporting (aligned to +/- strand):  ref base T (31/49);  new base A (4/5);  total (35/54)
Fisher's exact test for biased strand distribution p-value = 7.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.95e-02

AACAGAGCATGGGCTTTGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATTTAAACAATTTCAGAATAGACAAAAACTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCATTATGAATATCTTACATATATGTGTGACCTCAAAATGGTTCAATATTGAC  >  NC_000913/3887945‑3888482
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                        
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)
                                                                         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)
                                                                                                                            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)
                                                                                                                                                 agagTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATTTa                                                                                                                         <  1:82592/274‑1 (MQ=255)
                                                                                                                                                  gagTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGgcttcgcttc                                                                                                                                                                                                                                                     <  7:699577/149‑1 (MQ=255)
                                                                                                                                                  gagTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGgcttcgcttc                                                                                                                                                                                                                                                     >  7:403391/1‑149 (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                                                                                                                                                                                                                    <  6:770807/148‑1 (MQ=255)
                                                                                                                                                                                    tCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCttt                                                                                                                                                                                                                    <  5:749890/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:700488/149‑1 (MQ=255)
                                                                                                                                                                                                                                          gTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAata                                                                                                                                                             <  8:847879/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                                                                                                                                                >  7:879317/1‑140 (MQ=255)
                                                                                                                                                                                                                                                                aCGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAAc                                                                                                                                                <  8:879317/140‑1 (MQ=255)
                                                                                                                                                                                                                                                                     gCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTgat                                                                                                                                       >  6:765557/1‑144 (MQ=255)
                                                                                                                                                                                                                                                                     gCCGCATCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTgat                                                                                                                                       <  5:765557/144‑1 (MQ=255)
                                                                                                                                                                                                                                                                        gCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAGACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATTTAAACAATTTCAGAATAGACAAAAACTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCATTATGAATATCTTACATATATGTTTGACCTCAAAATGGTTCAATATTGac  >  2:156635/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                          catCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTa                                                                                                                                                                             >  6:204973/3‑102 (MQ=39)
                                                                                                                                                                                                                                                                          catCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTa                                                                                                                                                                             <  5:204973/100‑1 (MQ=39)
                                                                                                                                                                                                                                                                                aCTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATTTaaa                                                                                                                       <  8:245155/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                                aCTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATTTaaa                                                                                                                       <  8:373290/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                                aCTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATTTaaa                                                                                                                       >  8:133194/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                        
AACAGAGCATGGGCTTTGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATTTAAACAATTTCAGAATAGACAAAAACTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCATTATGAATATCTTACATATATGTGTGACCTCAAAATGGTTCAATATTGAC  >  NC_000913/3887945‑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: