Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,100,233 G→A intergenic (+39/‑63) ycdY → / → ycdZ redox enzyme maturation protein (REMP) chaperone for YcdX/DUF1097 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,100,2330GA93.0% 175.2 / 3.5 57intergenic (+39/‑63)ycdY/ycdZredox enzyme maturation protein (REMP) chaperone for YcdX/DUF1097 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base A (21/32);  total (23/34)
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.89e-01

CCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACA  >  NC_000913/1100139‑1100325
                                                                                              |                                                                                            
cccGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATttattgtt                                                                                         >  1:1144664/1‑100 (MQ=255)
   gcggtGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCa                                                                                      <  1:982751/96‑1 (MQ=255)
   gcgATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCa                                                                                      <  1:608695/100‑1 (MQ=255)
    cgATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCAt                                                                                     >  1:387544/1‑100 (MQ=255)
    cgATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCAt                                                                                     <  1:502192/100‑1 (MQ=255)
     gATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATa                                                                                    <  2:1037655/100‑1 (MQ=255)
         ccATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAAt                                                                                <  1:1102721/100‑1 (MQ=255)
         ccATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAAt                                                                                <  2:874412/100‑1 (MQ=255)
          cATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAAtg                                                                               >  2:977768/1‑100 (MQ=255)
          cATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAAtg                                                                               <  2:517179/100‑1 (MQ=255)
          cATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAAtg                                                                               >  1:828839/1‑100 (MQ=255)
             gatctcAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAAtgtgt                                                                            <  1:513657/95‑1 (MQ=255)
                tGCCATGAGGGACGATCTGGTAGCAGATTCTGTAGAATAAGTGTGATCTACGTGACTCATAACTGCGACGTATAATGTATTGTTGCGTgaaa                                                                                 >  2:634622/1‑92 (MQ=255)
                 gCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTc                                                                        >  2:356588/1‑100 (MQ=255)
                  cAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAAtgtgtagatc                                                                       >  2:513653/1‑95 (MQ=255)
                     tgtgGGACGAGCTGGAAGAATAGTCTGAATAATAAGTGTGATCTACGTCACGCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGa                                                                      >  2:760642/1‑98 (MQ=255)
                     tgtgGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGAtc                                                                    <  1:385638/100‑1 (MQ=255)
                     tgtgGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGAtc                                                                    <  2:375247/100‑1 (MQ=255)
                      gtgGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCCCTCATACCTGCACCGGATAATTTATTGTTGCATAAAATGTGTGCTCGAt                                                                     <  2:960855/98‑1 (MQ=255)
                      gtgGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATGATTTATTGTTGCATAAAATGTGTGCTCGAtct                                                                   <  1:1235866/100‑1 (MQ=255)
                      gtgGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGAtct                                                                   <  1:272246/100‑1 (MQ=255)
                      gtgGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGAtct                                                                   <  1:321051/100‑1 (MQ=255)
                       tgGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGAtctc                                                                  >  2:573997/1‑100 (MQ=255)
                         ggACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCtcat                                                                >  1:943603/1‑100 (MQ=255)
                          gACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCCACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCtcatt                                                               >  2:526662/1‑100 (MQ=255)
                          gACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCtcatt                                                               >  2:219953/1‑100 (MQ=255)
                                  ggAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCcg                                                       >  1:985039/1‑100 (MQ=255)
                                         gaTTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTc                                                <  2:250558/100‑1 (MQ=255)
                                         gaTTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTc                                                <  2:1206261/100‑1 (MQ=255)
                                           ttCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTg                                              <  2:255976/100‑1 (MQ=255)
                                            tCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGc                                             <  2:1113688/100‑1 (MQ=255)
                                             cTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCt                                            >  2:902142/1‑100 (MQ=255)
                                               gaagaaTAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTAt                                          <  2:1038336/100‑1 (MQ=255)
                                                    aTAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATgn                                     <  2:266459/100‑2 (MQ=255)
                                                      aaGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATgcgc                                   <  2:907230/100‑1 (MQ=255)
                                                        gtgtGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCgg                                 >  1:31387/1‑100 (MQ=255)
                                                        gtgtGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCgg                                 >  2:656804/1‑100 (MQ=255)
                                                         tgtgATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCgg                                 <  2:828833/99‑1 (MQ=255)
                                                          gtgATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCa                               <  1:356591/100‑1 (MQ=255)
                                                            gATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATg                             <  1:463447/100‑1 (MQ=255)
                                                               cTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAAc                          <  1:980302/100‑1 (MQ=255)
                                                                tACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACa                         >  1:72147/1‑100 (MQ=255)
                                                                tACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACa                         >  1:855377/1‑100 (MQ=255)
                                                                  cGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATa                       >  2:601373/1‑100 (MQ=255)
                                                                      aCTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTtc                   <  1:109982/100‑1 (MQ=255)
                                                                       cTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTtct                  >  2:1102945/1‑100 (MQ=255)
                                                                        tCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTtctc                 <  2:250067/100‑1 (MQ=255)
                                                                         cATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTtctct                >  2:1090208/1‑100 (MQ=255)
                                                                          aTAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTtctctc               <  2:387540/100‑1 (MQ=255)
                                                                           tAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTctctct              >  2:533103/1‑100 (MQ=255)
                                                                           tAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTctctct              <  2:943595/100‑1 (MQ=255)
                                                                            aaCTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTctctct              <  2:1118670/99‑1 (MQ=255)
                                                                            aaCTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTa             <  2:718494/100‑1 (MQ=255)
                                                                             aCTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTctctctct            <  1:885128/100‑3 (MQ=255)
                                                                              cTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTAtt           <  2:6518/100‑1 (MQ=255)
                                                                                gCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGc         <  2:166337/100‑1 (MQ=255)
                                                                                     ggATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATca    >  1:853657/1‑100 (MQ=255)
                                                                                       aTAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATcaca  >  1:664660/1‑100 (MQ=255)
                                                                                              |                                                                                            
CCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACA  >  NC_000913/1100139‑1100325

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