Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,888,224:1 +A 6.8% 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 6.8% 603.2 / 41.4 161intergenic (+32/‑211)mnmE/tnaCtRNA U34 5‑methylaminomethyl‑2‑thiouridine modification GTPase/tryptophanase leader peptide
Reads supporting (aligned to +/- strand):  ref base . (79/71);  new base A (8/3);  total (87/74)
Fisher's exact test for biased strand distribution p-value = 2.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.20e-02

Only 100 of 168 total aligned reads displayed.

TGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATTTAAACAATTTCA  >  NC_000913/3888075‑3888373
                                                                                                                                                      |                                                                                                                                                       
tGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAg                                                                                                                                                          <  2:1638388/150‑1 (MQ=255)
ggAACTGCTGGCGGAAGAGTGGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAg                                                                                                                                                          <  1:168347/149‑1 (MQ=255)
 gAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑a                                                                                                                                                        >  2:46765/1‑150 (MQ=255)
 gAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGCGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑a                                                                                                                                                        >  1:1518408/1‑150 (MQ=255)
   aCTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATg                                                                                                                                                      >  2:1074880/1‑150 (MQ=255)
     ggctgGCGGAAGATTTGCGTCTGGCTCAGCAGAACTTATGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCA‑TCTGACAGTCAG‑aatgcg                                                                                                                                                   <  1:1563952/149‑6 (MQ=255)
      gctgGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCgg                                                                                                                                                   >  1:1622867/1‑150 (MQ=255)
       ctgGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGCCGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgc                                                                                                                                                  >  1:889216/1‑150 (MQ=255)
        tgGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCA‑TCTGACAGTCAGAATGCGgct                                                                                                                                                 <  1:1420469/150‑1 (MQ=255)
         ggCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgctt                                                                                                                                                >  2:298613/1‑150 (MQ=255)
         ggCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCA‑TCTGACAGTCAGAATGCGgctt                                                                                                                                                <  2:1679073/150‑1 (MQ=255)
          gcgcAAGAGTTGCGTCTGGCACAGCAGAACTTACGCGAAAGCACCGGGTAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgcttc                                                                                                                                               <  2:633448/150‑1 (MQ=255)
          gCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgcttc                                                                                                                                               <  1:391283/150‑1 (MQ=255)
           cGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgcttcg                                                                                                                                              >  2:1186572/1‑150 (MQ=255)
            ggAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgcttcgc                                                                                                                                             >  1:922162/1‑150 (MQ=255)
            ggAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGAATGTCAG‑ATGCGgcttcgc                                                                                                                                             >  1:84657/1‑150 (MQ=255)
             gAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgcttcgct                                                                                                                                            <  1:45035/150‑1 (MQ=255)
               agagTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgcttcgcttc                                                                                                                                          <  2:1455747/150‑1 (MQ=255)
                   ttGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTg                                                                                                                                      <  2:510227/150‑1 (MQ=255)
                    tGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGt                                                                                                                                     >  1:1369850/1‑150 (MQ=255)
                     gCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGtt                                                                                                                                    <  1:337892/150‑1 (MQ=255)
                      cGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTa                                                                                                                                   >  1:1666189/1‑150 (MQ=255)
                       gTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTAc                                                                                                                                  >  1:1397491/1‑150 (MQ=255)
                             cacaGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTcc                                                                                                                            >  1:1016417/1‑150 (MQ=255)
                             cacaGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTcc                                                                                                                            >  2:1075826/1‑150 (MQ=255)
                              acaGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCt                                                                                                                           <  2:1270410/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)
                                          aaGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCaa                                                                                                               <  1:742553/150‑1 (MQ=255)
                                          aaGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCaa                                                                                                               <  1:60014/150‑1 (MQ=255)
                                          aaGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCaa                                                                                                               >  1:635835/1‑150 (MQ=255)
                                           aGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAc                                                                                                              <  2:890487/150‑1 (MQ=255)
                                           aGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAc                                                                                                              >  2:1687797/1‑150 (MQ=255)
                                            gCGAAATCACCGGGTAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACCAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAAcc                                                                                                             <  1:1287612/150‑1 (MQ=255)
                                              gAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCt                                                                                                           >  1:1002600/1‑150 (MQ=255)
                                                aaTCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCATCAACCCttt                                                                                                         >  1:57775/1‑150 (MQ=255)
                                                     ccGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTa                                                                                                    <  1:1327282/150‑1 (MQ=255)
                                                       ggggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTaaa                                                                                                  >  1:1227833/1‑150 (MQ=255)
                                                         ggAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACa                                                                                                <  1:669679/150‑1 (MQ=255)
                                                            aTTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTa                                                                                             >  1:1559717/1‑150 (MQ=255)
                                                                aCTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAAt                                                                                        <  1:699657/150‑1 (MQ=255)
                                                                 cTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAAGATTa‑taaat                                                                                       >  2:688076/1‑148 (MQ=255)
                                                                 cTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAAtt                                                                                       >  2:568926/1‑150 (MQ=255)
                                                                 cTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAAtt                                                                                       >  1:1226822/1‑150 (MQ=255)
                                                                   tCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑Ct                                                                                    <  1:968728/150‑1 (MQ=255)
                                                                      gacgacCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTAc                                                                                 <  1:1657709/150‑1 (MQ=255)
                                                                       acgacCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACg                                                                                >  1:1546052/1‑150 (MQ=255)
                                                                             tgctgGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAAttt                                                                          >  2:1631230/1‑150 (MQ=255)
                                                                              gctgGGGCTTATTTTCTCCAGCTTCTGTATTGGTAAGAAACCGCGGTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTta                                                                         <  2:1310433/150‑1 (MQ=255)
                                                                              gctgGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTta                                                                         <  2:1536698/150‑1 (MQ=255)
                                                                              gctgGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTta                                                                         <  1:298613/150‑1 (MQ=255)
                                                                                 ggggCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATaa                                                                      >  1:1568387/1‑150 (MQ=255)
                                                                                     gtgATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCTCTTACGAAGCCGCTTTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTTATTTATAATCtt                                                                  <  1:1529098/148‑1 (MQ=255)
                                                                                     gggATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCtt                                                                  <  2:1441219/149‑1 (MQ=255)
                                                                                      ggATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCttt                                                                 >  1:261019/1‑150 (MQ=255)
                                                                                         ttttCTACAGCTTCTGTATTTGTAATTAACCGCTATTACGAAGCCTCATTCTGACTGTAAT‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTaaa                                                              <  1:561351/150‑1 (MQ=255)
                                                                                            tctcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTaaaaaa                                                           <  2:1242299/150‑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)
                                                                                                 aGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAAGATTATAAATTACGTAAG‑AATTT‑TAATGTTTAAAAAAAgggt                                                      <  1:688076/150‑4 (MQ=25)
                                                                                                   cTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAAGATTATAAATTACGTAAG‑AATTT‑TAATGTTTAAAAAAAgggttg                                                    <  1:1076422/150‑6 (MQ=21)
                                                                                                   cTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCAttt                                                    >  1:1304496/1‑150 (MQ=255)
                                                                                                   cTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCAttt                                                    >  2:1656658/1‑150 (MQ=255)
                                                                                                           ttGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCACCCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTga                                            >  1:370920/1‑149 (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:219886/1‑150 (MQ=255)
                                                                                                             ggTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTc                                          >  1:714347/1‑150 (MQ=255)
                                                                                                                aagtaaCCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTcccc                                       <  2:1106456/150‑1 (MQ=255)
                                                                                                                   taaCCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAATCATTTAATATTGCTCCCCGaa                                    >  2:1547410/1‑150 (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‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTgat                          >  1:1486061/1‑150 (MQ=255)
                                                                                                                             agcGAAGCCGCA‑TCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTgat                          >  2:1160245/3‑150 (MQ=39)
                                                                                                                               cGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTgattc                        <  1:1075826/150‑1 (MQ=255)
                                                                                                                               cGAAGCCGCA‑TCTGACAGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTgattcg                       >  2:136670/1‑150 (MQ=255)
                                                                                                                                gAAGCCGCA‑TCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGAATTTAATATTGATCCACGAACGATTGTgattct                       >  2:1707433/1‑149 (MQ=255)
                                                                                                                                 aaGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTgattcga                      >  2:1529175/1‑150 (MQ=255)
                                                                                                                                  aGCCGCA‑TCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCACCCCTTTTTTTAAACATTT‑AAAGT‑CTTACGGAATTTATAATCTTTTAAAAAAGCATTTAATATTGCTCCCCGAACGATTGGgattcgat                     >  2:318144/1‑150 (MQ=255)
                                                                                                                                  aGCCGCA‑TCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTgattcgat                     >  1:1043085/1‑150 (MQ=255)
                                                                                                                                   gCCGCA‑TCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATAA‑AAAAT‑ATAACGTAATTTATAATCAATAAAAAAAGCATTAAATATTGATCAACGAACGATTGTgattccatt                    >  1:1716333/1‑150 (MQ=255)
                                                                                                                                    ccGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTgattcgattc                   <  2:922162/150‑1 (MQ=255)
                                                                                                                                    ccGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTgattcgattc                   <  2:1226822/150‑1 (MQ=255)
                                                                                                                                    ccGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑ATTACGTAATTTATAATCTTTAAAAAAAGTATTTAATATTGCTACCCTAACGATTGTgattcgattc                   >  2:805098/1‑150 (MQ=255)
                                                                                                                                    cagCATTCTGACTGTCAT‑ATGCGGCTTCTCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTgattcgattc                   <  2:84657/148‑1 (MQ=255)
                                                                                                                                       cATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTcaca                >  1:1712196/1‑150 (MQ=255)
                                                                                                                                       cATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTcaca                >  1:1705665/1‑150 (MQ=255)
                                                                                                                                       cATTATGACTGTCAG‑ATGCTGCTTCTCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTcaca                <  2:1129336/150‑1 (MQ=255)
                                                                                                                                          tCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACAttt             >  1:308158/1‑150 (MQ=255)
                                                                                                                                           cTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATTTa            <  1:519919/150‑1 (MQ=255)
                                                                                                                                             gACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATTTaaa          >  2:37616/1‑150 (MQ=255)
                                                                                                                                                tGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTACTGTTTTTCCTCAACCCTTTTTTTAAACATTA‑ATATT‑CTTACGTATTTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTAACATTTAAACaa       >  2:594542/1‑150 (MQ=255)
                                                                                                                                                 gTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTAACATTTAAACAAt      >  2:481017/1‑150 (MQ=255)
                                                                                                                                                     g‑ATGCGGCTTCGCTTCATTGTTACCACTCCAGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATTTAAACAATTTCa  >  1:356610/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                       
TGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATTTAAACAATTTCA  >  NC_000913/3888075‑3888373

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: ‑