Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,628,984 A→G D513D (GAT→GAC ettA ← energy‑dependent translational throttle protein EttA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,628,9840AG97.0% 96.2 / ‑4.2 33D513D (GAT→GACettAenergy‑dependent translational throttle protein EttA
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (15/17);  total (15/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.53e-01

GTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCTGGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGCGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCGGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAATCCAGAATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAGAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGT  >  NC_000913/4628734‑4629228
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    
gTGCTCAATGCCTGATGCGACGTTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCTGGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGCGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCGGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTaa                                                                                                                                                                                                                                                      <  1:201310/251‑1 (MQ=255)
        tGCCTGATGCGACTCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGGTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAggat                                                                                                                                                                                                                                              >  1:115034/1‑248 (MQ=255)
                     gCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATa                                                                                                                                                                                                                                 <  1:376689/251‑1 (MQ=255)
                                  tatcagtcatgtatattcaCAACGTATTAAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCACCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAAc                                                                                                                                                                                                                     <  1:374077/232‑1 (MQ=255)
                                            ggagattcccAACGTATTGAACTTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCg                                                                                                                                                                                                            <  2:83770/240‑1 (MQ=255)
                                              agattcaCAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCg                                                                                                                                                                                                        <  1:116342/245‑1 (MQ=255)
                                                   caCAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGAt                                                                                                                                                                                                    <  1:441187/249‑1 (MQ=255)
                                                     cAACGTATTGAAATTGCTGGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGCGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCGGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAATCCAGAATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAAc                                                                                                                                                                                                 <  2:135254/251‑1 (MQ=255)
                                                       aCGTACTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCa                                                                                                                                                                                               <  1:94078/251‑1 (MQ=255)
                                                                       tggtttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGa                                                                                                                                                                                <  2:387235/247‑1 (MQ=255)
                                                                         gtttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACt                                                                                                                                                                             >  1:168984/2‑251 (MQ=255)
                                                                                   cGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGc                                                                                                                                                                   >  1:225078/1‑251 (MQ=255)
                                                                                          ggCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTcc                                                                                                                                                            <  1:352352/251‑1 (MQ=255)
                                                                                          ggCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTc                                                                                                                                                             >  2:490259/1‑250 (MQ=255)
                                                                                           gCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCa                                                                                                                                                           <  1:339120/251‑1 (MQ=255)
                                                                                                             gCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACAAACGGTCGTGCGAGATAACCATCGCACCGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCGGGGTTTCg                                                                                                                                          >  1:67641/1‑250 (MQ=255)
                                                                                                                      gCATTACTTCGCATTACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCTAGCGTGCGCTTCTTGTACTCTTCGTACTCACTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGTGGAACCACCGGCCGTGCGAGATAACCATCGCACAGCCCGGCAACTCCCGCAGGGCGTTTTCCAgcgcccacggggttccgtcatgcaggc                                                                                                                                >  2:82129/1‑228 (MQ=255)
                                                                                                                                    tACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTtcg                                                                                                                  >  1:428842/1‑251 (MQ=255)
                                                                                                                                           gTACTTGATATGCGTCGGCTCCATTGCGTCTGCGCCCAGCGTGCGGTTCTTGTACTCTTCGTACTCATTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGTTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGTGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCATGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGagca                                                                                                           <  1:323252/251‑1 (MQ=255)
                                                                                                                                                  aTACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTc                                                                                                                                                                                                             >  2:69926/1‑146 (MQ=255)
                                                                                                                                                  aTACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCAt                                                                                                    >  1:461825/1‑251 (MQ=255)
                                                                                                                                                  aTACGCTTCGGCGCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTGACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTc                                                                                                                                                                                                             <  1:69926/146‑1 (MQ=255)
                                                                                                                                                      gCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTg                                                                                                >  1:187504/1‑251 (MQ=255)
                                                                                                                                                          cGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACGCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCAGCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTgccg                                                                                             <  2:34424/250‑1 (MQ=255)
                                                                                                                                                               ccagtccgTCTGCGCCCAGCCTGCGTTTCTTGTACTCTTCGTACTCCGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCCTCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGCGCAGCAGCATGTTGCCGCCAAcc                                                                                       >  2:81248/7‑251 (MQ=255)
                                                                                                                                                               ccAGTGCGGCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGTTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGCTATCCAGGGCGTTGGTTGGTTCGTCGCGCAGCAGCATGTTTCCGCCCAcc                                                                                       >  2:408763/1‑251 (MQ=255)
                                                                                                                                                                             cccAGCGTGCGTTTCTTGTACGCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTTCCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGc                                                                         <  2:47133/251‑1 (MQ=255)
                                                                                                                                                                                                                       ccTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAAc                               >  1:303339/1‑251 (MQ=255)
                                                                                                                                                                                                                             aagaaCTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGttt                         >  1:336928/1‑251 (MQ=255)
                                                                                                                                                                                                                                    cAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTg                   <  1:259889/250‑1 (MQ=255)
                                                                                                                                                                                                                                       cgttACCTTCATCCTGGTAGGCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCa               <  2:1834/249‑1 (MQ=255)
                                                                                                                                                                                                                                          tACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCACCCTGCAGCAGCTTCGCCAGATGCAGACGCCCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAAcc            >  2:139830/1‑251 (MQ=255)
                                                                                                                                                                                                                                           aCCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAAcc            <  2:187504/250‑1 (MQ=255)
                                                                                                                                                                                                                                                      ggTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCACCCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTttaaagt  <  1:42454/249‑1 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    
GTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCTGGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGCGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCGGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAATCCAGAATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAGAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGT  >  NC_000913/4628734‑4629228

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