Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,991,432 C→T 23.6% P94S (CCA→TCA)  cyaA → adenylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,991,4320CT23.6% 19.8 / 3.6 17P94S (CCA→TCA) cyaAadenylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (3/10);  new base T (2/2);  total (5/12)
Fisher's exact test for biased strand distribution p-value = 5.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.22e-01

CTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGA  >  NC_000913/3991294‑3991574
                                                                                                                                          |                                                                                                                                              
ctgaTGCCGTGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAAt                                                                                                                                            <  2:249458/140‑1 (MQ=255)
   aTGCCGGGGTACCTTGATGGGAACGTACCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCGGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTAc                                                                                                                                         <  2:9995/143‑1 (MQ=255)
              ccTTGATGTTAACTTTTGCAAAGGTATTTGTTTTTTCATGCTTGATTAAACTCAACGTCATTTCCTGATTGAGCTTGAACTGTATCGTGGAATGTAAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATaca                                                                                                                              <  2:191331/143‑1 (MQ=255)
               cTTCATGGTAACGTTCCCAAAGGCATTTGCCTGTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTCTCGTGGAACGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATacac                                                                                                                             <  1:47963/143‑1 (MQ=255)
                    tGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGTAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATgg                                                                                                                        >  2:177531/1‑143 (MQ=255)
                     ggTAACGTTCCCAAATGCATTTGCCTTTACACGCCTGATGAAATTCAACGCCACTACCTTAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCagtaa                                                                                                                       <  1:334571/143‑5 (MQ=255)
                             tCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCt                                                                                                               <  1:207021/143‑1 (MQ=255)
                                 aaaGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGATTGAACGGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTc                                                                                                           >  1:363328/1‑143 (MQ=255)
                                                            gaAACTCAACGCCACTACCTAAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCTTCCGTAGGGCAAAGTTGTTCCTCTGACCt                                                                                <  2:119506/143‑1 (MQ=255)
                                                                   aaCGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTTCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGTGAAGTTGTCCCGCTTACATgaaaaaa                                                                         >  1:148066/1‑137 (MQ=255)
                                                                   aaCGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTTCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATc                                                                         >  1:64880/1‑143 (MQ=255)
                                                                    aCGTCACTGCCTGAACGAGCTTGGACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCt                                                                        <  2:240709/143‑1 (MQ=255)
                                                                    aCGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGACGAAAGGAGAGCTTCCAATTACTGGTGTATACACCATGCACAGCACCTCGTCCGTAAGGCAACGCTGTTCCTCTGACCTGGAAATCt                                                                        >  2:46595/1‑143 (MQ=255)
                                                                                               acaccccggACTGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTAGACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATa                                             <  1:162699/136‑1 (MQ=255)
                                                                                                                            cgAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCGCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTaaa                <  2:81874/143‑1 (MQ=255)
                                                                                                                                tggggAGCTTTGAATTACTGGTGTATACACCATTGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTa            <  2:148066/139‑1 (MQ=255)
                                                                                                                                aGGTGAGCTTTCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTa            <  2:64880/143‑1 (MQ=255)
                                                                                                                                          ccAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGa  <  1:246259/143‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                              
CTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTCAGTACAGGATCCGCCGAAAGGTGAGCTTCCAATTACTGGTGTATACACCATGGGCAGCACCTCGTCCGTAGGGCAAAGTTGTTCCTCTGACCTGGATATCTGGGTCTGTCATCAATCCTGGCTCGATAGCGAAGAGCGCCAATTGCTACAACGTAAATGTAGCCTGCTGGA  >  NC_000913/3991294‑3991574

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