Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,017,222 Δ1 bp 9.2% intergenic (+30/‑111) udp → / → rmuC uridine phosphorylase/putative recombination limiting protein RmuC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,017,2220T. 9.2% 286.9 / 23.5 98intergenic (+30/‑111)udp/rmuCuridine phosphorylase/putative recombination limiting protein RmuC
Reads supporting (aligned to +/- strand):  ref base T (48/41);  new base . (2/7);  total (50/48)
Fisher's exact test for biased strand distribution p-value = 8.83e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

GGTAGCGGGTGTTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGT  >  NC_000913/4017075‑4017369
                                                                                                                                                    |                                                                                                                                                   
ggTAGCGGGTGTTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTc                                                                                                                                                   >  2:1239435/1‑150 (MQ=255)
ggTAGCGGGTGTTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTc                                                                                                                                                   <  2:1073157/150‑1 (MQ=255)
 gTAGCGGGTGTTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCg                                                                                                                                                  <  2:681402/150‑1 (MQ=255)
    gCGGGTGTTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGcc                                                                                                                                               <  2:608329/150‑1 (MQ=255)
       ggTGTTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttc                                                                                                                                            <  1:1098262/150‑2 (MQ=255)
       ggTGTTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCttt                                                                                                                                            >  2:1228954/1‑150 (MQ=255)
        gtgtTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttca                                                                                                                                           <  1:1498424/150‑3 (MQ=255)
        gtgtTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttca                                                                                                                                           <  2:1449183/150‑3 (MQ=255)
        gtgtTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttca                                                                                                                                           <  2:45703/150‑3 (MQ=255)
        gtgtTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttca                                                                                                                                           <  1:1036936/150‑3 (MQ=255)
         tgtTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttcag                                                                                                                                          <  1:917475/150‑4 (MQ=255)
           ttATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTa                                                                                                                                        <  1:336729/150‑1 (MQ=255)
           ttATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTAAAGGCCGAACGCGTCGGCCttcagac                                                                                                                                        <  2:1003562/150‑6 (MQ=255)
            tATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTAt                                                                                                                                       >  2:1531932/1‑150 (MQ=255)
            tATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTAt                                                                                                                                       <  1:656193/150‑1 (MQ=255)
                gTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGCAGCGGCGCGTCGTCTGGAGTAATTGTCTTCTCCTGTCAGAAGGCCGACGAGTTCGGAGTTTTTTAttttt                                                                                                                                   >  1:478522/1‑150 (MQ=255)
                gTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTAttttt                                                                                                                                   >  1:146301/1‑150 (MQ=255)
                 ttAACCGCACCCAGCAAGCGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGTGGCGCTTTGTCTGCTGTAATTCTGTTCTCTTGTCTTAAGGCCGACGTGTTCGTCCTTTTGTATTTTTg                                                                                                                                  <  2:214238/150‑1 (MQ=255)
                  tAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGc                                                                                                                                 >  2:314086/1‑150 (MQ=255)
                        cACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAgcg                                                                                                                           <  2:1156113/150‑1 (MQ=255)
                          cccAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGcc                                                                                                                         >  1:662845/1‑150 (MQ=255)
                            cAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTc                                                                                                                       >  2:1477858/1‑150 (MQ=255)
                              gCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGc                                                                                                                     >  1:736106/1‑150 (MQ=255)
                                aaGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAg                                                                                                                   <  2:1318524/150‑1 (MQ=255)
                                aaGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAg                                                                                                                   <  2:868975/150‑1 (MQ=255)
                                aaGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAg                                                                                                                   >  1:575444/1‑150 (MQ=255)
                                 agagATCCCGAATGCTGAGACGATGAAACAAACCTAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAgg                                                                                                                  <  1:821961/150‑1 (MQ=255)
                                  gagaTCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGa                                                                                                                 <  1:1183578/150‑1 (MQ=255)
                                  gagaTCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGa                                                                                                                 <  1:1343220/150‑1 (MQ=255)
                                  gagaTCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGa                                                                                                                 <  2:714839/150‑1 (MQ=255)
                                       cccGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGc                                                                                                            >  2:941650/1‑150 (MQ=255)
                                       cccGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGc                                                                                                            >  1:673085/1‑150 (MQ=255)
                                       cccGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGc                                                                                                            >  1:791654/1‑150 (MQ=255)
                                        ccGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGcc                                                                                                           <  1:1470342/150‑1 (MQ=255)
                                        cagAATGCTGAGACGATGAAACAAACCGAAACCCATGCGTTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGGTGTAATTCTCTTCTCCTGTCTGAATGCCGACTCGTTCGGCCTTTTGTATTTTTGCGTAGCGCATCGCAGGAAATGcc                                                                                                           <  2:499870/148‑1 (MQ=255)
                                         cGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCt                                                                                                          >  2:789254/1‑150 (MQ=255)
                                          gAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCtt                                                                                                         >  1:946987/1‑150 (MQ=255)
                                              aaggaggCGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCa                                                                                                     <  2:211511/147‑1 (MQ=255)
                                               cTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCaa                                                                                                    >  1:1075828/1‑150 (MQ=255)
                                                   gaCGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAAGGCATTTcctgcgag                                                                                                <  1:1106961/150‑7 (MQ=255)
                                                         gAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATCCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGttt                                                                                          <  2:330077/150‑1 (MQ=255)
                                                          aaacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTTGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTg                                                                                         >  2:1388658/1‑150 (MQ=255)
                                                          aaacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTg                                                                                         <  1:1228954/150‑1 (MQ=255)
                                                          aaacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTg                                                                                         <  1:1598094/150‑1 (MQ=255)
                                                          aaacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTg                                                                                         <  1:1531932/150‑1 (MQ=255)
                                                           aacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGTAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGt                                                                                        >  1:884893/1‑150 (MQ=255)
                                                           aacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGt                                                                                        >  1:124932/1‑150 (MQ=255)
                                                             caaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTAc                                                                                      <  1:1496313/150‑1 (MQ=255)
                                                             caaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTAc                                                                                      >  1:1631483/1‑150 (MQ=255)
                                                              aaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACa                                                                                     >  1:549692/1‑150 (MQ=255)
                                                                acCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGc                                                                                   <  2:147663/150‑1 (MQ=255)
                                                                   gAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGcaca                                                                                <  1:941650/150‑1 (MQ=255)
                                                                     aaGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAAt                                                                              >  2:249667/1‑150 (MQ=255)
                                                                      aGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAAtt                                                                             >  2:1355203/1‑150 (MQ=255)
                                                                          aTGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTAt                                                                         >  1:792955/1‑150 (MQ=255)
                                                                             cGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTAtttt                                                                      >  2:418722/1‑150 (MQ=255)
                                                                                  aaaaTCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCg                                                                 <  2:470989/150‑1 (MQ=255)
                                                                                  aaaaTCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCg                                                                 <  2:1028123/150‑1 (MQ=255)
                                                                                  aaaaTCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCg                                                                 <  1:225028/150‑1 (MQ=255)
                                                                                      tCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTa                                                             >  2:610756/1‑150 (MQ=255)
                                                                                       cGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTATAGCACAATTCTATTTTGTGCGGGTaa                                                            >  1:1360431/1‑150 (MQ=255)
                                                                                       cGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTaa                                                            >  2:667967/1‑150 (MQ=255)
                                                                                            tggAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAgttgt                                                       >  2:1421245/1‑150 (MQ=255)
                                                                                             ggAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAGTTCTATTTTGTGCGGGTAAgttgtt                                                      >  1:362090/1‑150 (MQ=255)
                                                                                              gAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAATCACGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGttgttg                                                     >  2:154955/1‑150 (MQ=255)
                                                                                                aGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCg                                                   >  2:1170323/1‑150 (MQ=255)
                                                                                                  cggcgCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTc                                                 <  2:662845/150‑1 (MQ=255)
                                                                                                   ggcgCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCa                                                <  1:163317/150‑1 (MQ=255)
                                                                                                   ggcgCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCa                                                >  1:875948/1‑150 (MQ=255)
                                                                                                     cgcgTCGGCTGCTGTATTTCTCTTCTCCTGTCTGAAGTCCGACGGGTTCGGACTTTTGTATTTTTGCGTAGGGCCTCGCAGGAAATGCCTTTCGAACTGGTCGTTTGTACAGCACTATTCTTTTTTGTGGGGGTAAGTTGTTGCGTCagg                                              <  2:673085/150‑1 (MQ=255)
                                                                                                      gcgTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTTGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCagga                                             <  1:1388658/150‑1 (MQ=255)
                                                                                                        gtcgtcTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCaggagg                                           >  2:1056157/1‑150 (MQ=255)
                                                                                                          cgtcTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCg                                         <  2:884893/150‑1 (MQ=255)
                                                                                                            tcTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGtt                                       >  2:1488461/1‑150 (MQ=255)
                                                                                                               gctgTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGTCGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTtgtg                                    <  1:1239435/150‑1 (MQ=255)
                                                                                                               gctgTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTtgtg                                    >  2:927786/1‑150 (MQ=255)
                                                                                                               gctgTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTtgtg                                    <  2:549692/150‑1 (MQ=255)
                                                                                                               gctgTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGGTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTtgtg                                    <  2:1544123/150‑1 (MQ=255)
                                                                                                                ctgTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTgg                                   >  1:478649/1‑150 (MQ=255)
                                                                                                                ctgTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTgg                                   >  1:1601578/1‑150 (MQ=255)
                                                                                                                  gTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGAt                                 >  1:733355/1‑150 (MQ=255)
                                                                                                                          cttctcCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAAt                         <  2:736106/150‑1 (MQ=255)
                                                                                                                           ttctcCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCACCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATc                        >  2:885123/1‑150 (MQ=255)
                                                                                                                             ctcCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCAt                      >  2:1129446/1‑150 (MQ=255)
                                                                                                                             ctcCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCAt                      <  2:820005/150‑1 (MQ=255)
                                                                                                                              tcCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATg                     <  2:977801/150‑1 (MQ=255)
                                                                                                                               ccTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATgg                    >  2:1444011/1‑150 (MQ=255)
                                                                                                                               ccTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAAATCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATgg                    >  1:969016/1‑150 (MQ=255)
                                                                                                                                 tGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGtt                  >  2:775659/1‑150 (MQ=255)
                                                                                                                                  gTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGttt                 >  1:570576/1‑150 (MQ=255)
                                                                                                                                  gTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGttt                 <  2:146301/150‑1 (MQ=255)
                                                                                                                                    caaaaGGCCGAACGCGTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTAc               >  1:281806/4‑150 (MQ=255)
                                                                                                                                    caaaaGGCCGAACGCGTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTAc               >  1:1518573/4‑150 (MQ=255)
                                                                                                                                     tGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACg              >  1:1407861/1‑150 (MQ=255)
                                                                                                                                     tGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACg              >  2:467236/1‑150 (MQ=255)
                                                                                                                                      gAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGACTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGc             >  1:976338/1‑150 (MQ=255)
                                                                                                                                         ggCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGt          <  1:1444011/150‑1 (MQ=255)
                                                                                                                                                  ggTTCGGCTTTTTGTATTTTTGCGTATCGCCTCGCAGTAAATGCCTTTCCAACTGGACGTTTTTACAGCAGAATTGTATTTTGTGCGGGTAAGTTGTTGCGTCAGGATGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGt  <  2:1631483/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGTAGCGGGTGTTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGT  >  NC_000913/4017075‑4017369

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: