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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,888,2120T. 8.2% 490.2 / 34.8 159intergenic (+20/‑223)mnmE/tnaCtRNA U34 5‑methylaminomethyl‑2‑thiouridine modification GTPase/tryptophanase leader peptide
Reads supporting (aligned to +/- strand):  ref base T (73/73);  new base . (9/4);  total (82/77)
Fisher's exact test for biased strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.23e-01

Only 100 of 161 total aligned reads displayed.

AGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACA  >  NC_000913/3888063‑3888359
                                                                                                                                                     |                                                                                                                                                      
 gCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCAtt                                                                                                                                                       >  1:1129336/1‑150 (MQ=255)
    tGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTg                                                                                                                                                    <  2:1581664/150‑1 (MQ=255)
      ggCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATATACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGAc                                                                                                                                                  <  2:397317/150‑1 (MQ=255)
            tGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAg                                                                                                                                            <  2:1638388/150‑1 (MQ=255)
            ggAACTGCTGGCGGAAGAGTGGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAg                                                                                                                                            <  1:168347/149‑1 (MQ=255)
             gAACTGCTGGCGGAAGAGTTGCGTCTGGCAAAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑a                                                                                                                                          >  1:885209/1‑150 (MQ=255)
               aCTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATg                                                                                                                                        >  2:1074880/1‑150 (MQ=255)
                 tgccgGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCg                                                                                                                                      <  2:1695781/150‑1 (MQ=255)
                    tgGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCA‑TCTGACAGTCAGAATGCGgct                                                                                                                                   <  1:1420469/150‑1 (MQ=255)
                     ggCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgctt                                                                                                                                  >  2:298613/1‑150 (MQ=255)
                       cGGAAGTGTTGCGTCTGGCACAGCAGAACTTAAGGAAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCTGATTTTCTCGAGCTTCTGTATTGTTAAGTAACCGCTCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgcttcg                                                                                                                                <  2:50455/150‑1 (MQ=255)
                       cGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgcttcg                                                                                                                                >  2:1186572/1‑150 (MQ=255)
                        ggAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgcttcgc                                                                                                                               >  1:922162/1‑150 (MQ=255)
                           agagTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTTACTGTCAG‑ATGCGgcttcgcttc                                                                                                                            <  2:1264052/150‑1 (MQ=255)
                               ttGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTg                                                                                                                        <  2:510227/150‑1 (MQ=255)
                                 gCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGtt                                                                                                                      <  1:337892/150‑1 (MQ=255)
                                  cggcTGGCACAGCAGATCTTAAGCGAAATCAGCGGGTATTTTACTTCAGAGGACCTGCTGGGGCTGATTTTCTCCAGCTTCTGTTTTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTa                                                                                                                     <  1:301330/147‑1 (MQ=255)
                                   gTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTAc                                                                                                                    >  1:1397491/1‑150 (MQ=255)
                                     cTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCa                                                                                                                  >  1:1449394/1‑150 (MQ=255)
                                       ggCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACt                                                                                                                <  1:274710/150‑1 (MQ=255)
                                         cacaGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTcc                                                                                                              >  1:1016417/1‑150 (MQ=255)
                                          acaGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCt                                                                                                             <  2:1200179/150‑1 (MQ=255)
                                           cagcagAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTg                                                                                                            <  2:777901/150‑1 (MQ=255)
                                             gcagAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGtt                                                                                                          <  1:384128/150‑1 (MQ=255)
                                                 aaCTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTc                                                                                                      <  1:835017/150‑1 (MQ=255)
                                                  aCTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTcc                                                                                                     <  2:60938/150‑1 (MQ=255)
                                                    ttAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACCAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTc                                                                                                   <  1:1061826/150‑1 (MQ=255)
                                                      aaGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCaa                                                                                                 <  1:60014/150‑1 (MQ=255)
                                                      aaGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCaa                                                                                                 >  1:635835/1‑150 (MQ=255)
                                                       aGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCCCTCCTGTTATTCCTCAAc                                                                                                >  1:1551944/1‑150 (MQ=255)
                                                       aGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAc                                                                                                >  2:1687797/1‑150 (MQ=255)
                                                       aGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAc                                                                                                <  2:890487/150‑1 (MQ=255)
                                                        gCGAAATCACCGGGTAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACCAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAcc                                                                                               <  1:1287612/150‑1 (MQ=255)
                                                         cGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAccc                                                                                              <  1:729488/150‑1 (MQ=255)
                                                          gAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCt                                                                                             >  1:1002600/1‑150 (MQ=255)
                                                            aaTCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCATCAACCCttt                                                                                           >  1:57775/1‑150 (MQ=255)
                                                               cACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCtttttt                                                                                        >  1:1713754/1‑150 (MQ=255)
                                                                 ccGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTa                                                                                      <  1:1327282/150‑1 (MQ=255)
                                                                   ggggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTaaa                                                                                    >  1:1227833/1‑150 (MQ=255)
                                                                    gggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAAc                                                                                   <  2:582520/150‑1 (MQ=255)
                                                                     ggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACa                                                                                  <  1:669679/150‑1 (MQ=255)
                                                                       attttACTTTAGACTATCTGCTGTGTCTGATTTTCTGCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCAGTCCTGTTATTCCTGAACCCTTTTTTTAAACAtt                                                                                <  2:295128/148‑1 (MQ=255)
                                                                        aTTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTa                                                                               >  1:787251/1‑150 (MQ=255)
                                                                           tACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTa‑aaa                                                                           >  1:890928/1‑150 (MQ=255)
                                                                            aCTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCTGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAAt                                                                          <  1:38808/150‑1 (MQ=255)
                                                                            aCTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAAt                                                                          <  1:478036/150‑1 (MQ=255)
                                                                             cTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAAGATTa‑taaat                                                                         >  2:688076/1‑148 (MQ=255)
                                                                             cTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAAtt                                                                         >  1:1226822/1‑150 (MQ=255)
                                                                             cTTCAGACGACCTGCTGGGGCGGATTTACTCCAGCTTCTGTATTTGTAAGTAACCGCGCTTACGAATCCGCATTCTGACTGACAG‑ATCCTGCTTCGCTTCATTTTTAACACTCCTGTTATTCCTCAACCCTTTTTTAAAACATTA‑AAAtt                                                                         >  2:320667/1‑150 (MQ=255)
                                                                                cAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑Ctt                                                                     >  1:1561615/1‑150 (MQ=255)
                                                                                 aGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCTCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTa                                                                    <  2:413664/150‑1 (MQ=255)
                                                                                  gacgacCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTAc                                                                   >  2:1642264/1‑150 (MQ=255)
                                                                                  gacgacCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTAc                                                                   <  1:1657709/150‑1 (MQ=255)
                                                                                   acgacCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACg                                                                  >  1:1546052/1‑150 (MQ=255)
                                                                                          gctgGGGCTTATTTTCTCCAGCTTCTGTATTGGTAAGAAACCGCGGTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTta                                                           <  2:1310433/150‑1 (MQ=255)
                                                                                          gctgGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTta                                                           <  2:1536698/150‑1 (MQ=255)
                                                                                             ggggCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATaa                                                        >  1:1568387/1‑150 (MQ=255)
                                                                                               ggCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATc                                                      <  2:1122591/150‑1 (MQ=255)
                                                                                                 gggATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCtt                                                    <  2:1441219/149‑1 (MQ=255)
                                                                                                 cGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCtt                                                    >  1:219468/1‑150 (MQ=255)
                                                                                                  ggATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCttt                                                   >  1:261019/1‑150 (MQ=255)
                                                                                                     ttttCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTaaa                                                <  2:1666189/150‑1 (MQ=255)
                                                                                                     ttttCTACAGCTTCTGTATTTGTAATTAACCGCTATTACGAAGCCTCATTCTGACTGTAAT‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTaaa                                                <  1:561351/150‑1 (MQ=255)
                                                                                                        tgtcCAGCTTCTGTATTGGTAAGTAACCGCTCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTaaaaaa                                             <  2:870061/148‑1 (MQ=255)
                                                                                                         ctcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTaaaaaaa                                            <  1:568926/150‑1 (MQ=255)
                                                                                                           ccAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑CTGCGGCTTCGCTTCATTGTTACCACTACTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGc                                          >  1:1297914/1‑150 (MQ=255)
                                                                                                           ccAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGc                                          >  1:1509454/1‑150 (MQ=255)
                                                                                                            cAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCa                                         <  2:1016417/150‑1 (MQ=255)
                                                                                                               cTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAAGATTATAAATTACGTAAG‑AATTT‑TAATGTTTAAAAAAAgggttg                                      <  1:1076422/150‑6 (MQ=21)
                                                                                                               cTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCAttt                                      >  1:1304496/1‑150 (MQ=255)
                                                                                                                 tCTGTTTTGGTATGTAACCGCGTTTACGAAGCCGCATTCTGACTGTCAT‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTACAATCTTTAAAAAAAGCATTTaa                                    <  2:1546052/150‑1 (MQ=255)
                                                                                                                   tGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAata                                  >  1:435972/1‑150 (MQ=255)
                                                                                                                     tATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATAtt                                >  2:1325185/1‑150 (MQ=255)
                                                                                                                         ggTGAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTc                            >  1:474216/1‑150 (MQ=255)
                                                                                                                         ggTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTc                            >  2:557898/1‑150 (MQ=255)
                                                                                                                         ggTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTc                            >  1:714347/1‑150 (MQ=255)
                                                                                                                         ggTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTc                            >  1:1513482/1‑150 (MQ=255)
                                                                                                                          gtaagtaaCCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTcc                           <  2:1369850/150‑1 (MQ=255)
                                                                                                                          gtaaataaCCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTcc                           >  1:366628/1‑150 (MQ=255)
                                                                                                                           taagtacCCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTccc                          >  1:670648/1‑150 (MQ=255)
                                                                                                                            aagtaaCCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTcccc                         <  2:1106456/150‑1 (MQ=255)
                                                                                                                             agtaaCCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATCCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCg                        <  1:1527056/150‑1 (MQ=255)
                                                                                                                                 aCCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACg                    <  1:1631230/150‑1 (MQ=255)
                                                                                                                                   cgcgCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGAt                  <  2:1227833/150‑1 (MQ=255)
                                                                                                                                    gcgcTTACTAAGCCTCATTCTTAATGTCAT‑TTGGGGTTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGAtt                 <  2:1551944/150‑1 (MQ=255)
                                                                                                                                      gcTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATtgt               <  2:1016925/150‑1 (MQ=255)
                                                                                                                                       cTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATtgtg              >  2:1563104/1‑150 (MQ=255)
                                                                                                                                       cTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATtgtg              <  1:1186572/150‑1 (MQ=255)
                                                                                                                                         tACGAAGCCGCATTCTGACTGTCAG‑CTGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCACCCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTgat            >  1:1626094/1‑150 (MQ=255)
                                                                                                                                         agcGAAGCCGCA‑TCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTgat            >  2:1160245/3‑150 (MQ=39)
                                                                                                                                          aCGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTTTAATCTTTAAAAAAAGCATTTAATATTGATCCCCGAACGATTGTgatt           >  1:170166/1‑150 (MQ=255)
                                                                                                                                           cGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTgattc          <  2:1509454/150‑1 (MQ=255)
                                                                                                                                            gAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGATCCCCGAACGATTGTgattcg         >  1:1360021/1‑150 (MQ=255)
                                                                                                                                              aGCCGCA‑TCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTgattcgat       >  1:1570208/1‑150 (MQ=255)
                                                                                                                                              aGCCGCA‑TCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTgattcgat       >  1:1043085/1‑150 (MQ=255)
                                                                                                                                                ccGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGTATTTAATATTGCTCCCCGAACGATTGTgattcgattc     >  2:1238347/1‑150 (MQ=255)
                                                                                                                                                ccGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTgattcgattc     <  2:922162/150‑1 (MQ=255)
                                                                                                                                                ccGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑ATTACGTAATTTATAATCTTTAAAAAAAGTATTTAATATTGCTACCCTAACGATTGTgattcgattc     >  2:805098/1‑150 (MQ=255)
                                                                                                                                                 cGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTca    >  1:693772/1‑150 (MQ=255)
                                                                                                                                                   cATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTcaca  >  1:1712196/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
AGCCTGGGCAGGTGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACA  >  NC_000913/3888063‑3888359

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑