Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	216,568	C→T	A146T (GCC→ACC)	yaeF ←	putative lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	216,568	0	C	T	94.3%	166.7 / ‑1.2	53	A146T (GCC→ACC)	yaeF	putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (32/18);  total (34/19)
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

GAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGCAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAA  >  NC_000913/216473‑216657                                                                                                |                                                                                          gaACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAAc                                                                                       >  1:45320/1‑100 (MQ=255)    cGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAAcgcg                                                                                    >  2:588352/1‑100 (MQ=255)    cGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACAGCTATCTTTGATTTTATTGGTAAAcgcg                                                                                    >  2:181811/1‑100 (MQ=255)      gATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGt                                                                                  >  1:372542/1‑100 (MQ=255)              gAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCtg                                                                          <  2:649627/100‑1 (MQ=255)                gCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCtgtg                                                                        <  2:967802/100‑1 (MQ=255)                gCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCtgtg                                                                        <  2:349177/100‑1 (MQ=255)                gCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCtgtg                                                                        >  2:7347/1‑100 (MQ=255)                 cacaTCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTgg                                                                       <  1:145804/100‑1 (MQ=255)                  acaTCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGc                                                                      <  2:1283453/100‑1 (MQ=255)                   caTCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCt                                                                     >  2:422557/1‑100 (MQ=255)                     tatGGCGAGTCACGATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCttg                                                                   >  1:879468/3‑100 (MQ=255)                     tCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCttg                                                                   >  2:237933/1‑100 (MQ=255)                      cTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCttgt                                                                  >  1:16053/1‑100 (MQ=255)                       tGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCttgt                                                                  <  2:1271104/99‑1 (MQ=255)                        ggCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCttgttg                                                                >  1:274223/1‑100 (MQ=255)                        ggCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCttgttg                                                                >  2:810997/1‑100 (MQ=255)                           gAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCgg                                                             <  1:678831/100‑1 (MQ=255)                                        gggAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCggg                                                >  1:498961/1‑100 (MQ=255)                                        gggAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCggg                                                <  1:157064/100‑1 (MQ=255)                                         ggAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGa                                               >  1:1100246/1‑100 (MQ=255)                                          gAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGCAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGAc                                              >  1:308248/1‑100 (MQ=255)                                             tAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCg                                           >  1:1132977/1‑100 (MQ=255)                                             tAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCg                                           >  1:710822/1‑100 (MQ=255)                                             tAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCg                                           >  1:975305/1‑100 (MQ=255)                                             tAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCg                                           >  1:1082314/1‑100 (MQ=255)                                               aaTTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTa                                         <  1:685165/100‑1 (MQ=255)                                                      cAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGaa                                  >  1:378188/1‑100 (MQ=255)                                                      cAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGaa                                  >  1:736145/1‑100 (MQ=255)                                                        aTGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGaaaa                                <  1:970354/100‑1 (MQ=255)                                                        aTGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGaaaa                                >  1:1155272/1‑100 (MQ=255)                                                                   gTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTa                     >  1:382066/1‑100 (MQ=255)                                                                    ttATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTAt                    <  2:417280/100‑1 (MQ=255)                                                                    ttATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTAt                    <  1:95259/100‑1 (MQ=255)                                                                     tataACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATg                   <  1:474286/100‑1 (MQ=255)                                                                       taACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTGTCACTATGCt                 >  2:1261869/1‑100 (MQ=255)                                                                       taACCGCTATCTTTGATTTTATTGGCAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCt                 <  2:869239/100‑1 (MQ=255)                                                                        aaCCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTTAGGGGTAAGATGCGGGACTCGTAAGACGAAAAGCTTATCACTATGCtt                >  2:880136/1‑100 (MQ=255)                                                                        aaCCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGTGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCt                 <  2:99083/99‑1 (MQ=255)                                                                        aaCCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCtt                >  2:390315/1‑100 (MQ=255)                                                                         aCCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCtt                <  2:631941/99‑1 (MQ=255)                                                                          ccGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCa              >  1:668030/1‑100 (MQ=255)                                                                          ccGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCa              >  1:1303034/1‑100 (MQ=255)                                                                          ccGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCa              >  2:1146500/1‑100 (MQ=255)                                                                            gCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATc            <  1:379298/100‑1 (MQ=255)                                                                               aTCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCt         <  1:422561/100‑1 (MQ=255)                                                                                 cTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCttt       <  1:962198/100‑1 (MQ=255)                                                                                 cTTTGATTTTATTGGCAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCttt       >  2:368585/1‑100 (MQ=255)                                                                                  tttGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCtttt      >  2:1261410/1‑100 (MQ=255)                                                                                   ttGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTACGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCttttt     >  2:29623/1‑100 (MQ=255)                                                                                    tGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCtttttt    >  2:66024/1‑100 (MQ=255)                                                                                    tGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCtttttt    >  1:113387/1‑100 (MQ=255)                                                                                      aTTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTaa  >  1:592103/1‑100 (MQ=255)                                                                                                |                                                                                          GAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGCAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAA  >  NC_000913/216473‑216657

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