Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,178,061	C→G	11.1%	intergenic (‑236/‑138)	ribB ← / → yqiC	3,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase/BMFP family putative fusogenic protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,178,061	0	C	G	11.1%	101.8 / 5.9	36	intergenic (‑236/‑138)	ribB/yqiC	3,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase/BMFP family putative fusogenic protein
Reads supporting (aligned to +/- strand):  ref base C (16/16);  new base G (2/2);  total (18/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.41e-01

AGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGTTACTCTCTCCCATCCGGACTCTAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCTTCACATAATTCTGGTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAA  >  CP009273/3177864‑3178304                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    aGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGTTACTCTCTCCCATCCGGACTCTAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCTTCACATAATTCTGGttt                                                                                                                                                                        <  1:602050/275‑1 (MQ=255)                  ggcggcATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGTTACTCTCTCCCATCCGGACTCTAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGAAGTTCGCCCCGCCCTGAGAATAAGCg                                                                                                                                                                                                                                    >  2:489982/1‑198 (MQ=255)                      gcATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGTTACTCTCTCCCATCCGGACTCTAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGAACTTCGCCCCGCCCTGAGAATAAGCGGAtt                                                                                                                                                                                                                                <  1:489990/198‑1 (MQ=255)                      gcATATAACCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGTTACTCTCTCCCATCCGGACTCTAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGAATTTCGTCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCt                                                                                                                                                                                         <  1:390818/237‑1 (MQ=255)                         tataGCCAAAATAACGTGAGCGGGTCCGTGCCCGACAGAATCGTTACTCTCTCCCATCCGGACTCTAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCTTCACATAATTATGGtta                                                                                                                                                                        >  1:397994/1‑250 (MQ=255)                                    tAACGTGAGCGGGTCCATGCCCGACAGAATCGTTACTCTCTCCCATCCGGACTCTAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGAATTTCGCCCCGCCCTGAGAATAAGCg                                                                                                                                                                                                                                    >  1:33121/1‑180 (MQ=255)                                        gTGAGCGGGTCCATGCCCGACAGAATCGTTACTCTCTCCCATCCGGACTCTAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGAtt                                                                                                                                                                                                                                <  2:33121/180‑1 (MQ=255)                                                      gCCCGCCAGATTCGTTACTCTCTCCCATTCGGACTCTAACCGACGGCCTCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGAATTTCGCCGCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCTTCACATAATTCTGGTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGa                                                                                                                   <  1:4736/275‑1 (MQ=255)                                                              gAATCGTTACTCTCTCCCATCCGGACTCTAACCGTCCGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGAATTTCGCCGCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCTGCCGCAACGCATAGCCTCAAATAATTCAGGTTTATGACTTCCCCTTACCGCACTACAATGGCACCCAACCCCCTACAGTACAGGGACCCTActa                                                                                                           >  2:4736/1‑273 (MQ=255)                                                                    ttACTCTCTCCCATCCGGACTCTAACCGTCGGCCCCGAAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCTGTGGGTAATTTCCCCCCGCCCTTAGAATAAGCGGATTCACTATAAc                                                                                                                                                                                                                       >  2:145474/1‑161 (MQ=255)                                                                        tctctcCCATCCGGACTCTAACCGTCGGCCCCGAAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTa                                                                                                                                                                                                                   <  1:145477/161‑1 (MQ=255)                                                                            tcCCATCCGGACTCTAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGAATTTCGCCCCGCCCTGAGAATAAg                                                                                                                                                                                                                                      >  1:403854/1‑138 (MQ=255)                                                                                aTCCGGACTCTAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGAATTTCGCCCCGCCCTGAGAATCAGCGGa                                                                                                                                                                                                                                  <  2:403849/138‑1 (MQ=255)                                                                                            aCCGTCGGCCCCGGAATTACACCGTATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGAATTTCGCCCCGCCATGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCTTCACATAATTCTGGTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAAtcg                                                                             <  1:345733/275‑1 (MQ=255)                                                                                                                    gATCTGCTGTCCTTTGAGTTCGCCCCCAGTGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCtt                                                                                                                                                                                        >  1:468568/1‑144 (MQ=255)                                                                                                                    gATCTGCTGTCCTTTGAGTTCGCACCCAAAGTGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGAATTTCGCCCCGCCCTGAGACTAAGAGGATTAACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCTTCACATAATTCTAGTTTATGCATTACCCTTATCGCACTACNATGGCACTCAacac                                                                                                                                  >  1:510106/1‑198 (MQ=255)                                                                                                                    gATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGAATTTCGCCCCGCCCTGAGAATAAGCg                                                                                                                                                                                                                                    >  2:209515/1‑100 (MQ=255)                                                                                                                        tgctgTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCATCTGAGTTGATTTACCGCCAGTGGGGAATGTCGCTCCGCCTTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCTTcaca                                                                                                                                                                                    <  2:468561/144‑1 (MQ=255)                                                                                                                        tgctgTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGAtt                                                                                                                                                                                                                                <  1:209519/100‑1 (MQ=255)                                                                                                                              ccTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGAATTTCGCCCCGCCCTGAGAATAAGCg                                                                                                                                                                                                                                    >  1:256479/1‑90 (MQ=255)                                                                                                                                  tGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGAtt                                                                                                                                                                                                                                <  2:256475/90‑1 (MQ=255)                                                                                                                                       tCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGAATTTcgccc                                                                                                                                                                                                                                                     >  1:91795/1‑64 (MQ=255)                                                                                                                                       tCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGAATTTCGCCCCGCCCTGAGaa                                                                                                                                                                                                                                          >  2:253390/1‑75 (MQ=255)                                                                                                                                           aCCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGAATTTcgccccgcc                                                                                                                                                                                                                                                 <  2:91795/64‑1 (MQ=255)                                                                                                                                           aCCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGAATTTCGCCCCGCCCTGAGAATAAg                                                                                                                                                                                                                                      <  1:253394/75‑1 (MQ=255)                                                                                                                                               aaaGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCTTCACATAATTCTg                                                                                                                                                                            >  2:60293/1‑129 (MQ=255)                                                                                                                                                   cgcTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCTTCACATAATTCTGGttt                                                                                                                                                                        <  1:60293/129‑1 (MQ=255)                                                                                                                                                     cTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGAATTTCGCCCCGCCCTGAGAATAAGc                                                                                                                                                                                                                                     >  1:568653/1‑66 (MQ=255)                                                                                                                                                       cgcgGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCTTCACATAATTCTGGTTTATGACTTACCCTTATCGCACTACAATgg                                                                                                                                            >  2:532280/1‑153 (MQ=255)                                                                                                                                                         cgGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGAt                                                                                                                                                                                                                                 <  2:568642/66‑1 (MQ=255)                                                                                                                                                           ggCTTTCAACTGAGTTGATTTACCGCCGGTGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCTTCACATAATTCTGGTTTATGACTTACCCTTATCGCACTACAATGGCACt                                                                                                                                        <  1:532291/153‑1 (MQ=255)                                                                                                                                                             cTTTCAACTGAGTTGATTTACCGCCGGTGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCTTCACATAATTCTGGTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGCAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTAGGAGAAGATGTGGAGaaaaaa            >  1:383702/1‑275 (MQ=255)                                                                                                                                                                      gAGTTGATTTACCGCCGGTGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCTTCACATAATTCTGGTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGGAAGATATGGAGAAAAAAATCCGc      >  1:403849/1‑272 (MQ=255)                                                                                                                                                                        gTTGATTTACCGCCGGTGGGGAATTTCGCCGCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCtt                                                                                                                                                                                        >  1:553911/1‑92 (MQ=255)                                                                                                                                                                         tggatttaagacaggTGGGGAANATCGCCCCGCCCTGAGAATAAAAGAATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCTTCACATAATTCTGGTTTATGACTTACCCTTATCGCACTACAATGGCACTAAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGGAAGATATGGAGAAAAAAATCAGCCaaa  <  2:403844/260‑1 (MQ=255)                                                                                                                                                                            aTTTACCGCCGGTGGGGAATTTCGCCGCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCTTcaca                                                                                                                                                                                    <  2:553900/92‑1 (MQ=255)                                                                                                                                                                                 ccgccgGTGGGGAATTTCGCCCCGCCATGAGACTAAGCGGCTTAACTAAAAc                                                                                                                                                                                                                       >  2:600873/1‑52 (MQ=255)                                                                                                                                                                                     cgGTGGGGAATTTCGCCCCGCCATGAGAATAAGCGGATTAACTATAACGCTa                                                                                                                                                                                                                   <  1:600884/52‑1 (MQ=255)                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    AGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTCCATGCCCGACAGAATCGTTACTCTCTCCCATCCGGACTCTAACCGTCGGCCCCGGAATTACACCGGATCTGCTGTCCTTTGAGTTCGCACCCAAAGCGCTCGCGGGCTTTCAACTGAGTTGATTTACCGCCGGTGGGGAATTTCGCCCCGCCCTGAGAATAAGCGGATTCACTATAACGCTAATGATTAGCGGCAGCAACGCATAGCTTCACATAATTCTGGTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAA  >  CP009273/3177864‑3178304

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑