Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,546,084	G→T	24.0%	H18Q (CAC→CAA)	yjiC ←	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,546,084	0	G	T	24.0%	37.5 / 9.8	25	H18Q (CAC→CAA)	yjiC	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (5/14);  new base T (2/4);  total (7/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 = 9.04e-01

TTTCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAATGA  >  CP009273/4545948‑4546221                                                                                                                                         |                                                                                                                                          tttCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGta                                                                                                                               >  1:284787/1‑149 (MQ=255)        aaCGCTTTGTCCACCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAAttt                                                                                                                        <  1:230026/149‑1 (MQ=255)                    aTCGTAATTATGGGTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAggg                                                                                                            <  1:26670/149‑1 (MQ=255)                       gTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGaa                                                                                                               <  2:165579/143‑1 (MQ=255)                            tgtGGTTATTGACACGAGTAAGAGGCGTAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCa                                                                                                          <  2:127459/141‑1 (MQ=255)                            tATGGTTAGTGCCACGAGTCACAGGCGGACGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTTAAGGGCa                                                                                                          <  2:48888/143‑1 (MQ=255)                                    ttGCCACGAAGAAGAGGCGGAAGATAGATGTTATCCGAATCCATGGTTGCATATTTTAATTTACGTGTAAACTCACGGTGGACTATTGGGTGTTGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTaa                                                                                                  <  2:116931/143‑1 (MQ=255)                                        cACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGGCTCTGAAGGGCATTGCTTAATGGc                                                                                              >  2:361805/1‑143 (MQ=255)                                             gTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACataat                                                                                   >  1:288259/1‑149 (MQ=255)                                              tAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACataata                                                                                  <  1:24667/149‑1 (MQ=255)                                                 gagGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTCAAGGGCATTGCTTAATGGCATAGACATAATAGtt                                                                               >  1:449960/1‑149 (MQ=255)                                                     cGGAATATGGATGGTTTCCGGATCCATGCGTTCAGATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCt                                                                           <  1:182923/149‑1 (MQ=255)                                                     cGGAAGATAGTTGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGCGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAg                                                                                 <  2:451652/143‑1 (MQ=255)                                                            tagatGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGa                                                                    >  1:331946/1‑149 (MQ=255)                                                                                    gCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTa                                                  <  2:421954/143‑1 (MQ=255)                                                                                    gCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGACGAAGGAATTGAGCATCTGTAATAATTTGACTCAGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTGTTTCACTGATGTTTa                                                  <  2:296305/143‑1 (MQ=255)                                                                                          tttAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGaaaaa                                            <  2:253273/143‑1 (MQ=255)                                                                                            tAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATaaa                                    <  1:315256/149‑1 (MQ=255)                                                                                               tttACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTACGAGTATATCACTGATGTTTAGAAAAATAGATAAAttt                                 <  1:128237/149‑1 (MQ=255)                                                                                                    gtgtAAACACACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTg                            <  1:403600/149‑1 (MQ=255)                                                                                                                tGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAagag                <  1:424141/149‑1 (MQ=255)                                                                                                                    gACTATTGGGTGATGAAGGAATTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGt            >  1:1243/1‑149 (MQ=255)                                                                                                                      ttgtggggTGATGAAGGAAGTGAGCATCTGTAATAATTTGGCTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAagag                <  2:35307/138‑1 (MQ=255)                                                                                                                            ggTGATGAAGGAAGTGGGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAAt    <  1:193066/149‑1 (MQ=255)                                                                                                                            ggTGAGGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTaa          >  2:290642/1‑143 (MQ=255)                                                                                                                                    agtaagTGAGCAGCTGAAATAAGTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAAGAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAATGa  <  2:187278/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                          TTTCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAATGA  >  CP009273/4545948‑4546221

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