Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,975,297 Δ1 bp 69.4% coding (153/1131 nt) wecE → TDP‑4‑oxo‑6‑deoxy‑D‑glucose transaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,975,2970G.69.4% 139.7 / 120.3 85coding (153/1131 nt)wecETDP‑4‑oxo‑6‑deoxy‑D‑glucose transaminase
Reads supporting (aligned to +/- strand):  ref base G (11/15);  new base . (32/27);  total (43/42)
Fisher's exact test for biased strand distribution p-value = 3.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.53e-01

CCATTTAACGCACCGCCGGTGGTGGGAACCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGT  >  NC_000913/3975151‑3975437
                                                                                                                                                  |                                                                                                                                            
ccATTTAACGCACCGCCGGTGGTGGGAACCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGAc                                                                                                                                            <  1:470143/149‑1 (MQ=255)
        cGCACCGCCGGCTGTCGGAACCGAACTCGACTCTATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTCCCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGc                                                                                                                                  <  1:2426/149‑1 (MQ=255)
         gCACCGCCGGTGGTGGGAGCCGAACTCGACTGTATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGGAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTGCGCGCCGTCCTGCa                                                                                                                                 >  1:663851/1‑149 (MQ=255)
            ccgccgGTGGGGGGAACCGAACTCGACTTTATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCg                                                                                                                              <  2:625873/149‑1 (MQ=255)
                  gtggtgGGAACCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGc                                                                                                                        >  1:304073/1‑149 (MQ=255)
                   tggtggGACCCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCg                                                                                                                         <  2:440314/149‑1 (MQ=255)
                    ggtggGAACCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTg                                                                                                                      <  1:830904/149‑1 (MQ=255)
                      tggGAACCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTg                                                                                                                      <  1:185966/149‑1 (MQ=255)
                      tggGAACCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGga                                                                                                                    <  1:155223/149‑1 (MQ=255)
                       gggAACCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTAC‑‑CCGCCGTCCTGCACCGCTTCGCTGgag                                                                                                                   >  2:668113/1‑149 (MQ=255)
                       gggAACCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACAGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGACCAACGTTTTGGCAGCGCCACAGTGTTACTGGCGCCGTCCTGCACCCCTTCGCTgg                                                                                                                     >  2:808214/1‑149 (MQ=255)
                        ggAACCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGgaga                                                                                                                  >  1:821595/1‑149 (MQ=255)
                        ggAACCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGACGGTTTTACCCGTCGCTGCCCGCAGTGGCTGGAGCAACGTTTTGGCAACGCCAAAGTGTTGCT‑‑CGCCGTCCCGCACCGCTTCGCGGgaga                                                                                                                  >  1:334978/1‑149 (MQ=255)
                        ggAACCGAACTCGACTATATGCAGTCGGCAATGGGCAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGgaga                                                                                                                  <  1:301580/149‑1 (MQ=255)
                         gAACCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTCCCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGCGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGGAGAt                                                                                                                 <  2:5695/149‑1 (MQ=255)
                          aaCCGAAGTCGACTATATGCAGTGGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGGAGATg                                                                                                                >  2:583427/1‑149 (MQ=255)
                          aaCCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGgaga                                                                                                                  <  2:175319/149‑1 (MQ=255)
                            ccGAACTCGACTAGATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGGAGATggc                                                                                                              <  1:56814/149‑1 (MQ=255)
                              gAACTCGACTATATGCGGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGGAGCTggaca                                                                                                            <  2:745004/149‑4 (MQ=255)
                                 cgcGACTATATGCAGTCGGCAACGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGGAGATGGCGgcgc                                                                                                         <  1:647359/147‑1 (MQ=255)
                                  tCGACTATATGCAGCCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCgct                                                                                                        >  2:820149/1‑149 (MQ=255)
                                  ccGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTCCCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGgcggcg                                                                                                          <  1:342415/148‑1 (MQ=255)
                                    gACTATATGCAGTCGGCGATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGACGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCgctgc                                                                                                      >  2:632063/1‑149 (MQ=255)
                                    gACTATATGCAGTCGGCAATGGGTAGAGGCAAACTGAGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTGCT‑‑CGCCGTCCTGCACCGCATCGCTGGAGATGGCGGCgctgc                                                                                                      >  2:268603/1‑149 (MQ=255)
                                    gACTATACGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACCCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCgct                                                                                                        <  2:536247/149‑1 (MQ=255)
                                     acttgatGCAGGCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCgctg                                                                                                       <  2:225975/144‑1 (MQ=255)
                                      cTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCggtgc                                                                                                      >  1:723400/1‑149 (MQ=255)
                                        atatGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGCCAGAGCCAAAATGTTACT‑‑CGCCGTCCTGCACCGCTTAGCTGGAGATGGCGGCgctgctgct                                                                                                  >  2:49052/1‑149 (MQ=255)
                                        acatGCAGTCGGCAATGGGTAGCGGCAACCTGTGTTGCGATGGCGGTTTTACCCGTCGCGGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTCCGCTGGAGATGGCGGCgctgctgct                                                                                                  <  1:340280/147‑1 (MQ=255)
                                         tacgCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTc                                                                                                 <  2:363107/146‑1 (MQ=255)
                                          atGCAGTCGGCAATGGGTAGCGGCTAACCGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCg                                                                                                >  1:639586/1‑149 (MQ=255)
                                          atGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCg                                                                                                >  2:330987/1‑149 (MQ=255)
                                          atGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCg                                                                                                >  2:730650/1‑149 (MQ=255)
                                                   gCAATGGGTAGCGGCAAACGGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGc                                                                                       <  2:766682/149‑1 (MQ=255)
                                                            aGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGGGGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATg                                                                              >  2:594692/1‑149 (MQ=255)
                                                               ggCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAg                                                                           >  2:324759/1‑149 (MQ=255)
                                                               ggCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGACTGGAGCACCGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCTATGCAg                                                                           >  2:556535/1‑149 (MQ=255)
                                                                  ataCTGTGTGGCGATGGCGGGTTTGCGCGTCGCTGCGAGCAGTGGCTGGAGCAACGTTTTGGCAGTGTGAAAGTGTTACT‑‑TGCCGTCCTGCATCGCTTCGCTGGAGGTGGCGGCGCTGCTGCTTGATTTCCAGCCTGGGGATGAAGTgg                                                                        >  2:340770/3‑148 (MQ=255)
                                                                  aaaCTGTGTGGCGATGGCGGTTTGACCCGTCGCGGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCCGCCTGGCGATGGAGCGa                                                                        >  2:656120/1‑149 (MQ=255)
                                                                   aaCTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGAt                                                                       <  1:176757/149‑1 (MQ=255)
                                                                   aaCTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGAt                                                                       <  1:656086/149‑1 (MQ=255)
                                                                    aCTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGa                                                                        <  1:171611/149‑1 (MQ=255)
                                                                      tgtgtgGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATc                                                                      <  2:733532/149‑1 (MQ=255)
                                                                          tgGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCg                                                                <  2:335023/149‑1 (MQ=255)
                                                                          tgGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGCGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCCCGATATCCAGCCTGGCGATGAAGTGATCATGCCg                                                                >  2:338523/1‑149 (MQ=255)
                                                                           ggCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCgctgc                                                                                                      >  1:570659/1‑110 (MQ=255)
                                                                           ggCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCgctgc                                                                                                      <  2:570659/110‑1 (MQ=255)
                                                                            tcGATGGCGGTTTTCCCCGCCGCTGCCAGCTGCGGCTGCAGCAACGTTTTGGCTGCGCCAAAGTGTTACCGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGTCTCTGCTGCTCGATCTCCAGCCTGGCGTGGAAGTGATCATGCCg                                                                <  2:104423/148‑1 (MQ=255)
                                                                             cGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGCTGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAGGAGATCATGCCGa                                                               >  2:489284/1‑149 (MQ=255)
                                                                               aTGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGc                                                             >  1:766666/1‑149 (MQ=255)
                                                                                 ggCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCGGGCGCTGAAGTGAGCATGCCGAGCTaca                                                         >  1:758905/1‑149 (MQ=255)
                                                                                  gCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTac                                                          >  1:242966/1‑149 (MQ=255)
                                                                                   cGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGCTGAAGTGATCATGCCGAGCTACAcc                                                       >  1:800633/1‑149 (MQ=255)
                                                                                    ggTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCt                                                      >  1:721721/1‑149 (MQ=255)
                                                                                     gTTTTACCCGTCGCGGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCtt                                                     <  1:11493/149‑1 (MQ=255)
                                                                                        ttccccGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGGTACT‑‑CGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGCTATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGt                                                  <  2:97990/146‑1 (MQ=255)
                                                                                        ttACCCGGCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGCCGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATTTTCTGTTTGGCGATGCAGTGATCATGGCGAGCTACACCttt                                                    <  2:32050/149‑1 (MQ=255)
                                                                                             cgttcgtGCCAGCGGGGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCa                                             <  2:466331/143‑1 (MQ=255)
                                                                                             cGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCa                                             <  2:607486/149‑1 (MQ=255)
                                                                                              gTCGCGTCCAGCAGGGGCTGGAGCACCGGTTTGGCGGCGCCAAGGTGTCACT‑‑CGCCGTCCTGCACCGCTTCGCGGGAGATGGCGGCGCTGCTGCTCGATTTCCAGCCTGACGATGAAGTGACCATGCCGAGCTACACCTTTGTCTCCAc                                            <  2:371812/149‑1 (MQ=255)
                                                                                                cGCTGCGCGCAGTGGCTGGGCCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCg                                          <  2:178613/149‑1 (MQ=255)
                                                                                                cGCTGCCAGCAGGGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCg                                          <  2:415304/149‑1 (MQ=255)
                                                                                                 gCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCAGGCGTTGGAGGGATGATGCCGAGCTACACCTTTGTCTCCACCGc                                         >  1:582218/1‑149 (MQ=255)
                                                                                                       agcagTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCCATGc                                   >  1:240808/1‑149 (MQ=255)
                                                                                                        gcagTGGCTGGAGCAACGTTTTGGCAGCGCCCAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGGTGAAGTGACTATGCAGAGCTGCACCTTCGTCTCCACCGCCGATGcc                                  >  2:96865/1‑149 (MQ=255)
                                                                                                        gcagTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATg                                    >  1:518516/1‑149 (MQ=255)
                                                                                                         cagTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCt                                 <  1:326769/149‑1 (MQ=255)
                                                                                                             ggCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACTGCCAATGCCttt                               <  1:617355/149‑1 (MQ=255)
                                                                                                             ggCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCttt                               <  1:440603/149‑1 (MQ=255)
                                                                                                             ggCTGGAGCAACGTTTTGGCAGCGCCAAAGTCTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCttt                               <  1:133325/149‑1 (MQ=255)
                                                                                                              gCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTGCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTtgtg                            >  2:630136/1‑149 (MQ=255)
                                                                                                              gCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCGTCGCTGGAGAAGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGGTCCACCTGGGTCTCCACGGCCAATGCCTTggtg                            >  2:289264/1‑149 (MQ=255)
                                                                                                              ccTGGAGCAACGTTTTGGCAGCGCCTAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGGGCTACACCTTTGTCTCCACCGCCAATGCCTTtgtg                            <  2:5639/148‑1 (MQ=255)
                                                                                                                   aGCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCg                       <  1:90254/149‑1 (MQ=255)
                                                                                                                    gCAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGGGATCATGCCGGGCTCCACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGt                      >  2:403611/1‑149 (MQ=255)
                                                                                                                     cAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGGGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCg                       >  2:204941/1‑149 (MQ=255)
                                                                                                                     cAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGACATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCACTGCCTTTGtgctgcc                       >  2:1749/1‑148 (MQ=255)
                                                                                                                     cAACGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCCCCTCGCTGGACCTGGCGGCGCTGCTGCCCGATATCCAGCCTGGCCATGAAGCGATCATGCCGAGCTACACCTTTGCCACCACCGCCAATGCCTGCGTGCcgcgcg                     >  2:507574/1‑147 (MQ=255)
                                                                                                                      aaCGTTTTGGCAGCGCCAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATGCAGACCGGTTACGAAGCGTTCATGAAGCGCCACACCTTTGTCTCCACCGTCAATGCCTTTGTGCTGCGTgg                    <  2:821595/149‑1 (MQ=255)
                                                                                                                               gCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCaaaa             >  1:113535/1‑149 (MQ=255)
                                                                                                                               gCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTAGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGGAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAGTGCCTTTGTTGTGCGTGGCGCaaaa             >  2:222015/1‑149 (MQ=255)
                                                                                                                               gCAGCGCCAAAGTGTCACGCGCCGTCCTGCACCGCTTCGCTGTAGATGGCGGCGCTGCTGCTCGAGCTCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACATTTGTCTCAATCGTCAATGCCTTTGTGCTGCGTGGCGCAAAAAt           <  2:187893/149‑1 (MQ=255)
                                                                                                                                     ccAAAGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTCCACCTTTGTCTCCACCGCCCATGCCTTTGTCCTGCGTGGCGCAAACATCGtctt     >  2:572160/1‑149 (MQ=255)
                                                                                                                                      cAAAGTGTTACT‑‑CGCCGTCCTGCACCCCTTCGCTGGAGCTGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCCCGCTCCGCCTTTGTCTCCACCGCCAATGCCTCTGTGCTGCGCGGCGCAAAAATCGttctt    >  1:257684/1‑149 (MQ=255)
                                                                                                                                        aaGTGTTACT‑‑CGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTtgt  >  1:829141/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
CCATTTAACGCACCGCCGGTGGTGGGAACCGAACTCGACTATATGCAGTCGGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGT  >  NC_000913/3975151‑3975437

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: