Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 111,617 T→G 5.4% intergenic (+184/+32) 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,6170TG 5.4% 240.9 / 3.7 74intergenic (+184/+32)mutT/yacG8‑oxo‑dGTP diphosphatase/DNA gyrase inhibitor YacG
Reads supporting (aligned to +/- strand):  ref base T (30/40);  new base G (2/2);  total (32/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTTTCTTCAGCAGCCCATTCTCCGAGGTCGATCAGCTGACAACGTTTGGAGCAA  >  NC_000913/111468‑111766
                                                                                                                                                     |                                                                                                                                                     
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                                                                                                                                                       >  1:770116/1‑133 (MQ=255)
                 cGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCt                                                                                                                                                       <  2:770116/133‑1 (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                                                                                                                      >  1:1430620/1‑150 (MQ=255)
                                 tATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTc                                                                                                                      <  2:1342354/150‑1 (MQ=255)
                                   tCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCAc                                                                                                                    <  1:1000172/150‑1 (MQ=255)
                                       gCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCt                                                                                                                <  1:584754/150‑1 (MQ=255)
                                        ccTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCtt                                                                                                               >  2:1124919/1‑150 (MQ=255)
                                        ccTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCtt                                                                                                               <  2:1829294/150‑1 (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                                                                                                        >  1:355792/1‑113 (MQ=37)
                                                                                    cGGATGCGGCGTGAACGCCTTATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGttct                                                                                                        <  2:355792/113‑1 (MQ=37)
                                                                                     ggATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCActgct                                                                  <  1:635135/150‑1 (MQ=255)
                                                                                                 cACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGttt                                                      <  2:172455/150‑1 (MQ=255)
                                                                                                   cgccgcATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGttttt                                                    >  2:715519/1‑150 (MQ=255)
                                                                                                     ccgcATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTttct                                                  <  2:110020/150‑1 (MQ=255)
                                                                                                     ccgcATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTttct                                                  <  1:560733/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                                 <  1:1500096/150‑1 (MQ=255)
                                                                                                                      acaAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTTTCTTCAGCAGCCCATTCTcc                                 >  2:1800562/1‑150 (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:346783/7‑150 (MQ=255)
                                                                                                                                tgataagACGCGGACAGCGTCGCATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTTTCTTCAGCAGCCCATTCTCCGAGGTCGATc                       >  2:353803/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)
                                                                                                                                              ccGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTTTCTTCAGCAGCCCATTCTCCGAGGTCGATCAGCTGACAACGttt         >  1:1524634/1‑150 (MQ=255)
                                                                                                                                               cgcgTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTTTCTTCAGCAGCCCATTCTCCGAGGTCGATCAGCTGACAACGTTTg        <  1:1133573/150‑1 (MQ=255)
                                                                                                                                                 cgTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTTTCTTCAGCAGCCCATTCTCCGAGGTCGATCAGCTGACAACGTTTGGa      <  1:1124919/150‑1 (MQ=255)
                                                                                                                                                     ttATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTTTCTTCAGCAGCCCATTCTCCGAGGTCGATCAGCTGACAACGTTTGGAGCaa  >  2:240349/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTTTCTTCAGCAGCCCATTCTCCGAGGTCGATCAGCTGACAACGTTTGGAGCAA  >  NC_000913/111468‑111766

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