Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,630,385 G→A S46S (TCC→TCT)  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,630,3850GA93.9% 87.7 / ‑3.9 34S46S (TCC→TCT) ettAenergy‑dependent translational throttle protein EttA
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (16/15);  total (17/17)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATAATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGCCGGTCGATCACCCTCGCAGTAAA  >  NC_000913/4630159‑4630625
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                
ccaccTTGTCAAACCCGGCACCCGGACCGGCGTACGGCGCATCCCCTCCACCCGGGCGTTCCAGGGGTTTACCCACTCCAGAACCCGCTTCTTCAAGGGCCTCACGCAGGGTGTGTTCCGGGTTCAGCTGCGTTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGCCGCGCTTCCCCTTCGATGTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGAccc                                                                                                                                                                                                                          <  1:531106/247‑1 (MQ=255)
       gTCAAAATCGGCCTCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCATTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAAACAATCTTGATGTCTGGCTGCGGCCGCGCCTCACCTTCGATGTCTTTCTCAATCCCC‑CCCATAATGCGACGCAGGGTAGAATTGCCCGCGCCGTTCAGAcaccgggccc                                                                                                                                                                                                                   >  2:202344/1‑242 (MQ=255)
                       ggtgcGGCGTACAGCGCATACCCTTCACCCAGGCTTTTCAGGGCTTTAACCACTTCAAAACCCGTTTCTTCAATGGACTCCCGCACGGTGTGTTCCGGGTTCAGCGGCGTTTCCTGCGGCAGTTAACCAATCTTTATGTCTGGCGGCGGCCCCGTTTCCCCTTCGATTTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAgg                                                                                                                                                                                                   <  1:202344/247‑1 (MQ=255)
                          ccGGCGTACAGCGCAACCCCTTCATCCAGGCGTTCCAGGGCGTTACCCCCTCCAAAACCCTCTTCTTCAATGGACCCCCGCAGGGTGTGTTCCGGGTTCAGCTGCGTTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCCCCTTCGATGTCTTTATCAAGGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGgaa                                                                                                                                                                                                <  1:618337/250‑1 (MQ=255)
                               gTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGAAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAAc                                                                                                                                                                                           >  1:55247/1‑251 (MQ=255)
                                tACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCATTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGCAACCCATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGACGTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGTCGACATTCCCCCACCCATCAAACCCAGGCCCCCAAATTTTACCCCAggggaaaaacc                                                                                                                                                                                          >  2:363447/1‑243 (MQ=255)
                                  cAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGGAACCAATCTTGATGTCTGGCTGCGGAAGCGCTTCACCTTCGATGTCTTTATTAATGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCCCCGTTCAGACCCAGGACACCCATTTTTTCCCCCGGGGAGAAACTcc                                                                                                                                                                                        >  2:329570/1‑250 (MQ=255)
                                        aTACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCCATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGGCTTTTTCAATGCCC‑GCCATAATGCGCACCAGGGTAGGCTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCAAGGGAAGAAACTCcgag                                                                                                                                                                                     >  1:275853/1‑248 (MQ=255)
                                        aTACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGGAAACAATCTTGATGTCTGGCTGCGGACGCGCCTCACCTTCGGTGGCTTTATCAATGCCC‑GCCATAATGCGCCGCAGGGTAGACTTACCCCCGCCGTTCAGGCCCAGGACACCAATTTTTGTCCCAGGGAAGaaaatccgagggg                                                                                                                                                                                  >  2:342870/1‑241 (MQ=255)
                                             cTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGttt                                                                                                                                                                             >  1:141965/1‑251 (MQ=255)
                                               ccATCCAGGCGTTTAAGGGCGTTAACCCCTACAGAACCCGCTTCTTCATGGGCCTCACGCGCGGTGTGTTCGGGTTTCGGCTGCGGTTCCTGCGGCAGATAACCATTCTTGAGGTCTGGCTGCGGACGCGCTTCACCTTGGATGTTTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGtttt                                                                                                                                                                            <  1:328420/249‑1 (MQ=255)
                                               ccATCCAGGCGTTACAGGGCGTAAACCGCTTAAGAAACCTCTTCTTCAGGGGCCTCACGCACGGTGTGTCGCGGGTTCAGCTGGGGTTCCTGGGGCAGGTAACCAATCTGTATGTCGGGCTGCGGCCGCGCTTCACCTTCGTTGTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGttt                                                                                                                                                                             <  1:367624/248‑1 (MQ=255)
                                                           ttcggggCGTAAACCACTTAAGAACCCGTTTCTTCAATGGACTCCCGCCGGGTTTGTTCCGGTTTCAGCTGCGTTCCCGGCGGCAGGTAACAAATTTTGGTGTCGGGCTGCGGACGCGCTCCCCCTCCGATGTCTTTATCAAGGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGAc                                                                                                                                                                <  1:631452/246‑1 (MQ=255)
                                                                             taaaaaaCCGCTTCTTCAATGGACTCCCCCCGGGTGTGTCCCGGGTTCAGCGGCGGTTCCTGCGGCAGTTACCCAATTTTGATGTCTGGCGGCGGCCGCCCTTCCCCTTCGATGTCTTAATCAAGGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACtt                                                                                                                                              <  1:613994/246‑1 (MQ=255)
                                                                                  aaCCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCg                                                                                                                                         >  1:598542/1‑250 (MQ=255)
                                                                                                     ccccgcgcggggtgtTCGGGGTTCCGCTGGGGTTCCGGCGGCAGGTACCCATTCTTGATGTCTGGCTGCGGACGCGCTTCCCCTGCGATGTCTTTA‑CAAAGCCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCACGGTATAAAc                                                                                                                     <  1:605224/240‑1 (MQ=255)
                                                                                                                ggtttCCGGTTTCATCTGCGTTCCCTGCGGCAGGTACCCAATCTTGTTGTCTGGCGGCGGCCCCTCTTCCCCTTCAATGTCTTTACCAATGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCa                                                                                                          <  1:656151/248‑1 (MQ=255)
                                                                                                                        gtttCAGCGGCGGTCCCGGGGGCAGTTACCCAATCTTGATGTCGGGCTGCGGCCGCTCTTCCCCTTCGATGTTTTTATAAAGCCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGAc                                                                                                  <  1:511637/249‑1 (MQ=255)
                                                                                                                          ttCAGCTGCGGTTTCTGCGGCAGGTAACCAATCTTGCTGTCTGGCTGAGGTCGCGCCTCACCTTCGATGTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCCCGCCCTTCAGACCCAAGACACCAATTTTTGCCCCAGGGCAGAAACTCAAAGAGATGTTTTTCAAAATATGACCTTTCGGCGGAACAACTTTGCCGACACGGTGAGTGATATTAACGCATTACGCCACGTTGGACtt                                                                                                >  2:360874/1‑251 (MQ=255)
                                                                                                                                             gCAGGTAACCAATCTTGATGTATGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGGGCCACGTTGGACTTCGCCTCTTTGTTTATCGTg                                                                             >  2:85201/1‑251 (MQ=255)
                                                                                                                                                 aTAACCAATCTTGTTGTCTGGCTGCGGACGCGCTTCCCCTTCGATGTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATaa                                                                         <  1:113674/251‑1 (MQ=255)
                                                                                                                                                          cTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTc                                                                >  1:416388/1‑251 (MQ=255)
                                                                                                                                                          cTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACCTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTc                                                                >  1:192713/1‑251 (MQ=255)
                                                                                                                                                            tGATGTCTGGCGGCGGACGCGCTTCCCCTTCGAGGTTTTTACAAATGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTc                                                                <  1:85201/249‑1 (MQ=255)
                                                                                                                                                                   gggtcgggagcccgcTTCCCCTTGGATTTCTTTACCAAGCCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGaa                                                       <  1:363447/239‑1 (MQ=255)
                                                                                                                                                                                    cccctgcGATGTCTTTCTCAAGGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGTAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTcccc                                      <  1:179906/245‑1 (MQ=255)
                                                                                                                                                                                           gATGTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGAAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGGCACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGCTGATGAGTTTTCAAAGGCGAAGTGTAGCATTTTTCCCCGCCTAAt                               >  2:496349/1‑251 (MQ=255)
                                                                                                                                                                                             tttCTTTATCAACCCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCTTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATc                              <  1:30381/248‑1 (MQ=255)
                                                                                                                                                                                             gttgtaaactaagcccc‑GCCATAATGTGCAGCAGGTGGGTCGACCCCGCGCCGTTCAGACCCAGGACCCCAATTTGTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGAGGCTTGGTATAAACAAATTGGGCCCCGTTGGACTTCGCCTCTATTTTTATCGTGATGATGAGTTGTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATc                              <  2:463205/236‑1 (MQ=255)
                                                                                                                                                                                                 tttATCAATGCCC‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAgt                         >  1:640634/1‑250 (MQ=255)
                                                                                                                                                                                                            cc‑GCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCAc              >  1:60789/1‑251 (MQ=255)
                                                                                                                                                                                                               gCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGCCCCATGACACCAATTTTTGCCCCAGTGACGAAACTCAGAGAGATGTTTTTCAACATATGACGTTTCGGCGGAACAACTTTGCCGACCCGATGCATGGTCTAAACGATTTGAGCCTCGTTGGCCTTCGCCTCAATTTTTATCGTGATGCTGAGTTTTCAAAGGCGAAGTGTCTCCTTTTTTCCCGCCTAATCCCAGTTGGTCGAT‑ACCCt           >  2:294110/1‑251 (MQ=255)
                                                                                                                                                                                                                         gcagcaGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCagtaaa  <  1:329570/251‑1 (MQ=255)
                                                                                                                                                                                                                         gcagcaGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCagtaa   >  1:589543/1‑250 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                
CCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCC‑GCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATAATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGCCGGTCGATCACCCTCGCAGTAAA  >  NC_000913/4630159‑4630625

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