Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,089,472 G→T 13.3% D199Y (GAC→TAC)  potF‑IV → putrescine‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,089,4720GT13.3% 77.2 / 4.1 30D199Y (GAC→TAC) potF‑IVputrescine‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (10/16);  new base T (2/2);  total (12/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.05e-01

TCGGCGACAACGCGCCGGTGGACAGCTGGGACCTGATCTTCAAGGAAGAGAACATCGCCAAGCTCAAGCAGTGCGGCGTGGCCCTGCTCGACTCGCCGTCGGAGATTCTACCGCTGGCGCTGCAGTACCTGGGCTTGCCGCCGAACAGCGACAAACCGGCGGACTACAAGAAGGCCGAGGAACTGATGCTGAAGATTCGTCCGCATATCACCTACTTCCACTCCTCGAAGTACATGGCAGACATCGCCAATGGCGACATCTGTGTGGCAGTGGGCTACAGCGGCAGCTTCTCGCAGGCA  >  NC_002947/6089323‑6089621
                                                                                                                                                     |                                                                                                                                                     
tCGGCGACAACGCGCCGGTGGACAGCTGGGACCTGATCTTCAAGGAAGAGAACATCGCCAAGCTCAAGCAGTGCGGCGTGGCCCTGCTCGACTCGCCGTCGGAGATTCTACCGCTGGCGCTGCAGTACCTGGGCTTGCCGCCGAACAGCg                                                                                                                                                       >  1:719663/1‑150 (MQ=255)
tCGGCGACAACGCGCCGGTGGACAGCTGGGACCTGATCTTCAAGGAAGAGAACATCGCCAAGCTCAAGCAGTGCGGCGTGGCCCTGCTCGACTCGCCGTCGGAGATTCTACCGCTGGCGCTGCAGTACCTGGGCTGGCCGCCGAACAGCg                                                                                                                                                       >  1:694956/1‑150 (MQ=255)
          cgcgCCGGTGGACAGCTGGGACCTGATCTTCAAGGAAGAGAACATCGCCAAGCTCAAGCAGTGCGGCGTGGCCCTGCTCGACTCGCCGTCGGAGATTCTACCGCTGGCGCTGCAGTACCTGGGCTTGCCGCCGAACAGCGACAAACcggc                                                                                                                                             <  2:21710/150‑1 (MQ=255)
             gcCGGTGGACAGCTGGGACCTGATCTTCAAGGAAGAGAACATCGCCAAGCTCAAGCAGTGCGGCGTGGCCCTGCTCGACTCGCCGTCGGAGATTCTACCGCTGGCGCTGCAGTACCTGGGCTTGCCGCCGAACAGCGACAAACCGGCGGa                                                                                                                                          <  1:304970/150‑1 (MQ=255)
                       aGCTGGGACCTGATCTTCAAGGAAGAGAACATCGCCAAGCTCAAGCAGTGCGGCGTGGCCCTGCTCGACTCGCCGTCGGAGATTCTACCGCTGGCGCTGCAGTACCTGGGCTTGCCGCCGAACAGCGACAAACCGGCGGACTACaagaag                                                                                                                                <  1:80300/150‑1 (MQ=255)
                                cTGATCTTCAAGGAAGAGAACATCGCCAAGCTCAAGCAGTGCGGCGTGGCCCTGCTCGACTCGCCGTCGGAGATTCTACCGCTGGCGCTGCAGTACCTGGGCTTGCCGCCGAACAGCGACAAACCGGCGGACTACAAGAAGGCCGAGGag                                                                                                                       >  1:356828/1‑149 (MQ=255)
                                    tCTTCAAGGAAGAGAACATCGCCAAGCTCAAGCAGTGCGGCGTGGCCCTGCTCGACTCGCCGTCGGAGATTCTACCGCTGGCGCTGCAGTACCTGGGCTTGCCGCCGAACAGCGACAAACCGGCGGACTACAAGAAGGCCGAGGAACTGa                                                                                                                   <  2:565268/150‑1 (MQ=255)
                                    tCTTCAAGGAAGAGAACATCGCCAAGCTCAAGCAGTGCGGCGTGGCCCTGCTCGACTCGCCGTCGGAGATTCTACCGCTGGCGCTGCAGTACCTGGGCTTGCCGCCGAACAGCGACAAACCGGCGGACTACAAGAAGGCCGAGGAACTGa                                                                                                                   <  2:356828/150‑1 (MQ=255)
                                    tCTTCAAGGAAGAGAACATCGCCAAGCTCAAGCAGTGCGGCGTGGCCCTGCTCGACTCGCCGTCGGAGATTCTACCGCTGGCGCTGCAGTACCTGGGCTTGCCGCCGAACAGCGACAAACCGGCGGACTACAAGAAGGCCGAGGAACTGa                                                                                                                   <  1:765133/150‑1 (MQ=255)
                                          aGGAAGAGAACATCGCCAAGCTCAAGCAGTGCGGCGTGGCCCTGCTCGACTCGCCGTCGGAGATTCTACCGCTGGCGCTGCAGTACCTGGGCTTGCCGCCGAACAGCGACAAACCGGCGGACTACAAGAAGGCCGAGGAACTGATGCTGa                                                                                                             <  2:311859/150‑1 (MQ=255)
                                                agaACATCGCCAAGCTCAAGCAGTGCGGCGTGGCCCTGCTCGACTCGCCGTCGGAGATTCTACCGCTGGCGCTGCAGTACCTGGGCTTGCCGCCGAACAGCGACAAACCGGCGGACTACAAGAAGGCCGAGGAACTGATGCTGAAGATTc                                                                                                       <  1:750377/150‑1 (MQ=255)
                                                       cGCCAAGCTCAAGCAGTGCGGCGTGGCCCTGCTCGACTCGCCGTCGGAGATTCTACCGCTGGCGCTGCAGTACCTGGGCTTGCCGCCGAACAGCGACAAACCGGCGGACTACAAGAAGGCCGAGGAACTGATGCTGAAGATTCGTCCGCa                                                                                                <  2:843222/150‑1 (MQ=255)
                                                       cGCCAAGCTCAAGCAGTGCGGCGTGGCCCTGCTCGACTCGCCGTCGGAGATTCTACCGCTGGCGCTGCAGTACCTGGGCTTGCCGCCGAACAGCGACAAACCGGCGGACTACAAGAAGGCCGAGGAACTGATGCTGAAGATTCGTCCGCa                                                                                                <  1:618206/150‑1 (MQ=255)
                                                          cAAGCTCAAGCAGTGCGGCGTGGCCCTGCTCGACTCGCCGTCGGAGATTCTACCGCTGGCGCTGCAGTACCTGGGCTTGCCGCCGAACAGCTACAAACCGGCGGACTACAAGAAGGCCGAGGAACTGATGCTGAAGATTCGTCCGCatat                                                                                             >  1:166219/1‑150 (MQ=255)
                                                          aaaGCTCAAGCAGTGCGGCGTGGCCCTGCTCGACTCGCCGTCGGAGATTCTACCGCTGGCGCTGCAGTACCTGGGCTTGCCGCCGAACAGCTACAAACCGGCGGACTACAAGAAGGCCGAGGAACTGATGCTGAAGATTCGTCCGCatat                                                                                             >  1:166224/2‑150 (MQ=255)
                                                           aaGCTCAAGCAGTGCGGCGTGGCCCTGCTCGACTCGCCGTCGGAGATTCTACCGCTGGCGCTGCAGTACCTGGGCTTGCCGCCGAACAGCTACAAACCGGCGGCCTACAAGAAGGCCGAGGAACTGATGCTGAAGATTCGTCCGCATATc                                                                                            <  2:166219/150‑1 (MQ=255)
                                                           aaGCTCAAGCAGTGCGGCGTGGCCCTGCTCGACTCGCCGTCGGAGATTCTACCGCTGGCGCTGCAGTACCTGGGCTTGCCGCCGAACAGCTACAAACCGGCGGAATACAAGAAGGCCGAGGAACTGATGCTGAAGATTCGTCCGCATATc                                                                                            <  2:166224/150‑1 (MQ=255)
                                                                                  ccTGCTCGACTCGCCGTCGGAGATTCTACCGCTGGCGCTGCAGTACCTGGGCTTGCCGCCGAACAGCGACAAACCGGCGGACTACAAGAAGGCCGAGGAACTGATGCTGAAGATTCGTCCGCATATCACCTACTTCCACTCCTCGAAGTa                                                                     >  2:126857/1‑150 (MQ=255)
                                                                                       tCGACTCGCCGTCGGAGATTCTACCGCTGGCGCTGCAGTACCTGGGCTTGCCGCCGAACAGCGACAAACCGGCGGACTACAAGAAGGCCGAGGAACTGATGCTGAAGATTCGTCCGCATATCACCTACTTCCACTCCTCGAAGTACATgg                                                                <  1:126857/150‑1 (MQ=255)
                                                                                        cGACTCGCCGTCGGAGATTCTACCTCTGGCGCTGCAGTACCTGGGCTTGCCGCCGAACAGCGACAAACCGGCGGACTACAAGAAGGCCGAGGAACTGATGCTGAAGATTCGTCCGCATATCACCTACTTCCACTCCTCGAAGTACATGGc                                                               <  2:749814/150‑1 (MQ=255)
                                                                                         gACTCGCCGTCGGAGATTCTACCGCTGGCGCTGCAGTACCTGGGCTTGCCGCCGAACAGCGACAAACCGGCGGACTACAAGAAGGCCGAGGAACTGATGCTGAAGATTCGTCCGCATATCACCTACTTCCACTCCTCGAAGTACATGGCa                                                              >  1:110171/1‑150 (MQ=255)
                                                                                         gACTCGCCGTCGGAGATTCTACCGCTGGCGCTGCAGTACCTGGGCTTGCCGCCGAACAGCGACAAACCGGCGGACTACAAGAAGGCCGAGGAACTGATGCTGAAGATTCGTCCGCATATCACCTACTTCCACTCCTCGAAGTACATGGCa                                                              <  2:719663/150‑1 (MQ=255)
                                                                                         gACTCGCCGTCGGAGATTCTACCGCTGGCGCTGCAGTACCTGGGCTTGCCGCCGAACAGCGACAAACCGGCGGACTACAAGAAGGCCGAGGAACTGATGCTGAAGATTCGTCCGCATATCACCTACTTCCACTCCTCGAAGTACATGGCa                                                              >  1:526559/1‑150 (MQ=255)
                                                                                                             aCCGCTGGCGCTGCAGTACCTGGGCTTGCCGCCGAACAGCGACAAACCGGCGGACTACAAGAAGGCCGAGGAACTGATGc                                                                                                                >  2:650162/1‑80 (MQ=255)
                                                                                                             aCCGCTGGCGCTGCAGTACCTGGGCTTGCCGCCGAACAGCGACAAACCGGCGGACTACAAGAAGGCCGAGGAACTGATGc                                                                                                                <  1:650162/80‑1 (MQ=255)
                                                                                                                    gcgcTGCAGTACCTGGGCTTGCCGCCGAACAGCGACAAACCGGCGGACTACAAGAAGGCCGAGGAACTGATGCTGAAGATTCGTCCGCATATCACCTACTTCCACTCCTCGAAGTACATGGCAGACATCGCCAATGGCGACATCtgtgtg                                   >  2:50128/1‑150 (MQ=255)
                                                                                                                      gcTGCAGTACCTGGGCTTGCCGCCGAACAGCGACAAACCGGCGGACTACAAGAAGGCCGAGGAACTGATGCTGAAGATTCGTCCGCATATCACCTACTTCCACTCCTCGAAGTACATGGCAGACATCGCCAATGGCGACATCTGTGTGGc                                 >  2:618601/1‑150 (MQ=255)
                                                                                                                      gcTGCAGTACCTGGGCTTGCCGCCGAACAGCGACAAACCGCCCGACGACAAGAAGGCCGAGGAACTGATGCTGAAGATTCGTCCGCATATCACCTACTTCCACTCCTCGAAGTACATGGCAGACATCGCCAATGGCGACATCTGTGTGGc                                 >  1:623306/1‑150 (MQ=255)
                                                                                                                          cAGTACCTGGGCTTGCCGCCGAACAGCGACAAACCGGCGGACTACAAGAAGGCCGAGGAACTGATGCTGAAGATTCGTCCGCATATCACCTACTTCCACTCCTCGAAGTACATGGCAGACATCGCCAATGGCGACATCTGTGTGGCAGTg                             <  2:624368/150‑1 (MQ=255)
                                                                                                                                cTGGGCTTGCCGCCGAACAGCGACAAACCGGCGGACTACAAGAAGGCCGAGGAACTGATGCTGAAGATTCGTCCGCATATCACCTACTTCCACTCCTCGAAGTACATGGCAGACATCGCCAATGGCGACATCTGTGTGGCAGTGGGCTAc                       >  1:635078/1‑150 (MQ=255)
                                                                                                                                                aCAGCGACAAACCGGCGGACTACAAGAAGGCCGAGGAGCTGAAGCTGAAGATACGCCCGCAAATCACCAACTTCAACTTCTCGAAGGACATGCAAGACATCGACAATGGCGTCATCTGCGGGGTATTGTGCTACAGCGGGAGATTCTCgc       >  2:358222/1‑150 (MQ=255)
                                                                                                                                                 cAGCGACAAACCGGCGGACTACAAGAAGGCCGAGGAACTGATGCTGAAGATTCGTCCGCATATCACCTACTTCCACTCCTCGAAGTACATGGCAGACATCGCCAATGGCGACATCTGTGTGGCAGTGGGCTACAGCGGCAGCTTCTCgca      <  1:618601/150‑1 (MQ=255)
                                                                                                                                                     gACAAACCGGCGGACTACAAGAAGGCCGAGGAGCTGATGCTGAAGATTCGTCCGCATATCACCTGCTTCCACTCCTCGAAGTACATGGCAGACATCGCCAATGGCGACATCTGTGTGGCAGTGGGCTACAGCGGCAGCTTCTCgcaggca  <  1:358222/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCGGCGACAACGCGCCGGTGGACAGCTGGGACCTGATCTTCAAGGAAGAGAACATCGCCAAGCTCAAGCAGTGCGGCGTGGCCCTGCTCGACTCGCCGTCGGAGATTCTACCGCTGGCGCTGCAGTACCTGGGCTTGCCGCCGAACAGCGACAAACCGGCGGACTACAAGAAGGCCGAGGAACTGATGCTGAAGATTCGTCCGCATATCACCTACTTCCACTCCTCGAAGTACATGGCAGACATCGCCAATGGCGACATCTGTGTGGCAGTGGGCTACAGCGGCAGCTTCTCGCAGGCA  >  NC_002947/6089323‑6089621

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