Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,984,802 G→A 7.4% A264V (GCA→GTA)  yqeF ← short chain acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,984,8020GA 7.4% 163.5 / 2.6 54A264V (GCA→GTA) yqeFshort chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (33/17);  new base A (2/2);  total (35/19)
Fisher's exact test for biased strand distribution p-value = 6.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

GATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCG  >  NC_000913/2984654‑2984948
                                                                                                                                                    |                                                                                                                                                  
gagaaGATCGCCTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGAGCAATTCCCATCAATGCCGGATCTACACCAACGCTGTCAAATGCGCGAATGCGGGCCAGCCCGGGTAAATTCAACGCTCGTg                                                                                                                                                    >  3:261201/4‑149 (MQ=255)
      aTCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTgg                                                                                                                                              <  1:23934/149‑1 (MQ=255)
       tCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCCTGCTTTGGc                                                                                                                                             >  5:262742/1‑149 (MQ=255)
          aCTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTc                                                                                                                                          >  6:138836/1‑149 (MQ=255)
                    aCTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCtcatcat                                                                                                                                 >  2:162075/1‑148 (MQ=255)
                        ccAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTa                                                                                                                            >  5:35760/1‑149 (MQ=255)
                            ccTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCg                                                                                                                         >  2:268729/1‑148 (MQ=255)
                             cTACACGCTCCAGGCAACGGCGGGTCTCATATATTTTATTTATTCGGATAAATGCGTGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTACTTTGGCTTCGCtcatc                                                                                                                                   <  8:120570/137‑1 (MQ=255)
                             cTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGGTCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTACTTTGGCTTCGCtcatc                                                                                                                                   <  2:61477/137‑1 (MQ=255)
                             cTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGGTCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTACTTTGGCTTCGCtcatc                                                                                                                                   >  1:61477/1‑137 (MQ=255)
                             cTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTACTTTGGCTTCGCtcatc                                                                                                                                   >  7:120570/1‑137 (MQ=255)
                              tACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAg                                                                                                                      >  2:116662/1‑149 (MQ=255)
                              tACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAg                                                                                                                      >  2:202286/1‑149 (MQ=255)
                                   gCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGAAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATAACCGCAGCTgcg                                                                                                                 >  5:205525/1‑149 (MQ=255)
                                                 cggGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGga                                                                                                   <  1:268729/149‑1 (MQ=255)
                                                 cggGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGga                                                                                                   >  6:181342/1‑149 (MQ=255)
                                                        cATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCa                                                                                            >  1:147869/1‑149 (MQ=255)
                                                          tACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGCTGATGCAt                                                                                           >  5:260592/1‑148 (MQ=255)
                                                            cacCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTAc                                                                                        <  5:181342/149‑1 (MQ=255)
                                                                 gcgcAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCt                                                                                   >  1:75526/1‑149 (MQ=255)
                                                                        tCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCtcatca                                                                                                                                  >  5:190734/1‑95 (MQ=255)
                                                                         cccATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAg                                                                            <  5:138836/148‑1 (MQ=255)
                                                                            aTCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAAcc                                                                        >  8:370081/1‑149 (MQ=255)
                                                                            aTCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAAcc                                                                        >  7:305596/1‑149 (MQ=255)
                                                                               aaTGCGTGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCga                                                                      <  2:272051/148‑1 (MQ=255)
                                                                                 tGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGACCCGAGAc                                                                   >  6:42739/1‑149 (MQ=255)
                                                                                    cGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTAt                                                                >  7:344620/1‑149 (MQ=255)
                                                                                      gATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATc                                                               >  7:316507/1‑148 (MQ=255)
                                                                                           acacCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATg                                                          >  1:86250/1‑148 (MQ=255)
                                                                                           acacCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATg                                                          >  8:173677/1‑148 (MQ=255)
                                                                                            cacCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGaa                                                        <  7:370081/149‑1 (MQ=255)
                                                                                             acCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAg                                                       >  6:154970/1‑149 (MQ=255)
                                                                                               cAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGa                                                     <  1:116662/149‑1 (MQ=255)
                                                                                                 aCGATGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCCTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGAtt                                                   <  6:262742/149‑1 (MQ=255)
                                                                                                     tGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTaaa                                               >  6:270365/1‑149 (MQ=255)
                                                                                                       gCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACg                                             <  5:270365/149‑1 (MQ=255)
                                                                                                          aaTGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGgatga                                                                                                >  1:60846/1‑95 (MQ=255)
                                                                                                          aaTGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGgatga                                                                                                <  2:60846/95‑1 (MQ=255)
                                                                                                              cgcgAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTaa                                       <  6:35760/148‑1 (MQ=255)
                                                                                                                    tGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCt                                >  1:187379/1‑149 (MQ=255)
                                                                                                                    tGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCt                                >  7:70275/1‑149 (MQ=255)
                                                                                                                          ccAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTg                          <  7:173677/149‑1 (MQ=255)
                                                                                                                           cAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGGCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTgg                         <  1:162075/149‑1 (MQ=255)
                                                                                                                              cACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGt                      >  4:315480/1‑149 (MQ=255)
                                                                                                                                 gggTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTa                                                                                                       >  4:63177/1‑65 (MQ=255)
                                                                                                                                 gggTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTa                                                                                                       <  3:63177/65‑1 (MQ=255)
                                                                                                                                 gggTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAg                   >  8:327792/1‑149 (MQ=255)
                                                                                                                                 gggTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAg                   >  3:336728/1‑149 (MQ=255)
                                                                                                                                   gTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTg                 >  8:186083/1‑149 (MQ=255)
                                                                                                                                   gTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTg                 >  4:114551/1‑149 (MQ=255)
                                                                                                                                                  gTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCg  >  3:950/1‑149 (MQ=255)
                                                                                                                                                  gTGCTTTGGCTTCGCTAATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCg  <  6:205525/149‑1 (MQ=255)
                                                                                                                                                  gTGCTTTGGCTTCGATCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCg  <  8:146693/149‑1 (MQ=255)
                                                                                                                                                   tGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCg  <  1:186291/148‑1 (MQ=255)
                                                                                                                                                   tGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCg  >  2:186291/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTGTCACAGAACCGAGACTATCAAATGAAGGATTTAAACGGGCTAAGCCTTCTGCGCTGGCGTCAGTGCGTGGCTGTTCATCG  >  NC_000913/2984654‑2984948

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