Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,888,224:1 +A 9.5% intergenic (+32/‑211) 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,2241.A 9.5% 302.0 / 33.2 84intergenic (+32/‑211)mnmE/tnaCtRNA U34 5‑methylaminomethyl‑2‑thiouridine modification GTPase/tryptophanase leader peptide
Reads supporting (aligned to +/- strand):  ref base . (30/46);  new base A (4/4);  total (34/50)
Fisher's exact test for biased strand distribution p-value = 7.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.64e-02

CTTTGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATTTAAACAATTTCAGAATAGACAAAAACTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCATTATGAATATCTTACATATATGTGTGACCTCAAAATGGTTCAATATTGAC  >  NC_000913/3887958‑3888482
                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                  
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)
                                                                  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)
                                                                                                                    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: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: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:133194/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                   aCTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATTTaaa                                                                                                                       <  8:245155/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                   aCTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATTTaaa                                                                                                                       <  8:373290/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                    cagTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATTTaa                                                                                                                        >  5:911911/3‑148 (MQ=255)
                                                                                                                                                                                                                                                                          gATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATCATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATTTAAACAATTTc                                                                                                                >  8:537170/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                  
CTTTGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAGATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATTTAAACAATTTCAGAATAGACAAAAACTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCATTATGAATATCTTACATATATGTGTGACCTCAAAATGGTTCAATATTGAC  >  NC_000913/3887958‑3888482

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