Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,336,360 A→G 7.7% L79P (CTG→CCG)  yghU ← GSH‑dependent disulfide bond oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,336,3600AG 7.7% 144.0 / 2.6 52L79P (CTG→CCG) yghUGSH‑dependent disulfide bond oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (20/28);  new base G (2/2);  total (22/30)
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.23e-01

TTCTCTGCCAGGTACAGCAGGATCGAGCCTGACTCGAACACCCGCACCGGCGGCTCCACGCTGCGGTCGAGCAGGGCAGGGATCTTCGAGTTGGGGTTGACCTGGACGAAGCCACTGGAAAACTGATCACCCTCGCCGATGCGGATCAGCCAGGCATCGTATTCGGCACCGGCGTGGCCGAGGGCAAGCAGCTCTTCAAGGGCAATGGTGACCTTGACGCCATTGGGCGTGGCCAGGGAATACAGCTGCAGCGGGTGCTTGCCCAGCGGCAGGTCCTTGTCGTGAGTCGGGCCGG  >  NC_002947/3336213‑3336507
                                                                                                                                                   |                                                                                                                                                   
ttCGCTGACAGGTACAGCAGGATCGAGCCTGACTCGAACACCCGCACCGGCGGCTCCACGCTGCGGTCGAGCAGGGCAGGGATCTTAGAGTTGGGGTTGACCTGGACGAAGCCACTGGAAAACTGATCACCCTCGCCGATGCGGATCAg                                                                                                                                                    >  2:70513/1‑149 (MQ=255)
  ctctGCCAGGTACAGCAGGATCGAGCCTGACT‑‑AACACCCGCACCGGCGGCTCCACGCTGCGGTCGAGCAGGGCAGGGATCTTCGAGTTGGGGTTGCCCTGGACGAAGCCACTGGAAAACTGATCACCCTCGCCGATGCGGATCAGCCAg                                                                                                                                                >  6:190410/1‑149 (MQ=255)
              cagcagGATCGAGCCTGACTCGAACACCCGCACCGGCGGCTCCACGCTGCGGTCGAGCAGGGCAGGGATCTTCGAGTTGGGGTTGACCTGGACGAAGCCACTGGAAAACTGATCACCCTCGCCGATGCGGATCAGCCAGGCATCGTAtt                                                                                                                                      >  4:157599/1‑149 (MQ=255)
               agcagGATCGAGCCTGACTCGAACACCCGCACCGGCGGCTCCACGCTGCGGTCGAGCAGGGCAGGGATCTTCGAGTTGGGGTTGACCTGGACGAAGCCACTGGAAAACTGATCACCCTCGCCGATGCGGATCAGCCAGGCATCGTATTc                                                                                                                                     <  1:160769/149‑1 (MQ=255)
               agcagGATCGAGCCTGACTCGAACACCCGCACCGGCGGCTCCACGCTGCGGTCGAGCAGGGCAGGGATCTTCGAGTTGGGGTTGACCTGGACGAAGCCACTGGAAAACTGATCACCCTCGCCGATGCGGATCAGCCAGGCATCGTATTc                                                                                                                                     >  1:144121/1‑149 (MQ=255)
                 cagGATCGAGCCTGACTCGAACACCCGCACCGGCGGCTCCACGCTGCGGTCGAGCAGGGCAGGGATCTTCGAGTTGGGGTTGACCTGGACGAAGCCACTGGAAAACTGATCACCCTCGCCGATGCGGATCAGCCAGGCATCGTATTCgg                                                                                                                                   >  1:217892/1‑149 (MQ=255)
                   ggATCGAGCCTGACTCGAACACCCGCACCGGCGGCTCCACGCTGCGGTCGAGCAGGGCAGGGATCTTCGAGTTGGGGTTGACCTGGACGAAGCCACTGGAAAACTGATCACCCTCGCCGATGCGGATCa                                                                                                                                                     <  8:65766/129‑1 (MQ=255)
                   ggATCGAGCCTGACTCGAACACCCGCACCGGCGGCTCCACGCTGCGGTCGAGCAGGGCAGGGATCTTCGAGTTGGGGTTGACCTGGACGAAGCCACTGGAAAACTGATCACCCTCGCCGATGCGGATCa                                                                                                                                                     >  7:65766/1‑129 (MQ=255)
                        gAGCCTTTCTCGAACACCCGCACCGGCGGCTCCAATCTGCTGTCGATCATGGCAGGGATATTCGAGTTGGGTTTGACCTTGACTAATACACTGGAAAACTGATCACCCTCTCCTATTCGTATCAGCCAGTCTTCTTATTCGGCACCGGc                                                                                                                            <  4:46578/149‑1 (MQ=255)
                         aGCCTGACTCGAACACCCGCACCGGCGGCTCCACGCTGCGGTCGAGCAGGGCAGGGATCTTCGAGTTGGGGTTGACCTGGACGAAGCCACTGGAAAACTGATCACCCTCGCCGATGCGGATCAGCCAGGCATCGTATTCGGCACCGGCg                                                                                                                           >  4:316475/1‑149 (MQ=255)
                                cTCGAACACCCGCACCGGCGGCTCCACGCTGCGGTCGAGCAGGGCAGGGATCTTCGAGTTGGGGTTGACCTGGACGAAGCCACTGGAAAACTGATCACCCTCGCCGATGCGGATCAGCCAGGCATCGTATTCGGCACCGGCGTGGCCGa                                                                                                                    <  2:217892/149‑1 (MQ=255)
                                 tCGAACACCCGCACCGGCGGCTCCACGCTGCGGTCGAGCAGGGCAGGGATCTTCGAGTTGGGGTTGACCTGGACGAAGCCACTGGAAAACTGATCACCCTCGCCGATGCGGATCAGCCAGGCATCGTATTCGGCACCGGCGt                                                                                                                          <  8:286246/142‑1 (MQ=255)
                                 tCGAACACCCGCACCGGCGGCTCCACGCTGCGGTCGAGCAGGGCAGGGATCTTCGAGTTGGGGTTGACCTGGACGAAGCCACTGGAAAACTGATCACCCTCGCCGATGCGGATCAGCCAGGCATCGTATTCGGCACCGGCGt                                                                                                                          >  7:286246/1‑142 (MQ=255)
                                 tCGAACACCCGCACCGGCGGCTCCACGCTGCGGTCGAGCAGGGCAGGGATCTTCGAGTTGGGGTTGACCTGGACGAAGCCACTGGAAAACTGATCACCCTCGCCGATGCGGATCAGCCAGGCATCGTATTCGGCACCGGCGTGGCCGAg                                                                                                                   >  1:47322/1‑149 (MQ=255)
                                  cGAACACCCGCACCGGCGGCTCCACGCTGCGGTCGAGCAGGGCAGGGATCTTCGAGTTGGGGTTGACCTGGACGAAGCCACTGGAAAACTGATCACCCTCGCCGATGCGGATCAGCCAGGCATCGTATTCGGCACCGGCGTGGCCGAgg                                                                                                                  <  3:122314/149‑1 (MQ=255)
                                   gAACACCCGCACCGGCGTCTCCACGCTGCGGTCGAGCAGGGCAGGTATCTTCGAGTTGGGGTTGACCTGGACGAAGCCACTGGAAAACTGATCACCCTCGCCGATGCGGATCAGCCAGGCATCGTATTCGGCACCGGCGTGGCCGAgg                                                                                                                  <  4:112965/148‑1 (MQ=255)
                                   gAACACCCGCACCGGCGGCTCCACGCTGCGGTCGAGCAGGGCAGGGATCTTCGAGTTGGGGTTGACCTGGACGAAGCCACTGGAAAACTGATCACCCTCGCCGATGCGGATCAGCCAGGCATCGTATTCGGCACCGGCGTGGCCGAggg                                                                                                                 <  4:210333/149‑1 (MQ=255)
                                   gAACACCCGCACCGGCGGCTCCACGCTGCGGTCGAGCAGGGCAGGGATCTTCGAGTTGGGGTTGACCTGGACGAAGCCACTGGAAAACTGATCACCCTCGCCGATGCGGATCAGCCAGGCATCGTATTCGGCACCGGCGTGGCCGAgg                                                                                                                  <  3:117572/148‑1 (MQ=255)
                                    aaCACCCGCACCGGCGGCTCCACGCTGCGGTCGAGCAGGGCAGGGATCTTCGAGTTGGGGTTGACCTGGACGAAGCCACTGGAAAACTGATCACCCTCGCCGATGCGGATCAGCCAGGCATCGTATTCGGCACCGGCGTGGCCGAGGGc                                                                                                                <  4:256016/149‑1 (MQ=255)
                                     acacCCGCACCGGCGGCTCCACGCTGCGGTCGAGCAGGGCAGGGATCTTCGAGTTGGGGTTGACCTGGACGAAGCCACTGGAAAACTGATCACCCTCGCCGATGCGGATCAGCCAGGCATCGTATTCGGCACCGGCGTGGCCGAGGGc                                                                                                                >  5:239489/1‑148 (MQ=255)
                                     acacCCGCACCGGCGGCTCCACGCTGCGGTCGAGCAGGGCAGGGATCTTCGAGTTGGGGTTGACCTGGACGAAGCCACTGGAAAACTGATCACCCTCGCCGATGCGGATCAGCCAGGCATCGTATTCGGCACCGGCGTGGCCGAGGGCa                                                                                                               <  3:157599/149‑1 (MQ=255)
                                        cccGCACCGGCGGCTCCACGCTGCGGTCGAGCAGGGCAGGGATCTTCGAGTTGGGGTTGACCTGGACGAAGCCACTGGAAAACTGATCACCCTCGCCGATGCGGATCAGCCAGGCATCGTATTCGGCACCGGCGTGGCCGAGGGCAagc                                                                                                            <  8:189678/149‑1 (MQ=255)
                                         ccGCACCGGCGGCTCCACGCTGCGGTCGAGCAGGGCAGGGATCTTCGAGTTGGGGTTGACCTGGACGAAGCCACTGGAAAACTGATAACCCTCGCCGATGCGGATCAGCCAGGCATCGTAGTCGGCACCGGCGTGGACGAGGGCCag                                                                                                             >  8:281548/1‑147 (MQ=255)
                                             aCCGGCGGCTCCACGCTGCGGTCGAGCAGGGCAGGGATCTTCGAGTTGGGGTTGACCTGGACGAAGCCACTGGAAAACTGATCACCCTCGCCGATGCGGATCAGCCAGGCATCGTATTCGGCACCGGCGTGGCCGAGGGCAAGCAGctc                                                                                                       <  3:316475/149‑1 (MQ=255)
                                                             tGCGGTCGAGCAGGGCAGGGATCTTCGAGTTGGGGTTGACCTGGACGAAGCCACTGGAAAACTGATCACCCTCGCCGATGCGGATCAGCCAGGCATCGTATTCGGCACCGGCGTGGCCGAGGGCAAGCAGCTCTTCAAGGGCAATGGt                                                                                        >  2:94375/1‑148 (MQ=255)
                                                               cGGTCGAGCAGGGCAGGGATCTTCGAGTTGGGGTTGACCTGGACGAAGCCACTGGAAAACTGATCACCCTCGCCGATGCGGATCAGCCAGGCATCGTATTCGGCACCGGCGTGGCCGAGGGCAAGCAGCTCTTCAAGGGCAATGGTGAc                                                                                     <  6:239489/149‑1 (MQ=255)
                                                                 gTCGAGCAGGGCAGGGATCTTCGAGTTGGGGTTGACCTGGACGAAGCCACTGGAAAACTGATCACCCTCGCCGATGCGGATCAGCCAGGCATCGTATTCGGCACCGGCGTGGCCGAGGGCAAGCAGCTCTTCAAGGGCAATGGTGACCt                                                                                   <  1:42440/149‑1 (MQ=255)
                                                                    gAGCAGGGCAGGGATCTTCGAGTTGGGGTTGACCTGGACGAAGCCACTGGAAAACTGATCACCCTCGCCGATGCGGATCAGCCAGGCATCGTATTCGGCACCGGCGTGGCCGAGGGCAAGCAGCTCTTCAAGGGCAATGGTGACCTTg                                                                                 <  2:144121/148‑1 (MQ=255)
                                                                          ggcagggATCTTCGAGTTGGGGTTGACCTGGACGAAGCCACTGGAAAACTGATCACCCTCGCCGATGCGGATCAGCCAGGCATCGTATTCGGCACCGGCGTGGCCGAGGGCAAGCAGCTCTTCAAGGGCAATGGTGACCTTGACGCCAt                                                                          >  3:127388/1‑149 (MQ=255)
                                                                                    ttCGAGTTGGGGTTGACCTGGACGAAGCCACTGGAAAACTGATCACCCTCGCCGATGCGGATCAGCCAGGCATCGTATTCGGCACCGGCGTGGCCGAGGGCAAGCAGCTCTTCAAGGGCAATGGTGACCTTGACGCCATTGGGCGTGGc                                                                >  5:31870/1‑149 (MQ=255)
                                                                                      cGAGTTGGGTGAGACCTGTACGAAGCCACTGGAAAACTCATCACCCTCGCCGATGCGGATCAGCCAGGCATCGAATTCGGCACCGGCGTGGCCGAGGGCAAGCAGCTCTTCAAGGGCAATGGTGACCTTGACGCCATTGGGCGTg                                                                  <  5:248035/145‑1 (MQ=255)
                                                                                       gAGTTGGGGTTGACCTGGACGAAGCCACTGGAAAACTGATCACCCTCGCCGATGCGGATCAGCCAGGCATCGTATTCGGCACCGGCGTGGCCGAGGGCAAGCAGCTCTTCAAGGGCAATGGTGACCTTGACGCCATTGGGCGTGGCCa                                                              <  8:47063/148‑1 (MQ=255)
                                                                                        aGTTGGGGTTGACCTGGACGAAGCCACTGGAAAACTGATCACCCTCGCCGATGCGGATCAGCCAGGCATCGTATTCGGCACCGGCGTGGCCGAGGGCAAGCAGCTCTTCAAGGGCAATGGTGACCTTGACGCCATTGGGCGTGGCCAgg                                                            <  8:234630/149‑1 (MQ=255)
                                                                                        aGTTGGGGTGGACCTGGACGAAGCCACTGGAAAACTGATCACCCTCGCCGATGCGGATCAGCCAGGCATCGTATTCGGCACCGGCGTGGCCGAGGGCAAGCAGCTCTTCAAGGGCAATGGTGACCTTGACGCCATTGGGCGTGGCCAgg                                                            <  4:234879/149‑1 (MQ=255)
                                                                                              ggTTGACCTGGACGAAGCCACTGGAAAACTGATCACCCTCGCCGATGCGGATCAGCCAGGCATCGTATTCGGCACCGGCGTGGCCGAGGGCAAGCAGCTCTTCAAGGGCAATGGTGACCTTGACGCCATTGGGCGTGGCCAGGGAATa                                                       >  1:251896/1‑148 (MQ=255)
                                                                                                ttGACCTGGACGAAGCCACTGGAAAACTGATCACCCTCGCCGATGCGGATCAGCCAGGCATCGTATTCGGCACCGGCGTGGCCGAGGGCAAGCAGCTCTTCAAGGGCAATGGTGACCTTGACGCCATTGGGCGTGGCCAgg                                                            <  1:74529/141‑1 (MQ=255)
                                                                                                ttGACCTGGACGAAGCCACTGGAAAACTGATCACCCTCGCCGATGCGGATCAGCCAGGCATCGTATTCGGCACCGGCGTGGCCGAGGGCAAGCAGCTCTTCAAGGGCAATGGTGACCTTGACGCCATTGGGCGTGGCCAgg                                                            >  2:74529/1‑141 (MQ=255)
                                                                                                ttGACCTGGACGAAGCCACTGGAAAACTGATCACCCTCGCCGATGCGGATCAGCCAGGCATCGTATTCGGCACCGGCGTGGCCGAGGGCAAGCAGATCTTCTAGGGCAATGGTGACCTTGACGCCATTGGGCGTGGCCAGGGAATACAg                                                    <  1:70513/149‑1 (MQ=255)
                                                                                                  gACCTGGACGAAGCCACTTGAAAACTGATCACCCTCGCCGATGCGGATCAGCCAGGCATCGTATTCGGCACCGGCGTGGCCGAGGGCAAGCAGCTCTTCAAGGGCAATGGTGACCTTGACGCCATTGGGCGTGGCCAGGGAATACAGc                                                   <  2:47322/148‑1 (MQ=255)
                                                                                                     cTGGACGAAGCCACTGGAAAACTGATCACCCTCGCCGATGCGGATCGGCCAGGCATCGTATTCGGCACCGGCGTGGCCGAGGGCAAGCAGCTCTTCAAGGGCAATGGTGACCTTGACGCCATTGGGCGTGGCCAGGGAATACAGCTg                                                 >  4:10417/1‑147 (MQ=255)
                                                                                                     cTGGACGAAGCCACTGGAAAACTGATCACCCTCGCCGATGCGGATCAGCCAGGCATCGTATTCGGCACCGGCGTGGCCGAGGGCAAGCAGCTCTTCAAGGGCAATGGTGACCTTGACGCCATTGGGCGTGGCCAGGGAATACAGcag                                                 >  7:290881/1‑145 (MQ=255)
                                                                                                     cTGGACGAAGCCACTGGAAAACTGATCACCCTCGCCGAAGCGGATCGGCCAGGCATCGTATTCGGCACCGGCGTGGCCGAGGGCAAGCAGCTCTTCAAGGGCAATGTTGACCTTGACGCCATTGGCCGTGGACAGGGAATACAGCTGCa                                               >  4:14906/1‑149 (MQ=255)
                                                                                                           gAAGCCACTGGAAAACTGATCACCCTCGCCGATGCGGATCAGCCAGGCATCGTATTCGGCACCGGCGTGGCCGAGGGCAAGCAGCTCTTCAAGGGCAATGGTGACCTTGACGCCATTGGGCGTGGCCAGGGAATACAGCTGCAGCGGGt                                         >  5:87147/1‑149 (MQ=255)
                                                                                                               ccACTGGAAAACTGATCACCCTCGCCGATGCGGATCAGCCAGGCATCGTATTCGGCACCGGCGTGGCCGAGGGCAAGCAGCTCTTCAAGGGCAATGGTGACCTTGACGCCATTGGGCGTGGCCAGGGAATACAGCTGCAGCGGGTGCtt                                     >  8:279740/1‑149 (MQ=255)
                                                                                                                     gAAAACTGATCACCCTCGCCGATGCGGATCAGCCAGGCATCGTATTCGGCACCGGCGTGGCCGAGGGCAAGCAGCTCTTCAAGGGCAATGGTGACCTTGACGCCATTGGGCGTGGCCAGGGAATACAGCTGCAGCGGGTGCTTGCCCAg                               >  6:148366/1‑149 (MQ=255)
                                                                                                                            gATCACCCTCGCCGATGCGGATCAGCCAGGCATCGTATTCGGCACCGGCGTGGCCGAGGGCAAGCAGCTCTTCAAGGGCAATGGTGACCTTGACGCCATTGGGCGTGGCCAGGGGATACAGCTGCAGCGGGTGCTTGCCCAGCGGCAgg                        <  8:305926/149‑1 (MQ=255)
                                                                                                                                aCCCTCGCCGATGCGGATCAGCCAGGCATCGTATTCGGCACCGGCGTGGCCGAGGGCAAGCAGCTCTTCAAGGGCAATGGTGACCTTGACGCCATTGGGCGTGGCCAGGGAATACAGCTGCAGCGGGTGCTTGCCCAGCGGCAGGTcc                     <  6:31870/148‑1 (MQ=255)
                                                                                                                                 cccTCGCCGATGCGGATCAGCCAGGCAGCGTATTCGGCACCGGCGTGGCCGAGGGCAAGCAGCTCTTCAAGGGCAATGGTGACCTTGACGCCATTGGGCGTGACCAGGGAATACAGCTGCAGCGGGTGCTTGCCCAGCGGCAGGTCCtt                   <  7:153242/149‑1 (MQ=255)
                                                                                                                                      gCCGATGCGGATCAGCCAGGCATCGTATTCGGCACCGGCGTGGCCGAGGGCAAGCAGCTCTTCAAGGGCAATGGTGACCTTGACGCCATTGGGCGTGGCCAGGGAATACAGCTGCAGCGGGTGCTTGCCCAGCGGCAGGTCCTTGTCGt              <  7:279740/149‑1 (MQ=255)
                                                                                                                                        cGATGCGGATCAGCCAGGCATCGTATTCGGCACCGGCGTGGCCGAGGGCAAGCAGCTCTTCAAGGGCAATGGTGACCTTGACTCCATTGGGCGTGGCCAGGGAATACAGCTGCAGCGGGTGCTTGCCCAGCGGCAGGTCCTTGTCGTGa            <  2:251896/149‑1 (MQ=255)
                                                                                                                                        cGATGCGGATCAGCCAGGCATCGTATTCGGCACCGGCGTGGCCGAGGGCAAGCAGCTCTTCAAGGGCAATGGTGACCTTGACGCCATTGGGCGTGGCCAGGGAATACAGCTGCAGCGGGTGCTTGCCCAGCGGCAGGTCCTTGTCGTGa            <  4:30770/149‑1 (MQ=255)
                                                                                                                                           tGCGGATCAGCCAGGCATCGTATTCGGCACCGGCGTGGCCGAGGGCAAGTAGCTCTTCAAGGGCAATGGTGACCTTGACGCCATTGGGCGTGGCCAGGGAATACAGCTGCAGCGGGTGCTTGCCCAGCGGCAGGTCCTTGTCGTGAGTc         >  4:271787/1‑149 (MQ=255)
                                                                                                                                              ggATCGGCCAGGCATCGTATTCGTCACCGGCGTGGCCGAGGGCAAGCAGCTCTTCAAGGGCAATGGTGACCTTGACGCCATTGGGCGTGGCCAGGGAATACAGCTGCAGCGGGTGCTTGCCCAGCGGCAGGTCCTTGTCGTGAGTCggg      <  3:14906/149‑1 (MQ=255)
                                                                                                                                              ggATCGGCCAGGCATCGTATTCGGCACCGGCGTGGCCGAGGGCAAGCAGCTCTTCAAGGGCAATGGTGACCTTGACGCCATTGGGCGTGGCCAGGGAATACAGCTGCAGCGGGTGCTTGCCCAGCGGCAGGTCCTTGTCGTGAGTCggg      <  3:10417/149‑1 (MQ=255)
                                                                                                                                                   aGCCAGGCATCGTATTCGGCACCGGCGTGGCCGAGGGCAAGCAGCTCTTCAAGGGCAATGGTGACCTTGACGCCATTGGGCGTGGCCAGGGAATACAGCTGCAGCGGGTGCTTGCCCAGCGGCAGGTCCTTGTCGTGAGTCGGGCCgg  <  4:127388/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TTCTCTGCCAGGTACAGCAGGATCGAGCCTGACTCGAACACCCGCACCGGCGGCTCCACGCTGCGGTCGAGCAGGGCAGGGATCTTCGAGTTGGGGTTGACCTGGACGAAGCCACTGGAAAACTGATCACCCTCGCCGATGCGGATCAGCCAGGCATCGTATTCGGCACCGGCGTGGCCGAGGGCAAGCAGCTCTTCAAGGGCAATGGTGACCTTGACGCCATTGGGCGTGGCCAGGGAATACAGCTGCAGCGGGTGCTTGCCCAGCGGCAGGTCCTTGTCGTGAGTCGGGCCGG  >  NC_002947/3336213‑3336507

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