Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	662,636	T→G	T3P (ACC→CCC)	ybeD ←	UPF0250 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	662,636	0	T	G	96.1%	176.4 / ‑4.7	52	T3P (ACC→CCC)	ybeD	UPF0250 family protein
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base G (24/25);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTA  >  NC_000913/662539‑662731                                                                                                  |                                                                                                aCCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGgtt                                                                                               <  2:967642/100‑1 (MQ=255)  ccACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGgttt                                                                                              >  2:1227153/1‑100 (MQ=255)   cACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTc                                                                                             <  2:688770/100‑1 (MQ=255)   cACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTc                                                                                             <  2:147117/100‑1 (MQ=255)    aCTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCa                                                                                            <  2:788304/100‑1 (MQ=255)     cTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCAt                                                                                           >  1:73818/1‑100 (MQ=255)        cAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTc                                                                                        >  2:343789/1‑100 (MQ=255)        cAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTc                                                                                        <  1:1563166/100‑1 (MQ=255)         aaCCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCa                                                                                       <  2:262972/100‑1 (MQ=255)          accaccTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAg                                                                                      >  1:646496/1‑100 (MQ=255)          accaccTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAg                                                                                      >  1:1449281/1‑100 (MQ=255)             accTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTc                                                                                   <  2:1241449/100‑1 (MQ=255)                 gATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTCGGTTTCATGTCAGCTCCGGc                                                                               <  2:963063/100‑1 (MQ=255)                    cAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGcgta                                                                            <  2:1575819/100‑1 (MQ=255)                     aaCCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGcgtaa                                                                           <  1:1527370/100‑1 (MQ=255)                         aGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGcgtaacgta                                                                       >  2:698033/1‑100 (MQ=255)                             cAGGTAACGGCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTa                                                                   <  2:695358/100‑1 (MQ=255)                              aGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTa                                                                   >  1:108280/1‑99 (MQ=255)                                gTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAAt                                                                <  2:902901/100‑1 (MQ=255)                                    cGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCa                                                            <  2:450214/100‑1 (MQ=255)                                     gCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCaa                                                           >  1:981270/1‑100 (MQ=255)                                      ccTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAAc                                                          >  2:863003/1‑100 (MQ=255)                                       cTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACt                                                         >  1:119412/1‑100 (MQ=255)                                         nccccATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTcc                                                       >  1:1188312/2‑100 (MQ=255)                                           cccATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCg                                                     >  1:980935/1‑100 (MQ=255)                                             cATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGcc                                                   >  1:369166/1‑100 (MQ=255)                                                aaCTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGa                                                <  1:550569/100‑1 (MQ=255)                                                  cTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAg                                              >  2:1502198/1‑100 (MQ=255)                                                   tttGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAgg                                             <  1:1404778/100‑1 (MQ=255)                                                    ttGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGc                                            <  1:897475/100‑1 (MQ=255)                                                     tGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCg                                           >  2:80576/1‑100 (MQ=255)                                                      tgaaagaaaaGGAGTAGGCAATTCAAGCAGTTCTTTAAGTTTGGTGTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGACCCTCCCGCCAGAAGGCgg                                          <  2:1476932/94‑1 (MQ=255)                                                      gtaagtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCgg                                          <  1:369550/100‑1 (MQ=255)                                                         agtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGGTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAg                                       >  2:451244/1‑100 (MQ=255)                                                             aaaGGAGAAGGGAATTCAAGCAGACCGTGAAGTTTGAATATCATGTCAGCTCCGGCGTGACGTAATTAAATAGCAACTCCTGATAGCAGGCgg                                          >  2:581371/1‑93 (MQ=255)                                                              aaGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGttttt                                  <  1:70506/100‑1 (MQ=255)                                                                      gggAATTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATg                          >  1:1052527/1‑100 (MQ=255)                                                                         aaTTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCtt                       >  1:1548776/1‑100 (MQ=255)                                                                          aTTCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTa                      >  1:1550449/1‑100 (MQ=255)                                                                           ttCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAg                     <  1:863008/100‑1 (MQ=255)                                                                           ttCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAg                     <  1:1096150/100‑1 (MQ=255)                                                                            tCAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGt                    >  1:1542593/1‑100 (MQ=255)                                                                             cAAGCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGta                   >  1:395639/1‑100 (MQ=255)                                                                                gCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTatat                >  1:1288477/1‑100 (MQ=255)                                                                                gCAGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTatat                <  2:73818/100‑1 (MQ=255)                                                                                  aGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATgg              >  2:17234/1‑100 (MQ=255)                                                                                  aGTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATgg              >  1:1015326/1‑100 (MQ=255)                                                                                   gTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGGGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATggg             <  2:1100151/100‑1 (MQ=255)                                                                                   gTTCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATggg             <  2:1038612/100‑1 (MQ=255)                                                                                     tCGTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGa           <  2:646494/100‑1 (MQ=255)                                                                                       gTTAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTGAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACg         <  2:727901/100‑1 (MQ=255)                                                                                         tAAGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGa       <  2:668523/100‑1 (MQ=255)                                                                                          aaGTTTGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGaa      >  1:370039/1‑100 (MQ=255)                                                                                              ttGGGTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTa  >  1:293738/1‑100 (MQ=255)                                                                                                  |                                                                                                ACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTA  >  NC_000913/662539‑662731

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