Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,856,453 C→T *118* (TAG→TAA ygbA ← uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,856,4530CT98.0% 155.6 / NA 50*118* (TAG→TAAygbAuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (25/24);  minor base A (1/0);  total (26/24)
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.99e-01

CCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCT  >  NC_000913/2856324‑2856583
                                                                                                                                 |                                                                                                                                  
ccTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGa                                                                                                                           >  1:1137473/1‑139 (MQ=255)
        cgcTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGtttttt                                                                                                                   <  1:1086100/139‑1 (MQ=255)
         gcTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTc                                                                                                                  >  2:657829/1‑139 (MQ=255)
         gcTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTc                                                                                                                  >  2:1035398/1‑139 (MQ=255)
          cTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCg                                                                                                                 >  2:1051640/1‑139 (MQ=255)
             cgcAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGAc                                                                                                              >  2:1384954/1‑139 (MQ=255)
                          gCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCgg                                                                                                 <  1:1035398/139‑1 (MQ=255)
                          gCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCgg                                                                                                 <  1:657829/139‑1 (MQ=255)
                            tGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCa                                                                                               >  1:1450579/1‑139 (MQ=255)
                              aaaCGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGt                                                                                             >  1:771165/1‑139 (MQ=255)
                                   cACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGa                                                                                                               >  2:525272/1‑116 (MQ=255)
                                   cACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGa                                                                                                               <  1:525272/116‑1 (MQ=255)
                                   cACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATAtttt                                                                                                      >  1:76032/1‑125 (MQ=255)
                                   cACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATAtttt                                                                                                      <  2:76032/125‑1 (MQ=255)
                                    aCGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGc                                                                                       >  2:323233/1‑139 (MQ=255)
                                          ctgctgTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGa                                                                                 <  2:771165/139‑1 (MQ=255)
                                          ctgctgTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGa                                                                                 <  1:1010433/139‑1 (MQ=255)
                                            gctgTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACg                                                                               <  2:1450579/139‑1 (MQ=255)
                                                  aCGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTc                                                                         >  2:1160188/1‑139 (MQ=255)
                                                  aCGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTc                                                                         <  1:846281/139‑1 (MQ=255)
                                                  aCGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTc                                                                         >  1:1282500/1‑139 (MQ=255)
                                                  aCGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTc                                                                         >  2:846281/1‑139 (MQ=255)
                                                  aCGGATGAAGCGACTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGa                                                                                 >  2:236015/1‑131 (MQ=255)
                                                  aCGGATGAAGCGACTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGa                                                                                 <  1:236015/131‑1 (MQ=255)
                                                    ggATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCAt                                                                       <  1:1384954/139‑1 (MQ=255)
                                                         aaGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCgatg                                                                  <  1:1051640/139‑1 (MQ=255)
                                                            cGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCgatgaga                                                               >  1:222878/1‑139 (MQ=255)
                                                                 gggAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACg                                                          <  2:1137473/139‑1 (MQ=255)
                                                                   gAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCa                                                        <  2:1282500/139‑1 (MQ=255)
                                                                            aaTCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGa                                               >  1:732074/1‑139 (MQ=255)
                                                                                   aTTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTg                                                                                         <  1:302566/90‑1 (MQ=255)
                                                                                   aTTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTg                                                                                         >  2:302566/1‑90 (MQ=255)
                                                                                   aTTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGc                                                         >  2:591924/1‑122 (MQ=255)
                                                                                   aTTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGc                                                         <  1:591924/122‑1 (MQ=255)
                                                                                                cAAATTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCAtt                           >  2:784815/1‑139 (MQ=255)
                                                                                                 aagttGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCATTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCcaca                                <  1:894914/130‑1 (MQ=255)
                                                                                                 aaaTTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCgatgagatg                                                             >  1:641208/1‑104 (MQ=255)
                                                                                                 aaaTTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCgatgagatg                                                             <  2:641208/104‑1 (MQ=255)
                                                                                                 aaaTTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCgat                                                                   >  2:1489038/1‑98 (MQ=255)
                                                                                                 aaaTTGCCGGACAGATCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCgat                                                                   <  1:1489038/98‑1 (MQ=255)
                                                                                                 aaaTTGCCGGACAGATCTGCCTGTCCGGCATAATATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCcaca                                >  2:894914/1‑133 (MQ=255)
                                                                                                               aTCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACg                                     >  1:604399/1‑114 (MQ=255)
                                                                                                               aTCTGCCTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACg                                     <  2:604399/114‑1 (MQ=255)
                                                                                                                     cTGTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATa      <  1:1160188/139‑1 (MQ=255)
                                                                                                                       gTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGa                                                                                 >  2:502844/1‑62 (MQ=255)
                                                                                                                       gTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGa                                                                                 <  1:502844/62‑1 (MQ=255)
                                                                                                                       gTCCGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAAt    <  1:323233/139‑1 (MQ=255)
                                                                                                                         ccGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGAt                                              >  2:1946/1‑95 (MQ=255)
                                                                                                                         ccGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGAt                                              <  1:1946/95‑1 (MQ=255)
                                                                                                                         ccGGCATATTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCt  >  2:95069/1‑139 (MQ=255)
                                                                                                                                 |                                                                                                                                  
CCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCT  >  NC_000913/2856324‑2856583

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