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

AACAGAGCATGGGCTTTGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATTTAAACAATTTCAGAATAGACAAAAACTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCATTATGAATATCTTACATATATGTGTGACCTCAAAATGGTTCAATATTGAC  >  NC_000913/3887945‑3888482
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                        
aaCAGGGCAATGGCTTGGACACAAACATGGGAGGCGGTTTCAGGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAAATTTGGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACt                                                                                                                                                                                                                                                                          <  4:478260/275‑1 (MQ=255)
 ccaagtattGGTCTTTGACCCCAACATGGAAGGCGGCGTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAGCAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGGCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACt                                                                                                                                                                                                                                                                          <  4:400297/266‑1 (MQ=255)
             ctttagaacCCACACGGTAGGGGGCTTCCTAGCGGGGCGTGCCGCCCTACAGGGGCTGGAACAGGCAGCGGAACATCNAAAGCGGGGCAAAGCGCAACTTTTGGGAGCCTGGGCAGGTGAACTGCTTGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTCAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgct                                                                                                                                                                                                                                                            <  3:176409/268‑1 (MQ=255)
             cgtgaacacAAACAGGGAGGGGGTCTCCCNAGGGCGTCGTCGAAACCTAGAGGGGNGGGAACAGGCAGCGGAAAATGTACAACAGGGGAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGNTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTGCTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGGGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgct                                                                                                                                                                                                                                                            <  3:137580/270‑1 (MQ=255)
                 gACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGc                                                                                                                                                                                                                                                                >  2:410128/1‑267 (MQ=255)
                     ccAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgct                                                                                                                                                                                                                                                            <  4:410126/267‑1 (MQ=255)
                      cAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGCAGAGTTGCGTCTGGCACAGCCGAACATAAGCGAAATCTCCGGGGAATTTACTTCAGACGAACTGCTGGGGCGGATTTTCTCCTGCTTCTGTATTTGTAGGTAACCGCACTTACGAAGCCGCATTCTGACTGTAAC‑ATACGGCTTCCCTTCAt                                                                                                                                                                                                                                                   >  4:226061/1‑275 (MQ=255)
                            ggagcgcgggctacTGGCGCGTCGTCGAAACCGACAGGACACGGAACAGGCTGCGGGACATCTACAACAAGGCAAGGAGCTACTGTGTGAAGCCAGGGTAGGTGAAAGGCTAGCGGAGCAGTTGAGTCTGGAAAAGCAGAACTTAAGCGAAATCACAGGGGTATTGACTGCGGATGACCTGCAGGGGCGGATTTTATCCAGCTTCGGTAGTGGTAAGTAACCGCGCTTGCGAAGCCGTATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATGGTTAc                                                                                                                                                                                                                                             <  4:140655/262‑1 (MQ=255)
                               aGGCGGGTTCCTGGCGGGGCGTGGCCACCTACCGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCAc                                                                                                                                                                                                                                          <  4:187116/275‑1 (MQ=255)
                               aGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCAc                                                                                                                                                                                                                                          <  2:473112/275‑1 (MQ=255)
                                ggcggcTTCCTGGGGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGGAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACt                                                                                                                                                                                                                                         <  2:385466/275‑1 (MQ=255)
                                     cTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTa                                                                                                                                                                                                                                              >  2:659786/1‑265 (MQ=255)
                                       ccccggCGGGTCGTCGCCCCCGCCAGGCGCTTGAACAGGCAGCGGAACATCTACACCAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGGATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTa                                                                                                                                                                                                                                  <  4:565978/271‑1 (MQ=255)
                                         cTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCAc                                                                                                                                                                                                                                          <  4:659784/265‑1 (MQ=255)
                                          tGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTAGCTTCATTTTTACCACTCCTGTTATTc                                                                                                                                                                                                                               >  1:139084/1‑275 (MQ=255)
                                                                         aCAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTg                                                                                                                                                                                                                                                                         <  3:106158/203‑1 (MQ=255)
                                                                               aGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGCTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGCGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTACCCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGGTTCATTGTTACCACTCCTGTTGTTCCTCAACCtttttt                                                                                                                                                                                                                  >  2:189314/1‑251 (MQ=255)
                                                                                          taaaacaGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTACAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGAGGCTTCGCTTCATTGTTACCAc                                                                                                                                                                                                                                          <  2:244074/213‑1 (MQ=255)
                                                                                                  ggCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCtt                                                                                                                                                                                                >  2:87925/1‑250 (MQ=255)
                                                                                                      aaGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGt                                                                                                                                                                                            <  4:87924/250‑1 (MQ=255)
                                                                                                            aaCTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCAtt                                                                                                                                                                                                                                                  >  1:104185/1‑190 (MQ=255)
                                                                                                                gTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTa                                                                                                                                                                                                                                              <  3:104184/190‑1 (MQ=255)
                                                                                                                           gggCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGNGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTTTTATTCCTCANCCCTTTTTTTAAACATTAAAATTTTTTATTAATTTATAATTTTTTNANACAAAATTTAATATTTCTTCCCaaaaaa                                                                                                                                              >  3:82591/1‑272 (MQ=255)
                                                                                                                           gggCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTg                                                                                                                                                                                                                                                 >  2:515669/1‑176 (MQ=255)
                                                                                                                            ggCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATGTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTAc                                                                                                                                                                                                                                             <  2:539106/179‑1 (MQ=255)
                                                                                                                               aGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTAc                                                                                                                                                                                                                                             <  4:515667/176‑1 (MQ=255)
                                                                                                                                 gTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTc                                                                                                                                                                                                                                                                       <  7:787038/149‑1 (MQ=255)
                                                                                                                                 gTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgct                                                                                                                                                                                                                                                            >  4:773549/1‑159 (MQ=255)
                                                                                                                                   gAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAg                                                                                                                                                                                                                                                                     <  7:586213/149‑1 (MQ=255)
                                                                                                                                     aCTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgcttcgc                                                                                                                                                                                                                                                        <  2:773551/159‑1 (MQ=255)
                                                                                                                                          tgGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCA‑TCTGACAGTCAGAATGCGgc                                                                                                                                                                                                                                                             <  7:107136/149‑1 (MQ=255)
                                                                                                                                            gCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgct                                                                                                                                                                                                                                                            >  6:301240/1‑148 (MQ=255)
                                                                                                                                              ggAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGCATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCA‑TCTGACAGTCAGAATGCGgctcc                                                                                                                                                                                                                                                          <  7:747457/148‑3 (MQ=255)
                                                                                                                                              ggAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgcttcg                                                                                                                                                                                                                                                         >  7:446869/1‑149 (MQ=255)
                                                                                                                                              ggAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCA‑TCTGACAGTCAGAATGCGgctcc                                                                                                                                                                                                                                                          >  8:747457/1‑146 (MQ=255)
                                                                                                                                                 agagTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATTTa                                                                                                                         <  1:82592/274‑1 (MQ=255)
                                                                                                                                                  gagTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgcttcgcttc                                                                                                                                                                                                                                                     <  7:699577/149‑1 (MQ=255)
                                                                                                                                                  gagTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgcttcgcttc                                                                                                                                                                                                                                                     >  7:403391/1‑149 (MQ=255)
                                                                                                                                                         gTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTAGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTa                                                                                                                                                                                                                                              <  7:289525/149‑1 (MQ=255)
                                                                                                                                                           cTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTAc                                                                                                                                                                                                                                             >  5:770807/1‑148 (MQ=255)
                                                                                                                                                             ggCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑NTGCGGCTTCGCTTCATTGTTACCNCTccn                                                                                                                                                                                                                                      >  4:722972/1‑152 (MQ=255)
                                                                                                                                                             ggCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCtttt                                                                                                                                                                                                                   >  4:73813/1‑172 (MQ=255)
                                                                                                                                                                acaGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTc                                                                                                                                                                                                                                        <  6:556704/148‑1 (MQ=255)
                                                                                                                                                                 cagcagAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTa                                                                                                                                                                                                                                  <  2:722974/153‑1 (MQ=255)
                                                                                                                                                                 cagcagAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTa                                                                                                                                                                                                               <  2:73814/172‑1 (MQ=255)
                                                                                                                                                                     agAACTTAAGCGAAATCACCGGGGAATTTACTTCAGCCGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTa                                                                                                                                                                                                                                              <  8:81264/137‑1 (MQ=255)
                                                                                                                                                                     agAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTa                                                                                                                                                                                                                                              >  7:81264/1‑137 (MQ=255)
                                                                                                                                                                        aCTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTAtt                                                                                                                                                                                                                                <  5:289474/148‑1 (MQ=255)
                                                                                                                                                                           tAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCA‑TCTGACAGTCAG‑aatg                                                                                                                                                                                                                                                                <  6:872652/112‑4 (MQ=255)
                                                                                                                                                                           tAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCA‑TCTGACAGTCAG‑aatg                                                                                                                                                                                                                                                                >  5:872652/1‑109 (MQ=255)
                                                                                                                                                                              gCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAc                                                                                                                                                                                                                         <  7:428262/149‑1 (MQ=255)
                                                                                                                                                                               cGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAc                                                                                                                                                                                                                         >  5:898433/1‑148 (MQ=255)
                                                                                                                                                                                    tCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCttt                                                                                                                                                                                                                    <  6:770807/148‑1 (MQ=255)
                                                                                                                                                                                    tCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCttt                                                                                                                                                                                                                    <  5:749890/148‑1 (MQ=255)
                                                                                                                                                                                       ccGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCttttttt                                                                                                                                                                                                                >  7:700488/1‑149 (MQ=255)
                                                                                                                                                                                          gggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTaa                                                                                                                                                                                                              <  5:301240/148‑1 (MQ=255)
                                                                                                                                                                                                ttACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTa                                                                                                                                                                                                        <  6:779999/148‑1 (MQ=255)
                                                                                                                                                                                                    ttCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAAtt                                                                                                                                                                                                   >  7:245155/1‑149 (MQ=255)
                                                                                                                                                                                                    ttCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCtt                                                                                                                                                                                                >  3:74473/1‑152 (MQ=255)
                                                                                                                                                                                                     ccAGACGACCTGCTGGGGCGGATTTTCTTCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACCTTAAAATTc                                                                                                                                                                                                  <  8:779670/148‑1 (MQ=255)
                                                                                                                                                                                                        gacgacCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGt                                                                                                                                                                                            <  1:74474/152‑1 (MQ=255)
                                                                                                                                                                                                                  tgGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTtat                                                                                                                                                                                     <  8:446869/149‑1 (MQ=255)
                                                                                                                                                                                                                       cGGTTTTTTTCCAGGTTTTGTATTTTGAANTATCCTCATTTACAAAACCGAATTANGACTGTCAG‑ATGCTGCTTCGCTTCATTGTTACCTCTCCTGTTATTTCTCAACCCTTTTTTTATACATTATAAATCTTACTNTATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATGTAAACAATTTCAGAATAGACAAAAACTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCAtt                                                  <  3:139083/275‑1 (MQ=255)
                                                                                                                                                                                                                        ggATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCtt                                                                                                                                                                               <  8:403391/149‑1 (MQ=255)
                                                                                                                                                                                                                            tttCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTaa                                                                                                                                                                            >  6:117227/1‑148 (MQ=255)
                                                                                                                                                                                                                                tcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTaaaaaaa                                                                                                                                                                       <  8:178024/149‑1 (MQ=255)
                                                                                                                                                                                                                                 ccAGCTTCTGTAGTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAAGTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATTTAAACAATTTCAGAATAGACAAAAACTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCATTATGAATATc                                         <  2:313123/274‑1 (MQ=255)
                                                                                                                                                                                                                                         tGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTaa                                                                                                                                                               <  6:898433/148‑1 (MQ=255)
                                                                                                                                                                                                                                          gTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAata                                                                                                                                                             <  8:700488/149‑1 (MQ=255)
                                                                                                                                                                                                                                          gTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAata                                                                                                                                                             <  8:847879/149‑1 (MQ=255)
                                                                                                                                                                                                                                             ttGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTg                                                                                                                                                          <  7:648139/149‑1 (MQ=255)
                                                                                                                                                                                                                                               ggTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGAGTCACATTTAAACAATTTCAGAATAGACAAAAACTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGCGGCTAAGTGCATTATGAATATCTTACATATAtgtgtg                          >  1:154770/1‑275 (MQ=255)
                                                                                                                                                                                                                                                 taagtaaCCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCTTTTTTTTAAACATTaaa                                                                                                                                                                                                      <  6:669369/101‑1 (MQ=255)
                                                                                                                                                                                                                                                 taagtaaCCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCTTTTTTTTAAACATTaaa                                                                                                                                                                                                      >  5:669369/1‑101 (MQ=255)
                                                                                                                                                                                                                                                   agtaaCCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAAc                                                                                                                                                                                                            >  2:762982/1‑93 (MQ=255)
                                                                                                                                                                                                                                                     taaCCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTg                                                                                                                                                                                                                                                 >  2:59565/1‑54 (MQ=255)
                                                                                                                                                                                                                                                     taaCCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTAt                                                                                                                                                                                                                                 >  4:371982/1‑70 (MQ=255)
                                                                                                                                                                                                                                                       aCCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTa                                                                                                                                                                                                        <  4:762980/93‑1 (MQ=255)
                                                                                                                                                                                                                                                         cgcgCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTAc                                                                                                                                                                                                                                             <  4:59564/54‑1 (MQ=255)
                                                                                                                                                                                                                                                         cgcgCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCt                                                                                                                                                                                                                             <  2:371984/70‑1 (MQ=255)
                                                                                                                                                                                                                                                          gcgcTTACGAAGCCGCATCCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAACATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATTTAAACAATTTCAGAATAGATAAAAACTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCAtt                                                  >  1:11052/1‑240 (MQ=255)
                                                                                                                                                                                                                                                                aCGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAAc                                                                                                                                                >  7:879317/1‑140 (MQ=255)
                                                                                                                                                                                                                                                                aCGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAAc                                                                                                                                                <  8:879317/140‑1 (MQ=255)
                                                                                                                                                                                                                                                                     gCCGCA‑TCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTgat                                                                                                                                       >  6:765557/1‑144 (MQ=255)
                                                                                                                                                                                                                                                                     gCCGCA‑TCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTgat                                                                                                                                       <  5:765557/144‑1 (MQ=255)
                                                                                                                                                                                                                                                                        gCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAGACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATTTAAACAATTTCAGAATAGACAAAAACTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCATTATGAATATCTTACATATATGTTTGACCTCAAAATGGTTCAATATTGac  >  2:156635/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                          catCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTa                                                                                                                                                                             >  6:204973/3‑102 (MQ=39)
                                                                                                                                                                                                                                                                          catCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTa                                                                                                                                                                             <  5:204973/100‑1 (MQ=39)
                                                                                                                                                                                                                                                                                aCTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATTTaaa                                                                                                                       <  8:245155/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                                aCTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATTTaaa                                                                                                                       <  8:373290/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                                aCTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATTTaaa                                                                                                                       >  8:133194/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                        
AACAGAGCATGGGCTTTGACACCAACATGGAAGGCGGCTTCCTGGCGCGTCGTCGCCACCTACAGGCGCTGGAACAGGCAGCGGAACATCTACAACAGGGCAAAGCGCAACTGTTGGGAGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTAAAATTCTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATTTAAACAATTTCAGAATAGACAAAAACTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCATTATGAATATCTTACATATATGTGTGACCTCAAAATGGTTCAATATTGAC  >  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: ‑