Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,154,782 A→G E224G (GAA→GGA)  yceG → septation protein, ampicillin sensitivity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,154,7820AG97.3% 124.5 / ‑3.8 37E224G (GAA→GGA) yceGseptation protein, ampicillin sensitivity
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (15/21);  total (16/21)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.35e-01

CTCAAGCGAGCGCACAAGAAAATGGTGAAAGCGGTCGATAGCGCCTGGGAAGGGCGTGCGGACGGTCTGCCTTATAAAGATAAAAACCAGTTGGTGACGATGGCATCAATTATCGAAAAAGAAACCGCCGTTGCCAGTGAACGCGATAAGGTTGCCTCAGTATTTATCAACCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGT  >  NC_000913/1154643‑1154922
                                                                                                                                           |                                                                                                                                            
cTCAAGCGAGCGCACAAGAAAATGGTGAAAGCGGTCGATAGCGCCTGGGAAGGGCGTGCGGACGGTCTGCCTTATAAAGATAAAAACCAGTTGGTGACGATGGCATCAATTATCGAAAAAGAAACCGCCGTTGCCAGTGGACGCGATAAgg                                                                                                                                   >  1:363537/1‑151 (MQ=255)
  cAAGCGAGCGCACAAGAAAATGGTGAAAGCGGTCGATAGCGCCTGGGAAGGGCGTGCGGACGGTCTGCCTTATAAAGATAAAAACCAGTTGGTGACGATGGCATCAATTATCGAAAAAGAAACCGCCGTTGCCAGTGGACGCGATAAGGtt                                                                                                                                 >  2:20630/1‑151 (MQ=255)
          cgcACAAGAAAATGGTGAAAGCGGTCGATAGCGCCTGGGAAGGGCGTGCGGACGGTCTGCCTTATAAAGATAAAAACCAGTTGGTGACGATGGCATCAATTATCGAAAAAGAAACCGCCGTTGCCAGTGGACGCGATCAGGTTGCCTCAGt                                                                                                                         >  2:471188/1‑151 (MQ=255)
             acaAAAAAATGGTGAAAGCGGTCGATAGCGCCGGGAAAGGGCGGGCGGACGGTCTGCCTTAAAAAGAAAAAACCCAGTTGGTGACGATGCCACCAATTATCAAAAAAGAAACCGCCGTTGCCAGTGGACGCGATAAGGTTGCCTCAGTa                                                                                                                        <  2:17805/149‑1 (MQ=255)
               aaGAAAATGGTGAAAGCGGTCGATAGCGCCTGGGAAGGGCGTGCGGACGGTCTGCCTTATAAAGATAAAAACCAGTTGGTGACGATGGCATCAATTATCGAAAAAGAAACCGCCGTTGCCAGTGGACGCGATAAGGTTGCCTCAGTATTTa                                                                                                                    <  1:209632/151‑1 (MQ=255)
               aaGAAAATGGTGAAAGCGGTCGATAGCGCCTGGGAAGGGCGTGCGGACGGTCTGCCTTAAAAAGATAAAAACCAGTTGGTGACGATGGCATCAATTATCGAAAAAGAAACCGCCGTTGCCAGTGGACGCGATAAGGTTGCCTCAGTATTTa                                                                                                                    <  1:44806/151‑1 (MQ=255)
                aGAAAATGGTGAAAGCGGTCGATAGCGCCTGGGAAGGGCGTGCGGACGGTCTGCCTTATAAAGATAAAAACCAGTTGTTGACGATGGCATCAATTATCGAAAAAGAAACCGCCGTTGCCAGTGGACGCGATAAGGTTGCCTCAGTATTTAt                                                                                                                   <  2:91966/151‑1 (MQ=255)
                  aaaaTGGTGAAAGCGGTCGATAGCGCTTGGGAAGGGCGTGCGGACGGTCTGCATTATAAAGATAAAAACCATTTGGTGACGATGGCATAAATTATCGAAAAAGAAACCGCCGTTGCCAGTGGACGCGATAAGGTTGCCTCAGTATTTATCa                                                                                                                 <  1:247020/151‑1 (MQ=255)
                  aaaaTGGTGAAAGCGGTCGATAGCGCCTGGGAAGGGCGTGCGGACGGTCTGCCTTATAAAGATAAAAACCAGTTGGTGACGATGGCATCAATTATCGAAAAAGAAACCGCCGTTGCCAGTGGACGCGATAAGGTTGCCTCAGTATTTATc                                                                                                                  >  2:202102/1‑150 (MQ=255)
                    aagggTGAAAGCGGTCGATAGCGCCTGGGAAGGGCGTGCGGACGGTCTGCCTTATAAAGATAAAAACCAGTTGGTGACGATGGCATCAATTATCGAAAAAGAAACCGCCGTTGCCAGTGGACGCGATAAGGTTGCCTCAGTATTTATCAAc                                                                                                               <  1:596017/148‑1 (MQ=255)
                         tGAAAGCGGTCGATAGCGCCTGGGAAGGGCGTGCGGACGGTCTGCCTTATAAAGATAAAAACCAGTTGGTGACGATGGCATCAATTATCGAAAAAGAAACCGCCGTTGCCAGTGGACGCGATAAGGTTGCCTCAGTATTTATCAACCGtt                                                                                                           <  1:73700/150‑1 (MQ=255)
                           aaaGCGGTCGATAGCGCCTGGGAAGGGCGTGCGGACGTTCTCCCTTAAAAAAATAAAACCCAGTTGGTGACAATGGCATCAATTATCGAAAAAGAAACCGCCGTTGCCAGTGGACGCGATAAGGTTGCCTCAGTATTTATCAACCGTTTAc                                                                                                        <  1:44270/151‑1 (MQ=255)
                            aaGCGGTCGATAGCGCCTGGGAAGGGCGTGCGGACGGTCTGCCTTATAAAGATAAAAACCAGTTGGTGACGATGGCATCAATTATCGAAAAAGAAACCGCCGTTGCCAGTGGACGCGATAAGGTTGCCTCAGTATTTATCAACCGTTTACg                                                                                                       >  1:327564/1‑151 (MQ=255)
                              gcgtgcaaTAGCGCAGGGGAAGGGCGTGCGGACGGTCTCCCTTATAAAGAAAAAAACCAGTTGGTGACGTTGGCATCAATTATCGAAAAAGAATCCGCCGTTGCCAGTGGACGCGATAAGGTTGCCTCAGTATTTATCAACCGTTTACGc                                                                                                      <  1:550004/143‑1 (MQ=255)
                                     aTAGCGCCTGGGAAGGGCGTGCGGACGGTCTGCCTTATAAAGATAAAAACCAGTTGGTGACGATGGCATCAATTATCGAAAAAGAAACCGCCGTTGCCAGTGGACGCGATAAGGTTGCCTCAGTATTTATCAACCGTTTACGCATTGGTAt                                                                                              >  1:554951/1‑151 (MQ=255)
                                         cgcCTGGGAAGGGCGGGCGGACGGTCTGCCTTAAAAAGAAAAAAACCAGTTGGTGACGATGGCATCATTTTTCGAAAAAGAAACCGCCGTTGCCAGTGGACGCGATAAGGTTGCCTCAGTATTTATCAACCGTTTACGCATTGGTATg                                                                                             <  2:9742/148‑1 (MQ=255)
                                           ccgggggAGGGCGTGCGGACGGTCTGCCTTAAAAAGATAAAACCCAGTTGGTGACGATGGCATAAATTATCGAAAAAGAAACCGCCGTTGCCAGTGGACGCGATAAGGTTGCCTCAGTATTTATCAACCGTTTACGCATTGGTATGCGCCt                                                                                        <  2:21282/148‑1 (MQ=255)
                                                   gggCGTGCGGACGGTCTGCCTTATAAAGATAAAAACCAGTTGGTGACGATGGCATCAATTATCGAAAAAGAAACCGCCGTTGCCAGTGGACGCGATAAGGTTGCCTCAGTATTTATCAACCGTTTACGCATTGGTATGCGCCTGCAgaccg                                                                                <  1:517357/151‑1 (MQ=255)
                                                               ggTCTGCCTTATAAAGATAAAAACCAGTTGGTGACGATGGCATCAATTATCGAAAAAGAAACCGCCGTTGCCAGTGGACGCGATAAGGTTGCCTCAGTATTTATCAACCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTg                                                                     >  2:265472/1‑150 (MQ=255)
                                                                      cTTATAAAGATAAAAACCAGTTGGTGACGATGGCATCAATTATCGAAAAAGAAACCGCCGTTGCCAGTGGACGCGATAAGGTTGCCTCAGTATTTATCAACCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACgg                                                             >  1:397340/1‑151 (MQ=255)
                                                                               aTAAAAACCAGTTGGTGACGATGGCATCAATTATCGAAAAAGAAACCGCCGTTGCCAGTGGACGCGATAAGGTTGCCTCAGTATTTATCAACCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATggg                                                       >  2:41365/1‑148 (MQ=255)
                                                                                 aaaaaCCAGTGGGTGACGATGGCATAAATTATCGAAAAAGAAACCGCCGTTGCCAGTGGACGCGATAAGGTTGCCTCAGTATTTATCAACCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGagag                                                   <  2:554940/150‑1 (MQ=255)
                                                                                      ccAGTTGGGGACGATGGCATAAATTATCAAAAAAGAAACCGCCGTTGCCAGTGGACGCGATAAGGTTGCCTCAGTATTTATCAACCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTta                                              <  2:230208/150‑1 (MQ=255)
                                                                                       cAGTTGGTGACGATGGCATCAATTATCGAAAAAGAAACCGCCGTTGCCAGTGGACGCGATAAGGTTGCCTCAGTATTTATCAACCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTtat                                             >  1:65407/1‑150 (MQ=255)
                                                                                         gTTGGTGACGATGGCATCAATTATCGAAAAAGAAACCGCCGTTGCCAGTGAACGCGATAAGGTTGCCTCAGTATTTATCAACCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAAt                                          >  1:99192/1‑151 (MQ=255)
                                                                                               gACGATGGCATCAATTATCGAAAAAGAAACCGCCGTTGCCAGTGGACGCGATAAGGTTGCCTCAGTATTTATCAACCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCaa                                     >  2:593775/1‑150 (MQ=255)
                                                                                                aCGATGGCATCAATTATCGAAAAAGAAACCGCCGTTGCCAGTGGACGCGATAAGGTTGCCTCAGTATTTATCAACCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCaa                                     <  1:433853/149‑1 (MQ=255)
                                                                                                aCGATGCCATCAATTATCGAAAAAAAAACCGCCGTTGCCAGTGGACGCGATAAGGTTGCCTCAGTATTTATCAACCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAAc                                   <  2:626483/151‑1 (MQ=255)
                                                                                                     ggCATAAATTATCGAAAAAGAAACCGCCGTTGCCAGTGGACGCGATAAGGTTGCCTCAGTATTTATCAACCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACttt                                <  2:594694/149‑1 (MQ=255)
                                                                                                        aTCAATTATCGAAAAAGAAACCGCCGTTGCCAGTGGACGCGATAAGGTTGCCTCAGTATTTATCAACCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGt                           >  1:41529/1‑151 (MQ=255)
                                                                                                           aaTTATCGAAAAAGAAACCGCCGTTGCCAGTGGACGCGATAAGGTTGCCTCAGTATTTATCAACCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGc                         >  1:112636/1‑150 (MQ=255)
                                                                                                                    aaaaGAAACCGCCGTTGCCAGTGGACGCGATAAGGTTGCCTCAGTATTTATCAACCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGaa               <  2:363530/151‑1 (MQ=255)
                                                                                                                     aaaGAAACCGCCGTTGCCAGTGGACGCGATAAGGTTGCCTCAGTATTTATCAACCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGaaa              <  2:193354/151‑1 (MQ=255)
                                                                                                                         aaaCCGCCGTTGCCAGTGGACGCGATAAGGTTGCCTCAGTATTTATCAACCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGcc          >  2:584311/1‑151 (MQ=255)
                                                                                                                         aaaCCGCCGTTGCCAGTGGACGCGATAAGGTTGCCTCAGTATTTATCAACCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGcc          <  1:44370/151‑1 (MQ=255)
                                                                                                                               ccgTTGCCAGTGGACGCGATAAGGTTGCCTCAGTATTTATCAACCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGCCAGc    <  1:384116/151‑1 (MQ=255)
                                                                                                                                 gTTGCCAGTGGACGCGATAAGGTTGCCTCAGTATTTATCAACCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGt  >  1:334019/1‑151 (MQ=255)
                                                                                                                                           |                                                                                                                                            
CTCAAGCGAGCGCACAAGAAAATGGTGAAAGCGGTCGATAGCGCCTGGGAAGGGCGTGCGGACGGTCTGCCTTATAAAGATAAAAACCAGTTGGTGACGATGGCATCAATTATCGAAAAAGAAACCGCCGTTGCCAGTGAACGCGATAAGGTTGCCTCAGTATTTATCAACCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGT  >  NC_000913/1154643‑1154922

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