Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,388,753 G→A G40D (GGT→GAT)  yfbO → uncharacterized protein YfbO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,388,7530GA96.8% 91.8 / ‑4.0 31G40D (GGT→GAT) yfbOuncharacterized protein YfbO
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (16/14);  total (16/15)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.23e-01

TGCTTAACATTCATGTTATGATGCTCTTGAAACTAATAATACGTATGGAATTGTCACTATTTATCTGCTAATGAATTCAAACAATGTCATTCGTTTGATGAATATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGATTTTTTTATTGACAATAATATCCATGGTTCAATATGTTGCAATGTGATTCCAGAACAATCTCCGCAAGCTATCTATCACCATTTTCTGAAAATACGTGAGCGTAACAATGTTAGTGATGTTTTGGTAGAAATTACGATGTTCGATGACCCAGACTGGCCCTTTAGTGAAAGTATTTTAGTTATCACCACTGCTTCACCAGAAGAAGTTCAGAGCTGGTTTGTAGAAGAAATTGCCCCTGATGAGTGCTGGGAAGGATGG  >  NC_000913/2388511‑2388985
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                        
tGCTTAACATTCATGTTATGATGCTCTTGAAACTAATAATACGTATGGAATTGTCACTATTTATCTGCTAATGAATTCAAACAATGTCATTCGTTTGATGAATATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGATTTTTTTATTGACAATAATATCCATGATTCAatat                                                                                                                                                                                                                                  >  2:72875/1‑251 (MQ=255)
    tAACATTCATGTTATGATGCTCTTGAAACTAATAATACGTATGGAATTGTCACTATTTATCTGCTAATGAATTCAAACAATGTCATTCGTTTGATGAATATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGATTTTTTTATTGACAATAATCTCCCTGATTCAGTATGt                                                                                                                                                                                                                                <  1:450051/249‑1 (MQ=255)
                      gCTCTTGAAACTAATAATACGTATGGAATTGTCACTATTTATCTGCTAATGAATTCAAACAATGTCATTCGTTTGATGAATATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGATTTTTTTATTGACAATAATATCCATGATTCAATATGTTGCAATGTGATTCCAGAAc                                                                                                                                                                                                             <  1:133284/250‑1 (MQ=255)
                      gCTCTTGAAACTAATAATACGTATGGAATTGTCACTATTTATCTGCTAATGAATTCAAACAATGTCATTCGTTTGATGAATATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGATTTTTTTATTGACAATAATATCCATGATTCAATATGTTGCAATGTGATTCCAGAAc                                                                                                                                                                                                             >  1:405823/1‑250 (MQ=255)
                       ctctTGAAACTAATAATACGTATGGAATTGTCACTATTTATCTGCTAATGAATTCAAACAATGTCATTCGTTTGATGAATATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGATTTTTTTATTGACAATAATATCCATGATTCAATATGTTGCAATGTGATTCCAGAACa                                                                                                                                                                                                            >  1:265095/1‑250 (MQ=255)
                              aaCTAATAATACGTATGGAATTGTCACTATTTATCTGCTAATGAATTCAAACAATGTCATTCGTTTGATGAATATATTATTTCAGGGAATAATTATGACACCTTTAGAACGACTTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGATTTTTTTATTGACAATAATATCCATGATTCAATATGTTGCAATGTGATTCCAGAACAATCTCCGc                                                                                                                                                                                                    >  1:290364/1‑251 (MQ=255)
                                                                                      tCATTCGTTTGATGAATATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAGTGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGATTTTTTTATTGACAATAGTATCCATGATTCAATATGTTGCAATGTGATTCCAGAACAATCTCCGCa                                                                                                                                                                                                   >  1:297874/1‑196 (MQ=255)
                                                                                      tCATTCGTTTGATGAAGATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATTGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGATTTTTTTATTGACAATAATAGCCATTATTCAATATGTTGCAATTTGATTTCAGAACAATCTCCGCa                                                                                                                                                                                                   <  2:297874/196‑1 (MQ=255)
                                                                                       cATTCGTTTGATGAATATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGATTTTTTTATTGACAATAATATCCATGATTCAATATGTTGCAATGTGATTCCAGAACAATCTCCGCAAGCTATCTATCACCATTTTCTGAAAATACGTGAGCGTAACAATGTTAGTGATGtt                                                                                                                                            <  1:2826/250‑1 (MQ=255)
                                                                                                      tatatTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGATTTTTTTATTGACAATACTATCCATGATTCAATATGTTGCAATGTGATTCCAGAACAATCTACGCAAGCTATCTATCACCATTTTCTGAAAATACGTGAGCGTAACAATGTTAGTGATGTTTTGGTAGAAATTACg                                                                                                                             >  2:370895/1‑250 (MQ=255)
                                                                                                            attTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGATTTTTTTATTGACAATAATATCCATGATTCAata                                                                                                                                                                                                                                   >  1:337676/1‑142 (MQ=255)
                                                                                                            attTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGATTTTTTTATTGACAATAATATCCATGATTCAata                                                                                                                                                                                                                                   <  2:337676/142‑1 (MQ=255)
                                                                                                            attTCAGGGAATAATTATGAAACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGATTTTTTTATTGACAATAATATCCATGATTCAATATGTTGCAATGTGATTCCAGAACAATCTCCGCAAGCTATCTATCACCATTTTCTGAAAATACGTGAGCGTAACAATGTTAGTGATGTTTTGGTAGAAATTACGATGTTCg                                                                                                                      <  2:8348/251‑1 (MQ=255)
                                                                                                                         atTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGATTTTTTTATTGACAATAATATCCATGATTCAATATGTTGCAATGTGATTCAAGAACAATCTCCGCAAGCTATCTATCACCATTTTCTGAAAATACGTGAGCGTAACAATGTTAGTGATGTTTTGGTAGAAATTACGATGTTCGATGACCAAGACTg                                                                                                         >  2:148511/1‑251 (MQ=255)
                                                                                                                              gACACCTTTTGAACGCATTACCCAACTTGTCAATATCAGTGGTGATGTTAATAATCCTGATAATCCTCGTACATTGTTATCTCTTGGGGATTTTTTTATTGACAATAATATCCATGATTCAATATTTTGCAATGTGATTCAAGAACAATCTCCGCAAGCTATCTATCACCATTTTCTGAAAATACGTGAGCGTCACAATGTTAGTGATGTTTTGTAAGAAATTACGATGTTCGTTGACCCAGACTGGCCCt                                                                                                    >  2:467076/1‑251 (MQ=255)
                                                                                                                                  ccTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGATTTTTTTATTGACAATAATATCCATGATTCAATATGTTGCAATGTGATTCCAGAACAATCTCCGCAAGCTATCTATCACCATTTTCTGAAAATACGTGAGCGTAACAATGTTAGTGATGTTTTGGTAGAAATTACGATGTTCGATGACCCAGACTGGCCCttt                                                                                                  >  2:465942/1‑249 (MQ=255)
                                                                                                                                            gCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGATTTTTTTATTGACAATAATATCCATGATTCAATATGTTGCAATGTGATTCCAGAACAATCTCCGCAAGCTATCTATCACCATTTTCTGAAAATACGTGAGCGTAACAATGTTAGTGATGTTTTGGTAGAAATTACGATGTTCGATGACCCAGACTGGCCCTTTAGTGAATGTAtt                                                                                      >  2:207879/1‑251 (MQ=255)
                                                                                                                                              aTTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGATTTTTTTATTGACAATAATATCCATGATTCAATATGTTGCAATGTGATTCCAGAACAATCTCCGCAAGCTATCTATCACCATTTTCTGAAAATACGTGAGCGTAACAATGTTAGTGATGTTTTGGTAGAAATTACGATGTTCGATGACCCAGACTGGCCCTTTAGTGAAAGTAtttt                                                                                    >  2:263940/1‑251 (MQ=255)
                                                                                                                                               ttACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGATTTTTTTATTGACAATAATATCCATGATTCAATATGTTGCAATGTGATTCCAGAACAATCTCCGCAAGCTATCTATCACCATTTTCTGAAAATACGTGAGCGTAACAATGTTAGTGATGTTTTGGTAGAAATTACGATGTTCGATGACCCAGACTGGCCCTTTAGTTAAAGTAttt                                                                                     <  1:136701/249‑1 (MQ=255)
                                                                                                                                               ttACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATATCTTGAGGATTTTTTTATTGACAATAATATCCATGATTCAATATGTTGCAATGTGATTCCAGAACAATCTCCGCAAGCTATCTATCACCATTTTATGAAAATACGTGAGCGTAACAATGTTAGTGATGTTTTGGTAGAAATTACGATGTTCGATGACCAAGACTGGCCCTTTAGTGAAAGTAttg                                                                                     >  2:136701/1‑248 (MQ=255)
                                                                                                                                                    cAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGATTTTTTTATTGACAATAATATCCATGATTCAATATGTTGCAATGTGATTCCAGAACAATCTCCGCAAGCTATCTAGCACCATTTTCTGAAAATACGTGAGCGTAACAATGTTAGTGATGTTTTGGTAGAAATTACGATGTTCGATGACCCAGACTGGCCCTTTAGTGAAAGTATTTTAGTTAt                                                                              <  1:453614/251‑1 (MQ=255)
                                                                                                                                                          gTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGATTTTTTTATTGACAATAATATCCATGGTTCAATATGTTGCAATGTGATTCCAGAACAATCTCCGCAAGCTATCTATCACCATTTTCTGAAAATACGTGAGCGTAACAATGTTAGTGATGTTTTGGTAGAAATTACGATGTTCGATGACCCAGACTGGCCCTTTAGTGAAAGTATTTTAGTTATcacca                                                                         <  1:260567/250‑1 (MQ=255)
                                                                                                                                                                                        gATACTCCTCGTCCATTGTTATCTCTTGAGGATTTTTTTATTGACAATAATATCCATGATTCAATATGTTGCAATGTGATTCCAGAACAATCTCCGCAAGCTATCTATCACCATTTTCTGAAAATACGTGAGCGTAACAATGTTAGTGATGTTTTGGTAGAAATTACGATGTTCGATGACCCAGACTGGCCCTTTAGTGAAAGTATTTTAGTTATCACCACTGCTTCACCAGAAGAAGTTCAGAGCTGGtt                                          <  1:72875/251‑1 (MQ=255)
                                                                                                                                                                                                  gTCCATTGTTATCTCTTGAGGATTTTTTTATTGACAATAATATCCATGATTCAATATGTTGCAATGTGATTCCAGAACAATCTCCGCAAGCTATCTATCACCATTTTCTGAAAATACGTGAGCGTAACAATGTTAGTGATGTTTTGGTAGAAATTACGATGTTCGATGACCCAGACTGGCCCTTTAGTGAAAGTATTTTAGTTATCACCACTGCTTCACCAGAAGAAGTTCAGAGCTGGTTTGTAgaagaa                                <  1:360897/251‑1 (MQ=255)
                                                                                                                                                                                                    ccATTGTTATCTCTTGAGGATTTTTTTATTGACAATAATATCCATGATTCAATATGTTGCAATGTGATTCCAGAACAATCTCCGCAAGCTATCTATCACCATTTTCTGAAAATACGTGAGCGTAACAATGTTAGTGATGTTTTGGTAGAAATTACGATGTTCGATGACCCAGACTGGCCCTTTAGTGAAAGTATTTTAGTTATCACCACTGCTTCACCAGAAGAAGTTCAGAGCTGGTTTGTAGAAGAAAt                              >  1:210279/1‑251 (MQ=255)
                                                                                                                                                                                                    ccATTGTTATCTCTTGAGGATTTTTTTATTGACAATAATATCCATGATTCAATATGTTGCAATGTGATTCCAGAACAATCTCCGCAAGCTATCTATCACCATTTTCTGAAAATACGTGAGCGTAACAATGTTAGTGATGTTTTGGTAGAAATTACGATGTTCGATGACCCAGACTGGCCCTTTAGTGAAAGTATTTTAGTTATCACCACTGCTTCACCAGAAGAAGTTCAGAGCTGGTTTGTAGAAGAAAt                              <  1:263940/251‑1 (MQ=255)
                                                                                                                                                                                                      aTTGTTATCTCTTGAGGATTTTTTTATTGACAATAATATCCATGATTCAATATGTTGCAATGTGATTCCAGAACAATCTCCGCAAGCTATCTATCACCATTTTCTGAAAATACGTGAGCGTAACAATGTTAGTGATGTTTTGGTAGAAATTACGATGTTCGATGACCCAGACTGGCCCTTTAGTGAAAGTATTTTAGTTATCACCACTGCTTCACCAGAAGAAGTTCAGAGCTGGTTTGTAGAAGAAATTg                            <  1:465942/251‑1 (MQ=255)
                                                                                                                                                                                                         gTTATCTCGTGAGGATTTTTTTATTGACAATAATATCCATGATTCAATATGTTGCAATGTGATTCCAGAACAATCTCCGCAAGCTATCTATCACCATTTTCTGAAAATACGTGAGCGTAACAATGTTAGTGATGTTTTGGTAGAAATTACGATGTTCGATGACCCAGACTGGCCCTTTAGTGAAAGTATTTTAGTTATCACCACTGCTTCACCAGAAGAAGTTCAGAGCTGGTTTGTAGAAGAAATTGccc                         <  2:265095/251‑1 (MQ=255)
                                                                                                                                                                                                                    aGGATTTTTTTATTGACAATAATATCCATGATTCAATATGTTGCAATGTGATTCCAGAACAATCTCCGCAAGCTATCTATCACCATTTTCTGAAAATACGTGAGCGTAACAATGTTAGTGATGTTTTGGTAGAAATTACGATGTTCGATGACCCAGACTGGCCCTTTAGTGAAAGTATTTTAGTTATCACCACTGCTTCACCAGAAGAAGTTCAGAGCTGGTTTGTAGAAGAAATTGCCCCTGATGAGTGc              <  1:131487/251‑1 (MQ=255)
                                                                                                                                                                                                                         ttttttATTGACAATAATATCCATGATTCAATATGTTGCAATGTGATTCCAGAACAATCTCCGCAAGCTATCTATCACCATTTTCTGAAAATACGTGAGCGTAACAATGTTAGTGATGTTTTGGTAGAAATTACGATGTTCGATGACCCAGACTGGCCCTTTAGTGAAAGTATTTTAGTTATCACCACTGCTTCACCAGAAGAAGTTCAGAGCTGGTTTGTAGAAGAAATTGCCCCTGATGAGTGCTggg          >  1:91558/1‑250 (MQ=255)
                                                                                                                                                                                                                                ttGACAATAATATCCATGATTCAATATGTTGCAATGTGATTCCAGAACAATCTCCGCAAGCTATCTATCACCATTTTCTGAAAATACGTGAGCGTAACAATGTTAGTGATGTTTTGGTAGAAATTACGATGTTCGATGACCCAGACTGGCCCTTTAGTGAAAGTATTTTAGTTATCACCACTGCTTCACCAGAAGAAGTTCAGAGCTGGTTTGTAGAAGAAATTGCCCCTGATGAGTGCTGGGAAGGATgg  >  2:84176/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                        
TGCTTAACATTCATGTTATGATGCTCTTGAAACTAATAATACGTATGGAATTGTCACTATTTATCTGCTAATGAATTCAAACAATGTCATTCGTTTGATGAATATATTATTTCAGGGAATAATTATGACACCTTTAGAACGCATTACCCAACTTGTTAATATCAATGGTGATGTTAATAATCCTGATACTCCTCGTCCATTGTTATCTCTTGAGGATTTTTTTATTGACAATAATATCCATGGTTCAATATGTTGCAATGTGATTCCAGAACAATCTCCGCAAGCTATCTATCACCATTTTCTGAAAATACGTGAGCGTAACAATGTTAGTGATGTTTTGGTAGAAATTACGATGTTCGATGACCCAGACTGGCCCTTTAGTGAAAGTATTTTAGTTATCACCACTGCTTCACCAGAAGAAGTTCAGAGCTGGTTTGTAGAAGAAATTGCCCCTGATGAGTGCTGGGAAGGATGG  >  NC_000913/2388511‑2388985

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