Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,612,010	C→T	A106T (GCG→ACG)	yneG ←	conserved protein YneG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,612,010	0	C	T	95.3%	129.6 / ‑4.1	43	A106T (GCG→ACG)	yneG	conserved protein YneG
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (19/22);  total (21/22)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.02e-01

ACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGC  >  NC_000913/1611879‑1612146                                                                                                                                    |                                                                                                                                         aCTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGTCACGatat                                                                                                                                  <  2:843783/140‑1 (MQ=255)   tGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGTCACGATATaa                                                                                                                                <  2:62595/140‑1 (MQ=255)     aCCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGTCACGATATAACg                                                                                                                              <  1:247264/140‑1 (MQ=255)     aCCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGTCACGATATAACg                                                                                                                              >  1:661081/1‑140 (MQ=255)     aCCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAACAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGTCACGATATAACg                                                                                                                              <  2:657160/140‑1 (MQ=255)      ccGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGTCACGATATAACgt                                                                                                                             >  2:539869/1‑140 (MQ=255)      ccGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGGACAACCAACCAGTGATGGATTACCGTCACGATATAACgt                                                                                                                             >  2:398847/1‑140 (MQ=255)          aCCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCCGATTCTTCTGTACAGCCAGCCAATGATAGATCACCGTCACGATATAACgttgtt                                                                                                                         >  2:306600/1‑140 (MQ=255)               gATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGttct                                                                                                                    >  1:344507/1‑140 (MQ=255)                  aGCATGTTCGCTGACAGAGGCGGCAAATAAACGTATCCTGGGAATGCTCTTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGTCACGATATAACGTTGTTGTTCTTCg                                                                                                                 <  2:413034/140‑1 (MQ=255)                    cATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGTCACGATATAACGTTGTTGTTCTTCGCt                                                                                                               <  1:1020070/140‑1 (MQ=255)                         tCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCGAGTGATAGGTTAGCGTAACGATATATCGTTGTTGTTCTTAGCTTAACg                                                                                                          >  2:489144/1‑140 (MQ=255)                          cGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGTCACGATATAACGTTGTTGTTCTTCGCTTAACGa                                                                                                         >  2:697851/1‑140 (MQ=255)                           gCTGACAGCGGCGGCCAATAAAAGTCCCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCCGATTCATCTGTACAACCAGCCAGTGATAGATTACCGTCACGATATAACGTTGTTGTTCTTCGCTTAACGaa                                                                                                        <  1:306600/140‑1 (MQ=255)                               aCAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGTCACGATATAACGTTGTTGTTCTTCGCTTAACGAAacac                                                                                                    <  2:661081/140‑1 (MQ=255)                                cAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGTCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACAcc                                                                                                   <  2:207422/140‑1 (MQ=255)                                     gcggcAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGTCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCg                                                                                              <  1:856254/140‑1 (MQ=255)                                          aaaTAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGTCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATg                                                                                         <  1:539869/140‑1 (MQ=255)                                                     cccTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGTCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACtt                                                                              >  2:718403/1‑140 (MQ=255)                                                     cccTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGTCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACtt                                                                              >  1:125013/1‑140 (MQ=255)                                                        tGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGTCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGc                                                                           >  2:562841/1‑140 (MQ=255)                                                          ggAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGTCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTAt                                                                                     <  2:956771/128‑1 (MQ=255)                                                          ggAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGTCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTAt                                                                                     >  1:956771/1‑128 (MQ=255)                                                              tGCTATTTGGTCGTTTGGCGCAATTCTCATGGCTGGTTCCTCTGTACACCCAGCCAGTGATAGATTACCGTCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACa                                                                     <  2:457598/140‑1 (MQ=255)                                                                       gTCGTTTTCCGCAATTCTCATGCCTTATTTATCTGTACAACCAGCCAGTGATAGATTACCGTCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGAca                                                            <  1:398847/140‑1 (MQ=255)                                                                        tCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGTCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGAcag                                                           >  2:175196/1‑140 (MQ=255)                                                                          gTTTTCCGCAATTCTCACGGCTAATTCATCTTTACAACCAGCCAGTGATAGATTACCGTCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGAcagca                                                         <  2:344507/140‑1 (MQ=255)                                                                                      tctcATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGTCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGc                                             <  2:595483/140‑1 (MQ=255)                                                                                           tGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGTCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCt                                        <  1:235100/140‑1 (MQ=255)                                                                                           tGGCTGATTCATCTGTACAACCAGCCAGTGACAGATTACTGTCACGATTTAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCt                                        <  2:295743/140‑1 (MQ=255)                                                                                              cTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGTCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAg                                     <  1:235174/140‑1 (MQ=255)                                                                                                  ttCATCTGTACAACCAGCCAGTGATAGATTACCGTCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGa                                  >  2:216221/1‑139 (MQ=255)                                                                                                  ttCATCTGTACAACCAGCCAGTGATAGATTACCGTCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGa                                  <  1:216221/139‑1 (MQ=255)                                                                                                       cTGTACAACCAGCCAGTGATAGATTACCGTCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCACGACAGCAGGTGGCGGTGGCGTGCTGAGCGATAAaca                            >  1:290196/1‑140 (MQ=255)                                                                                                          tACAACCAGCCAGTGATAGATTACCGTCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCg                         >  1:725608/1‑140 (MQ=255)                                                                                                                 agccagTGATAGATTACCGTCAAGATATAACCTTGCTGTCCTTCGCTTCACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGccc                  <  2:290196/140‑1 (MQ=255)                                                                                                                  gccagTGATAGATTACCGTCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCcg                 >  2:132146/1‑140 (MQ=255)                                                                                                                   ccagTGATAGATTACCGTCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCcgc                >  2:512026/1‑140 (MQ=255)                                                                                                                    cagTGATAGATTACCGTCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCcgcg               >  1:999496/1‑140 (MQ=255)                                                                                                                       tGATAGATTACCGTCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGATCGATAAACAACGGATGCCCGCGCAt            >  1:406148/1‑140 (MQ=255)                                                                                                                            gATTACCGTCACGATTTAACGTCGTGGTTTTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATggggg       <  1:371318/140‑1 (MQ=255)                                                                                                                             aTTACCGTCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGt      <  1:562841/140‑1 (MQ=255)                                                                                                                               tACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGttt    >  1:944403/1‑140 (MQ=255)                                                                                                                                 ccGTCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGc  >  1:544675/1‑140 (MQ=255)                                                                                                                                    |                                                                                                                                         ACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGC  >  NC_000913/1611879‑1612146

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