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.1% 179.8 / 3.5 58intergenic (+39/‑63)ycdY/ycdZredox enzyme maturation protein (REMP) chaperone for YcdX/DUF1097 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (3/1);  new base A (30/24);  total (33/25)
Fisher's exact test for biased strand distribution p-value = 6.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

GCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACA  >  NC_000913/1100142‑1100328
                                                                                           |                                                                                               
gcgATGCCATTAGTGCAATGTGGGACGAGCTGAAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGc                                                                                          <  1:864047/99‑1 (MQ=255)
    tGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATaaa                                                                                     >  2:660565/1‑100 (MQ=255)
      ccATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAAt                                                                                   <  1:307948/100‑1 (MQ=255)
       cATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAAtg                                                                                  >  2:381872/1‑100 (MQ=255)
       cATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAAtg                                                                                  >  1:645135/1‑100 (MQ=255)
         ttAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAAtgtg                                                                                <  2:1451484/100‑1 (MQ=255)
         ttAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAAtgtg                                                                                <  2:1449468/100‑1 (MQ=255)
         ttAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAAtgtg                                                                                <  1:758120/100‑1 (MQ=255)
         tcaGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAAtgtg                                                                                <  1:812240/98‑1 (MQ=255)
          tAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAAtgtgt                                                                               <  1:1342862/100‑1 (MQ=255)
             gggaaTGTGGGAGGGGCTGGAAGAAGATTTTGAAGAATAAGTGTGGTATACGTCACTCATAGCTGCGACTTATGATTTGTTGGGGCATAAAAtgtg                                                                                <  2:732948/93‑1 (MQ=255)
              gCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTc                                                                           >  2:913480/1‑100 (MQ=255)
              gCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTc                                                                           >  1:1083867/1‑100 (MQ=255)
                 aTGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGAt                                                                        >  1:1358369/1‑100 (MQ=255)
                  tgtgGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGAtc                                                                       <  2:148973/100‑1 (MQ=255)
                  tgtgGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGAtc                                                                       <  1:620146/100‑1 (MQ=255)
                   gtgGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGAtct                                                                      <  1:1273634/100‑1 (MQ=255)
                   gtgGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGAtct                                                                      >  2:1063459/1‑100 (MQ=255)
                      ggACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCtcat                                                                   >  2:1207167/1‑100 (MQ=255)
                      ggACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCtcat                                                                   <  1:1159320/100‑1 (MQ=255)
                      ggACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCtcat                                                                   <  2:191474/100‑1 (MQ=255)
                            gCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATgg                                                             <  2:504414/100‑1 (MQ=255)
                                  agaagaTTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGt                                                       >  1:787880/1‑100 (MQ=255)
                                  agaagaTTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGt                                                       >  2:1462198/1‑100 (MQ=255)
                                   gaagaTTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGtt                                                      >  2:440171/1‑100 (MQ=255)
                                     agaTTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGtttt                                                    >  2:564031/1‑100 (MQ=255)
                                      gaTTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTc                                                   <  1:564034/100‑1 (MQ=255)
                                         tCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGc                                                <  1:1093235/100‑1 (MQ=255)
                                          cTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCt                                               >  2:1246762/1‑100 (MQ=255)
                                          cTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCt                                               >  2:154698/1‑100 (MQ=255)
                                          cTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCt                                               >  1:412976/1‑100 (MQ=255)
                                            gaagaaTAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTAt                                             <  1:904381/100‑1 (MQ=255)
                                            gaagaaTAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTAt                                             >  2:177900/1‑100 (MQ=255)
                                            gaagaaTAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTAt                                             >  2:104840/1‑100 (MQ=255)
                                            gaagaaTAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTAt                                             >  2:1070967/1‑100 (MQ=255)
                                            gaagaaTAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTAt                                             >  2:855417/1‑100 (MQ=255)
                                               gaaTAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCAt                                          >  1:760874/1‑100 (MQ=255)
                                                     gtgtGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCgg                                    <  1:1545119/100‑1 (MQ=255)
                                                         gATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATg                                >  2:1411266/1‑100 (MQ=255)
                                                          aTCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGa                               <  1:1078138/100‑1 (MQ=255)
                                                           tCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGaa                              >  1:9227/1‑100 (MQ=255)
                                                            cTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAAc                             <  1:148947/100‑1 (MQ=255)
                                                             tACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACa                            >  1:438570/1‑100 (MQ=255)
                                                                 tCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACt                        >  2:617939/1‑100 (MQ=255)
                                                                 tCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACt                        >  1:587301/1‑100 (MQ=255)
                                                                    cTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTtct                     >  2:782093/1‑100 (MQ=255)
                                                                      cATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTtctct                   >  1:310050/1‑100 (MQ=255)
                                                                           cTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTAtt              >  1:383917/1‑100 (MQ=255)
                                                                            tGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTg             <  1:1403914/100‑1 (MQ=255)
                                                                            tGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTg             <  1:1127353/100‑1 (MQ=255)
                                                                             gCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGc            <  1:514460/100‑1 (MQ=255)
                                                                              cAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCa           >  1:1194659/1‑100 (MQ=255)
                                                                              cAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCa           >  1:1048094/1‑100 (MQ=255)
                                                                                  ggATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATca       >  1:18763/1‑100 (MQ=255)
                                                                                    aTAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATcaca     <  1:1357104/100‑1 (MQ=255)
                                                                                      aaTTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCacaac   >  1:40762/1‑100 (MQ=255)
                                                                                      aaTTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCacaac   <  2:79701/100‑1 (MQ=255)
                                                                                       aTTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCacaaca  >  2:793930/1‑100 (MQ=255)
                                                                                           |                                                                                               
GCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACA  >  NC_000913/1100142‑1100328

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: