Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	2,420,353	A→G	57.8%	intergenic (‑86/+135)	hisJ ← / ← argT	histidine/lysine/arginine/ornithine transporter subunit/lysine/arginine/ornithine transporter subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	2,420,353	0	A	G	57.8%	10.0 / 27.3	29	intergenic (‑86/+135)	hisJ/argT	histidine/lysine/arginine/ornithine transporter subunit/lysine/arginine/ornithine transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (8/4);  new base G (10/6);  total (18/11)
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.94e-01

GCAGTTGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGACGTATCAGT  >  CP009273/2420212‑2420492                                                                                                                                              |                                                                                                                                            gcAGTTGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCAACAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGGGCCATAAAATTGACCATCAAGGCAAGATGTcca                                                                                                                                      >  1:114071/1‑148 (MQ=255)                aaGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCACGGCACAGTATCCCGCCCCGTTAGACGta                                                                                                                       >  1:305070/1‑149 (MQ=255)                     cAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATCTTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTattatt                                                                                                                  >  1:215745/1‑149 (MQ=255)                              agagagCGATAGCACCTGTTTTTTCATAATGTCCTCAACTGACAGTTTTTTTAGAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAATATATCCCGCCCCGTTAGACGTATTATTTTTc                                                                                                              <  2:150280/144‑1 (MQ=255)                                            ccAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACATGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTa                                                                                                <  2:144293/144‑1 (MQ=255)                                               gTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGGATATTTACAGGTTTGTTACACTTCTCTGTCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCa                                                                                        >  1:272165/1‑149 (MQ=255)                                                    ttCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTCTGCa                                                                                        >  2:40193/1‑144 (MQ=255)                                                              ccTCAACTGACAGTCTTTTTATAGGGAACATTTACAGTTCTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCACGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGca                                                                              <  2:17488/144‑1 (MQ=255)                                                               cTCAACTGACAGTCTTTTTATAGGGAATATTTACACGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGa                                                                        >  1:376630/1‑149 (MQ=255)                                                                tCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGAt                                                                       >  1:368271/1‑149 (MQ=255)                                                                tCAACTGACAGTCTTTTTATAGGGAATATTTACAGGGTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGAt                                                                       >  1:226804/1‑149 (MQ=255)                                                                 tatctgtctgtttttttATAGGATATATTTACAGGGTTGTTACACTTATCGGGCCATAAAATTCACCATCCAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATa                                                                      <  1:33482/137‑1 (MQ=255)                                                                      gACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGACATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGca                                                                 <  1:112717/149‑1 (MQ=255)                                                                         aGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGcaaca                                                              >  1:33508/1‑149 (MQ=255)                                                                                 tataGGGAATATTTACAGGATTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAAc                                                           >  2:194311/1‑144 (MQ=255)                                                                                      ggAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAGGTATGCATTTCAATGCACCACGCTAGTCCTACAAACCTGGc                                                      >  2:137718/1‑144 (MQ=255)                                                                                            tttACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCatga                                           <  1:44003/149‑1 (MQ=255)                                                                                                 aGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCACGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCCCCatgg                                           >  2:228897/1‑143 (MQ=255)                                                                                                   gTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTCCAACAAACCTggacgcgaccataata                                         >  2:208769/1‑130 (MQ=255)                                                                                                     tttttACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATgatgat                                       <  2:242333/141‑1 (MQ=255)                                                                                                        ttACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCAt                                    >  2:47785/1‑144 (MQ=255)                                                                                                         aacacTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCAtt                                   <  2:8766/143‑1 (MQ=255)                                                                                                             cTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTACAATTGCTTCCGTATTCTTTGCAATGCATAAACAGAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCAt                               <  2:343564/144‑1 (MQ=255)                                                                                                              ttATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTa                         >  1:277184/1‑149 (MQ=255)                                                                                                               tATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTAc                        <  1:382311/149‑1 (MQ=255)                                                                                                               tATCGTGCCATAAAATTGACCATCAAGGCAAGATAGCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCa                             >  2:219180/1‑144 (MQ=255)                                                                                                                aTCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTAAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAg                            <  2:227287/144‑1 (MQ=255)                                                                                                                                    aTCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGACGTATCAg   >  1:238742/1‑149 (MQ=255)                                                                                                                                     tCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCAATTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGACGTATCAGt  >  1:265623/1‑149 (MQ=255)                                                                                                                                        agtgaaaaTATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGACGTATCa    <  2:173261/140‑1 (MQ=255)                                                                                                                                              |                                                                                                                                            GCAGTTGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGACGTATCAGT  >  CP009273/2420212‑2420492

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: ‑