Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 111,609 T→A 5.0% intergenic (+176/+40) mutT → / ← yacG 8‑oxo‑dGTP diphosphatase/DNA gyrase inhibitor YacG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913111,6090TA 5.0% 253.2 / 2.4 79intergenic (+176/+40)mutT/yacG8‑oxo‑dGTP diphosphatase/DNA gyrase inhibitor YacG
Reads supporting (aligned to +/- strand):  ref base T (32/43);  new base A (2/2);  total (34/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

CGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTTTCTTCAGCAGCCCATTCTCCGAGGTCGATCAGCTGACAA  >  NC_000913/111460‑111754
                                                                                                                                                     |                                                                                                                                                 
cGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGt                                                                                                                                                   >  1:1342354/1‑150 (MQ=255)
 gACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTc                                                                                                                                                  <  1:1837940/150‑1 (MQ=255)
 gACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTc                                                                                                                                                  <  1:183421/150‑1 (MQ=255)
     ttCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCcgcg                                                                                                                                              <  1:1615072/150‑1 (MQ=255)
        gCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCt                                                                                                                                           <  1:1007759/150‑1 (MQ=255)
        gCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCt                                                                                                                                           >  2:717443/1‑150 (MQ=255)
         cacaCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCCAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCtt                                                                                                                                          >  1:123831/1‑150 (MQ=255)
           cacGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGTGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTAt                                                                                                                                        >  2:1828264/1‑150 (MQ=255)
           cacGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTAt                                                                                                                                        >  2:551390/1‑150 (MQ=255)
                  ccTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCt                                                                                                                                 >  2:1244799/1‑150 (MQ=255)
                    tGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTAc                                                                                                                               <  2:474138/150‑1 (MQ=255)
                    tGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATAAGACGCGGACAGCGTCGCATCAGGcatct                                                                                                                               <  1:742813/150‑5 (MQ=38)
                      aTGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGt                                                                                                                             >  1:1669227/1‑150 (MQ=255)
                         cGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCt                                                                                                                                           <  2:770116/133‑1 (MQ=255)
                         cGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCt                                                                                                                                           >  1:770116/1‑133 (MQ=255)
                         cGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCg                                                                                                                          >  1:846503/1‑150 (MQ=255)
                          gACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCgg                                                                                                                         >  1:1998394/1‑150 (MQ=255)
                            cgcgccagcgTCTTATAAGCCATAAAGCGATTTCTAACTAATAGATAAATTTGTTTTAGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCa                                                                                                                       >  2:1887444/8‑150 (MQ=25)
                               tGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCa                                                                                                                    <  1:1909366/150‑1 (MQ=255)
                                ggCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAg                                                                                                                   >  1:1720611/1‑150 (MQ=255)
                                     gTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACaaa                                                                                                              <  1:102167/150‑1 (MQ=255)
                                        ttATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGt                                                                                                           >  2:1408498/1‑150 (MQ=255)
                                         tATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTc                                                                                                          <  2:1342354/150‑1 (MQ=255)
                                         tATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTc                                                                                                          >  1:1430620/1‑150 (MQ=255)
                                           tCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCAc                                                                                                        <  1:1000172/150‑1 (MQ=255)
                                               gCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCt                                                                                                    <  1:584754/150‑1 (MQ=255)
                                                ccTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCtt                                                                                                   <  2:1829294/150‑1 (MQ=255)
                                                ccTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCtt                                                                                                   >  2:1124919/1‑150 (MQ=255)
                                                 cTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCttt                                                                                                  <  2:443416/150‑1 (MQ=255)
                                                  tacaGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTg                                                                                                 >  1:329978/4‑150 (MQ=255)
                                                  tAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTg                                                                                                 <  2:591059/150‑1 (MQ=255)
                                                   aaaGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTgg                                                                                                <  1:799401/150‑1 (MQ=255)
                                                    aaGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGt                                                                                               <  1:536682/150‑1 (MQ=255)
                                                    aaGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGt                                                                                               <  1:551390/150‑1 (MQ=255)
                                                     aGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGtt                                                                                              <  1:1392588/150‑1 (MQ=255)
                                                         aTTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGttcttc                                                                                          <  1:1408498/150‑1 (MQ=255)
                                                            tCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCt                                                                                       <  1:1552900/150‑1 (MQ=255)
                                                              tAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTcc                                                                                     <  2:1720611/150‑1 (MQ=255)
                                                                 cTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGGTcattt                                                                                  >  2:1294902/1‑146 (MQ=255)
                                                                 cTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGt                                                                                  <  2:516528/150‑1 (MQ=255)
                                                                 cTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTGCTTCGCTCCAGt                                                                                  <  1:1634419/150‑1 (MQ=255)
                                                                      tGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCg                                                                             >  1:1505485/1‑150 (MQ=255)
                                                                       gATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGc                                                                            >  2:350957/1‑150 (MQ=255)
                                                                          ataTTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCTCAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCttt                                                                         <  2:1635337/148‑1 (MQ=255)
                                                                           aaTTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTc                                                                        <  2:329978/150‑1 (MQ=255)
                                                                           aaTTTGTTTTTGTAGGTCGGATAAGGCGTTAACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCCGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTc                                                                        >  2:1340545/1‑150 (MQ=255)
                                                                             tttGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCgg                                                                      <  1:648018/150‑1 (MQ=255)
                                                                               tGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGaa                                                                    <  1:1460271/150‑1 (MQ=255)
                                                                                 tttttGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAg                                                                  >  1:110020/1‑150 (MQ=255)
                                                                                  ttttGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGa                                                                 <  1:620239/150‑1 (MQ=255)
                                                                                       tAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGcc                                                            >  1:793907/1‑150 (MQ=255)
                                                                                        aGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCa                                                           <  1:717443/150‑1 (MQ=255)
                                                                                            cGGATGCGGCGTGAACGCCTTATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGttct                                                                                            <  2:355792/113‑1 (MQ=37)
                                                                                            cGGATGCGGCGTGAACGCCTTATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGttct                                                                                            >  1:355792/1‑113 (MQ=37)
                                                                                             ggATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCActgct                                                      <  1:635135/150‑1 (MQ=255)
                                                                                                         cACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGttt                                          <  2:172455/150‑1 (MQ=255)
                                                                                                           cgccgcATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGttttt                                        >  2:715519/1‑150 (MQ=255)
                                                                                                             ccgcATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTttct                                      <  1:560733/150‑1 (MQ=255)
                                                                                                             ccgcATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTttct                                      <  2:110020/150‑1 (MQ=255)
                                                                                                             ccgcATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGCAAAGATCGCCACTGCTGGGGATTCGTTTttct                                      >  1:1302593/1‑150 (MQ=255)
                                                                                                                  tCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCTTCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTTTCTTcagc                                 <  1:1828264/150‑1 (MQ=255)
                                                                                                                  tCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTTTCTTcagc                                 <  2:793907/150‑1 (MQ=255)
                                                                                                                    cGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTTTCTTcagcag                               >  1:299038/1‑150 (MQ=255)
                                                                                                                    cGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTTTCTTcagcag                               >  1:1861382/1‑150 (MQ=255)
                                                                                                                      aCATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTTTCTTCAGCAGcc                             <  2:846503/150‑1 (MQ=255)
                                                                                                                         tttGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTTTCTTCAGCAGCCCAt                          >  2:596793/1‑150 (MQ=255)
                                                                                                                              acaAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTTTCTTCAGCAGCCCATTCTcc                     >  2:1800562/1‑150 (MQ=255)
                                                                                                                              acaAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTTTCTTCAGCAGCCCATTCTcc                     <  1:1500096/150‑1 (MQ=255)
                                                                                                                               caAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTTTCTTCAGCAGCCCATTCTCCg                    >  2:533283/1‑150 (MQ=255)
                                                                                                                                    tGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTTTCTTCAGCAGCCCATTCTCCGAGGTc               >  2:229672/1‑150 (MQ=255)
                                                                                                                                    tGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTTTCTTCAGCAGCCCATTCTCCGAGGTc               <  2:1998394/150‑1 (MQ=255)
                                                                                                                                        tgataagACGCGGACAGCGTCGCATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTTTCTTCAGCAGCCCATTCTCCGAGGTCGATc           >  2:353803/7‑150 (MQ=255)
                                                                                                                                        tgataagACGCGGACAGCGTCGCATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTTTCTTCAGCAGCCCATTCTCCGAGGTCGATc           >  2:346783/7‑150 (MQ=255)
                                                                                                                                        tGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTTTCTTCAGCAGCCCATTCTCCGAGGTCGATc           <  2:1430620/150‑1 (MQ=255)
                                                                                                                                         gATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTTTCTTCAGCAGCCCATTCTCCGAGGTCGATCa          >  2:88122/1‑150 (MQ=255)
                                                                                                                                          ataagACGCGGACAGCGTCGCATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTTTCTTCAGCAGCCCATTCTCCGAGGTCGATCAg         <  1:1424978/146‑1 (MQ=255)
                                                                                                                                          aTGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTTTCTTCAGCAGCCCATTCTCCGAGGTCGATCAg         <  1:1942124/150‑1 (MQ=255)
                                                                                                                                             cGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTTTCTTCAGCAGCCCATTCTCCGAGGTCGATCAGCTg      <  1:715519/150‑1 (MQ=255)
                                                                                                                                                cGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTTTCTTCAGCAGCCCATTCTCCGAGGTCGATCAGCTGACa   <  1:229672/150‑1 (MQ=255)
                                                                                                                                                 gCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTTTCTTCAGCAGCCCATTCTCCGAGGTCGATCAGCTGACaa  >  1:292481/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTTTCTTCAGCAGCCCATTCTCCGAGGTCGATCAGCTGACAA  >  NC_000913/111460‑111754

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