Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,888,219 T→A 7.7% 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,2190TA 7.7% 474.8 / 20.0 164intergenic (+27/‑216)mnmE/tnaCtRNA U34 5‑methylaminomethyl‑2‑thiouridine modification GTPase/tryptophanase leader peptide
Reads supporting (aligned to +/- strand):  ref base T (78/73);  new base A (9/4);  total (87/77)
Fisher's exact test for biased strand distribution p-value = 2.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 165 total aligned reads displayed.

TGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATTTAAACAA  >  NC_000913/3888075‑3888368
                                                                                                                                                |                                                                                                                                                        
tGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAg                                                                                                                                                     <  2:1638388/150‑1 (MQ=255)
 gAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGCGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑a                                                                                                                                                   >  1:1518408/1‑150 (MQ=255)
 gAACTGCTGGCGGAAGAGTTGCGTCTGGCAAAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑a                                                                                                                                                   >  1:885209/1‑150 (MQ=255)
     tgccgGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCg                                                                                                                                               <  2:1695781/150‑1 (MQ=255)
        tgGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCA‑TCTGACAGTCAGAATGCGgct                                                                                                                                            <  1:1420469/150‑1 (MQ=255)
          gcgcAAGAGTTGCGTCTGGCACAGCAGAACTTACGCGAAAGCACCGGGTAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgcttc                                                                                                                                          <  2:633448/150‑1 (MQ=255)
          gCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgcttc                                                                                                                                          <  1:391283/150‑1 (MQ=255)
           cGGAAGTGTTGCGTCTGGCACAGCAGAACTTAAGGAAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCTGATTTTCTCGAGCTTCTGTATTGTTAAGTAACCGCTCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgcttcg                                                                                                                                         <  2:50455/150‑1 (MQ=255)
           cGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgcttcg                                                                                                                                         >  2:1186572/1‑150 (MQ=255)
            ggAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGAATGTCAG‑ATGCGgcttcgc                                                                                                                                        >  1:84657/1‑150 (MQ=255)
             gAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgcttcgct                                                                                                                                       <  1:45035/150‑1 (MQ=255)
              aaGAGTTGCGTCTGGCACAGGAGAACTTAAGCGTAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgcttcgctt                                                                                                                                      <  1:342208/150‑1 (MQ=255)
               agagTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGgcttcgcttc                                                                                                                                     <  2:1455747/150‑1 (MQ=255)
                   ttGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTg                                                                                                                                 <  2:510227/150‑1 (MQ=255)
                      cggcTGGCACAGCAGATCTTAAGCGAAATCAGCGGGTATTTTACTTCAGAGGACCTGCTGGGGCTGATTTTCTCCAGCTTCTGTTTTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTa                                                                                                                              <  1:301330/147‑1 (MQ=255)
                       gTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTAc                                                                                                                             >  1:1397491/1‑150 (MQ=255)
                        tCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCTGCTTCGCTTCATTGTTAcc                                                                                                                            >  2:1044890/1‑150 (MQ=255)
                         cTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCa                                                                                                                           >  1:1449394/1‑150 (MQ=255)
                           ggCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACt                                                                                                                         <  1:274710/150‑1 (MQ=255)
                             cacaGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTcc                                                                                                                       >  2:1075826/1‑150 (MQ=255)
                             cacaGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTcc                                                                                                                       >  1:1016417/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)
                                 gcagAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGtt                                                                                                                   >  1:583501/1‑150 (MQ=255)
                                        ttAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACCAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTc                                                                                                            <  1:1061826/150‑1 (MQ=255)
                                          aaGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCaa                                                                                                          <  1:1210064/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)
                                             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)
                                                            aTTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTa                                                                                        >  1:1559717/1‑150 (MQ=255)
                                                               tACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTa‑aaa                                                                                    >  1:890928/1‑150 (MQ=255)
                                                                 cTTCAGACGACCTGCTGGGGCGGATTTACTCCAGCTTCTGTATTTGTAAGTAACCGCGCTTACGAATCCGCATTCTGACTGACAG‑ATCCTGCTTCGCTTCATTTTTAACACTCCTGTTATTCCTCAACCCTTTTTTAAAACATTA‑AAAtt                                                                                  >  2:320667/1‑150 (MQ=255)
                                                                   tCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑Ct                                                                               <  1:968728/150‑1 (MQ=255)
                                                                    cAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑Ctt                                                                              >  1:1561615/1‑150 (MQ=255)
                                                                     aGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTTTTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAATGTTTa‑aaaaaaggg                                                                              <  1:1645945/150‑7 (MQ=255)
                                                                     aGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCTCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTa                                                                             <  2:413664/150‑1 (MQ=255)
                                                                      gacgacCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTAc                                                                            <  1:1657709/150‑1 (MQ=255)
                                                                       acgacCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACg                                                                           >  1:1546052/1‑150 (MQ=255)
                                                                         gacCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGAATTCTTACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTa                                                                         <  1:1289271/150‑1 (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)
                                                                                     gggATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCtt                                                             <  2:1441219/149‑1 (MQ=255)
                                                                                     cGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCtt                                                             >  1:219468/1‑150 (MQ=255)
                                                                                         ttttCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTaaa                                                         <  2:1666189/150‑1 (MQ=255)
                                                                                            tgtcCAGCTTCTGTATTGGTAAGTAACCGCTCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTaaaaaa                                                      <  2:870061/148‑1 (MQ=255)
                                                                                            tctcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTaaaaaa                                                      <  2:1242299/150‑1 (MQ=255)
                                                                                             ctcCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTaaaaaaa                                                     <  1:568926/150‑1 (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                                               >  2:1656658/1‑150 (MQ=255)
                                                                                                   cTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCAttt                                               >  1:1304496/1‑150 (MQ=255)
                                                                                                     tCTGTTTTGGTATGTAACCGCGTTTACGAAGCCGCATTCTGACTGTCAT‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTACAATCTTTAAAAAAAGCATTTaa                                             <  2:1546052/150‑1 (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                                     >  1:219886/1‑150 (MQ=255)
                                                                                                             ggTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTc                                     >  2:557898/1‑150 (MQ=255)
                                                                                                             ggTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTc                                     >  1:714347/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)
                                                                                                                       cgcgCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGAt                           <  2:1227833/150‑1 (MQ=255)
                                                                                                                        gcgcTTACTAAGCCTCATTCTTAATGTCAT‑TTGGGGTTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGAtt                          <  2:1551944/150‑1 (MQ=255)
                                                                                                                           cTTACGTAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATtgtg                       <  2:1568387/150‑1 (MQ=255)
                                                                                                                           cTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATtgtg                       <  1:1186572/150‑1 (MQ=255)
                                                                                                                           cTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATtgtg                       >  2:1563104/1‑150 (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)
                                                                                                                               cGAAGCCGCA‑TCTGACAGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTgattcg                  >  2:136670/1‑150 (MQ=255)
                                                                                                                                gAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGATCCCCGAACGATTGTgattcg                  >  1:1360021/1‑150 (MQ=255)
                                                                                                                                gAAGCCGCA‑TCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTgattcg                  >  1:549910/1‑150 (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:1570208/1‑150 (MQ=255)
                                                                                                                                  aGCCGCA‑TCTGACAGTCAGAATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTgattcgat                >  1:1043085/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)
                                                                                                                                     cGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTca             >  1:693772/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       >  2:726703/1‑150 (MQ=255)
                                                                                                                                           cTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATTTa       <  1:519919/150‑1 (MQ=255)
                                                                                                                                            tGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCTATTCACATTTaa      >  1:1676615/1‑150 (MQ=255)
                                                                                                                                                tGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTACTGTTTTTCCTCAACCCTTTTTTTAAACATTA‑ATATT‑CTTACGTATTTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTAACATTTAAACaa  >  2:594542/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                        
TGAACTGCTGGCGGAAGAGTTGCGTCTGGCACAGCAGAACTTAAGCGAAATCACCGGGGAATTTACTTCAGACGACCTGCTGGGGCGGATTTTCTCCAGCTTCTGTATTGGTAAGTAACCGCGCTTACGAAGCCGCATTCTGACTGTCAG‑ATGCGGCTTCGCTTCATTGTTACCACTCCTGTTATTCCTCAACCCTTTTTTTAAACATTA‑AAATT‑CTTACGTAATTTATAATCTTTAAAAAAAGCATTTAATATTGCTCCCCGAACGATTGTGATTCGATTCACATTTAAACAA  >  NC_000913/3888075‑3888368

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