Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	355,454	G→A	G178D (GGC→GAC)	codB →	cytosine transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	355,454	0	G	A	97.0%	95.0 / ‑4.0	33	G178D (GGC→GAC)	codB	cytosine transporter
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (19/13);  total (19/14)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.05e-01

GTTTCCGGTTTACTGATGACCGTCACCGTCTTTTTTGGCATTTCGGCGCTGACGGTTCTTTCGGTGATTGCGGTTCCGGCTATCGCCTGCCTGGGCGGTTATTCCGTGTGGCTGGCTGTTAACGGCATGGGCGGCCTGGACGCATTAAAAGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTCATTTATCAGTGCGGGTACGCTCACCGCTGACTTTGTCCGGTTTGGTCGCAATGCCAAACTG  >  NC_000913/355321‑355590                                                                                                                                      |                                                                                                                                         gTTTCCGGTTTACTGATGACCGTCACCGTCTTTTTTGGCATTTCGGCGCTGACGGTTCTTTCGGTGATTGCGGTTCCGGCTATCGCCTGCCTGGGCGGTTATTCCGTGTGGCTGGCTGTTAACGGCATGGGCGGCCTGGa                                                                                                                                    <  2:117933/140‑1 (MQ=255)  tttCCGGTTTACTGATGACCGCCACCGTCTTTTTTGGCATTTCGGCGCTGACGGTTCTTTCGGTGATTGCGGTTCCGGGTATCGCCTGCCTGGGCGGTTATTCCGTGTGGCTGGCTGTTAACGGCCTGGGTGCCCTGGAc                                                                                                                                   >  2:56769/1‑140 (MQ=255)           tACTGATGACCGTCACCGTCTTTTTTGGCATTTCGGCGCTGACGGTTCTTTCGGTGATTGCGGTTCCGGCTATCGCCTGCCTGGGCGGTTATTCCGTGTGGCGGGCTGTTTACGGCATGGGCGAGCTGGACGCATTaaaa                                                                                                                          >  1:215716/1‑140 (MQ=255)               gatgaCCGTCACCGTCTTTTTTGGCATTTCGGCGCTGACGGTTCTTTCGGTGATTGCGGTTCCGGCTATCGCCTGCCTGGGCGGTTATTCCCTGTGGCTGACTGTTAAAGGCATGGGGGACCTGGACGCATTAAAAGCgg                                                                                                                      >  2:199804/1‑140 (MQ=255)                   accgtcaccgtcTTTTTGGGCATTTCGGCGCAGACGGTTGTTTCGGTGCTTGCGGTTCCGGCTATCGCCTGCCTGGGCGGTTATTCCGTGTGGCTGGCTGTTAACGGCATGGGCGACCTGGACGCATTAAAAGCGGTCGt                                                                                                                  <  2:223863/140‑1 (MQ=255)                     cgtcaccgtcTTTTTTGGCATTTCGGCGCTGACGGTTCTTTCGGTGATTGCGGTTCCGGCTATCGCCTGCCTGGGCGGTTATTCCGTGTGGCTGGCTGTTAACGGCATGGGCGACCTGGACGCATTAAAAGCGGTCGTTc                                                                                                                >  2:637446/1‑140 (MQ=255)                       tcaccgtcTTTTTTGGCATTTCGGCGCTGACGGTTCTTTCGGTGATTGCGGTTCCGGCTATCGCCTGCCTGGGCGGTTATTCCGTGTGGCTGGCTGTTAACGGCATGGGCGACCTGGACGCATTAAAAGCGGTCGTTccc                                                                                                              >  1:600962/1‑140 (MQ=255)                                                        ttCTTTCGGTGATTGCGGTTCCGGCTATCGCCTGCCTGGGCGGTTATTCCGTGTGGCTGGCTGTTAACGGCATGGGCGACCTGGACGCATTAAAAGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCgcgctggcg                                                                             <  1:442512/140‑1 (MQ=255)                                                         tCTTTCGGTGATTGCGGTTCCGGCTATCGCCTGCCTGGGCGGTTATTCCGTGTGGCTGGCTGTTAACGGCATGGGCGACCTGGACGCATTAAAAGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCgcgctggcgc                                                                            <  1:33491/140‑1 (MQ=255)                                                         tCTTTCGGTGATTGCGGTTCCGGCTATCGCCTGCCTGGGCGGTTATTCCGTGTGGCTGGCTGTTAACGGCATGGGCGACCTGGACGCATTAAAAGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCgcgctggcgc                                                                            <  1:47922/140‑1 (MQ=255)                                                              cGGTGATTGCGGTTCCGGCTATCGCCTGCCTGGGCGGTTATTCCGTGTGGCTGGCTGTTAACGGCATGGGCGACCTGGACGCATTAAAAGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGtt                                                                       >  2:594326/1‑140 (MQ=255)                                                                gTGATTGCGGTTCCGGCTATCGCCTGCCTGGGCGGTTATTCCGTGTGGCTGGCTGTTAACGGCATGGGCGACCTGGACGCATTAAAAGCGGTCGTTCCCGCACAACCGTTTGATATCCATGTCGCGCTAGCGCTGGTtgt                                                                     >  1:295147/1‑140 (MQ=255)                                                                  gATTGTGGTTCCGGCTATCGCCTGCCTGGGCGGTTATTCCGTGTGGCTGGCTGTTAACGGCATGGGCGACCTGGACGCATTAAAAGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGTTGTgg                                                                   <  2:350342/140‑1 (MQ=255)                                                                  gATTGCGGTTCCGGCTACCGCCTGCCTGGGCAGTTATTCTGTGTGGCTTGCTGTTAACGGCATGGGCGACCTGGACGCATTAAAAGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTAGCGCTGGTTGTgg                                                                   <  2:295147/140‑1 (MQ=255)                                                                    ttGCGGTTCCGGCTATCGCCTGCCTGGGCGGTTATTCCGTGTGGCTGGCTGTTAACGGCATGGGCGACCTGGACGCATTAAAAGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGTTGTgggg                                                                 >  1:580293/1‑140 (MQ=255)                                                                    ttGCGGTTCCGGCTATCGCCTGCCTGGGCGGTTATTCCGTGTGGCTGGCTGTTAACGGCATGGGCGACCTGGACGCATTAAAAGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGTTGTgggg                                                                 <  2:654020/140‑1 (MQ=255)                                                                            ccGGCTATCGCCTGCTTGGGCGATTCTTCAGTGTGGCTGGCTGTTAACGGCATGGGCGACCTGGACGCATTAAAAGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTCATTTAt                                                         <  1:199804/140‑1 (MQ=255)                                                                            ccGGCTATCGCCTGCCTGGGCGGTTATTCCGTGTGGCTGGCTGTTAACGGCATGGGCGACCTGGACGCATTACAAGCGGTCGTTCCCGCACAACCGTTAGACTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTCATTTAt                                                         <  2:359795/140‑1 (MQ=255)                                                                            ccGGCTATCGCCTGCCTGGGCGGTTATTCCGTGTGGCTGGCTGTTAACGGCATGGGCGACCTGGACGCATTAAAAGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTCATTTAt                                                         <  1:479055/140‑1 (MQ=255)                                                                               gCTATCGCCTGCCTGGGCGGTTATTCCGTGTGGCTGGCTGTTAACGGCATGGGCGACCTGGACGCATTAAAAGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTCATTTATCAg                                                      <  1:56769/140‑1 (MQ=255)                                                                               gCTATCGCCTGCCTGGGCGGTTATTCCGTGTGGCTGGCTGTTAACGGCATGGGCGACCTGGACGCATTAAAAGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTCATTTATCAg                                                      >  1:583073/1‑140 (MQ=255)                                                                                cTATCGCCTGCCTGGGCGGTTATTCCGTGTGGCTGGCTGTTAACGGCATGGGCGACCTGGACGCATTAAAAGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTCATTTATCAGt                                                     <  1:637446/140‑1 (MQ=255)                                                                                          cctgGGCGGATATTCCGTGTGGCTGGCTGTTAACGGCATGGGCGACCTGGACGCATTAAAAACGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTCATTTATCAGTGCCGGGACGc                                           >  2:40680/1‑140 (MQ=255)                                                                                             gggCGGTTATTCCGTGTGGCTGGCTGTTAACGGCATGGGCGACCTGGACGCATTAAAAGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTCATTTATCAGTGCGGGTACGCTCa                                        >  2:584535/1‑140 (MQ=255)                                                                                               gCGGTTATTCCGTGTGGCTGGCTGTTAACGGCATGGGCGACCTGGACGGATTAAAAGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCGTGCTGGCGCGGGTTGTGGGGTGATTTATCAGAGCGGGTACGCTCAcc                                      >  1:292669/1‑140 (MQ=255)                                                                                                    tATTCCGTGTGGCTGGCTGTTAACGGCATGGGCGACCTGGACGCATTAAAAGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTCATTTATCAGTGCGGGTACGCTCACCGCTGa                                 >  1:454969/1‑140 (MQ=255)                                                                                                    tATTCCGTGTGGCTGGCTGTTAACGGCATGGGCGACCTGGACGCATTAAAAGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTCATTTATCAGTGCGGGTACGCTCACCGCTGa                                 >  1:30633/1‑140 (MQ=255)                                                                                                    tATTCCGTGTGGCTGGCTGTTAACGGCATGGGCGACCTGGACGCATTAAAAGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTCATTTATCAGTGCGGGTACGCTCACCGCTGa                                 <  1:459133/140‑1 (MQ=255)                                                                                                    tATTCCGTGTGGCTGGCTGTTAACGGCAGGGATGATTTAGGCGCATTAAAAGCGGTCGTTCCCGCCCAATCTTTAGATTTCAATGTCGCGCTCGCTTTAGTAGTGGGGTCATTTATCAGTGCGGGTACGCTCACCGCTGa                                 <  2:485024/140‑1 (MQ=255)                                                                                                            gtgGCTGGCTGTTAACGGCATGGGCGACCTGGACGCATTAAAAGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTCATTTATCAGTGCGGGTACGCTCACCGCTGACTTTGTcc                         >  1:457590/1‑140 (MQ=255)                                                                                                             tggctggcTGTTAACGGCATGGGCGACCTGGACGCATTAAAAGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTCATTTATCAGTGCGGGTACGCTCACCGCTGACTTTGTCCg                        >  1:117075/1‑140 (MQ=255)                                                                                                                              cATGGGCGACCTGGACGCATTAAAAGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTCATTTATCAGTGCGGGTGCGCGCACCGCTGACTTTGTCCGGTTTGGTCGCAATGCCa       >  2:266840/1‑140 (MQ=255)                                                                                                                               aTGGGCGACCTGGACGCATTAAAAGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTCATTTATCAGTGCGGGTACGCTCACCGCTGACTTTTTCCGGTTTGGTCGCAATGCCaa      >  2:211549/1‑140 (MQ=255)                                                                                                                                 gggCGACCTGGACGCATTAAAAGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTCATTTATCAGTGCGGGTACGCTCACCGCTGACTTTGTCCGGTTTGGTCGCAATGCCAAAc    >  2:565625/1‑140 (MQ=255)                                                                                                                                 gggCGACCTGGACGCATTAAAAGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTCATTTATCAGTGCGGGTACGCTCACCGCTGACTTTGTCCGGTTTGGTCGCAATGCCAAAc    >  1:364025/1‑140 (MQ=255)                                                                                                                                 gggCGACCTGGACGCATTAAAAGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTCATTTATCAGTGCGGGTACGCTCACCGCTGACTTTGTCCGGTTTGGTCGCAATGCCAAAc    >  1:518843/1‑140 (MQ=255)                                                                                                                                 gggCGACCTGGACGCATTAAAAGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTCATTTATCAGTGCGGGTACGCTCACCGCTGACTTTGTCCGGTTTGGTCGCAATGCCAAAc    >  1:604335/1‑140 (MQ=255)                                                                                                                                   gcgacCTGGACGCATTAAAAGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTCATTTATCAGTGCGGGTACGCTCACCGCTGACTTTGTCCGGTTTGGTCGCAATACCAAACTg  >  2:626252/1‑140 (MQ=255)                                                                                                                                      |                                                                                                                                         GTTTCCGGTTTACTGATGACCGTCACCGTCTTTTTTGGCATTTCGGCGCTGACGGTTCTTTCGGTGATTGCGGTTCCGGCTATCGCCTGCCTGGGCGGTTATTCCGTGTGGCTGGCTGTTAACGGCATGGGCGGCCTGGACGCATTAAAAGCGGTCGTTCCCGCACAACCGTTAGATTTCAATGTCGCGCTGGCGCTGGTTGTGGGGTCATTTATCAGTGCGGGTACGCTCACCGCTGACTTTGTCCGGTTTGGTCGCAATGCCAAACTG  >  NC_000913/355321‑355590

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