Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	432,771	T→G	R27R (CGA→CGC)	yajI ←	putative lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	432,771	0	T	G	98.2%	192.3 / ‑3.6	54	R27R (CGA→CGC)	yajI	putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (25/28);  total (25/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

ATTTTGCATTCAGCCTGTTTTGCTGAGTGATTTTGACTGTTTCTTGATTGAGTTGCGTCATCTCCTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGTCGCACTTCACTTTGCTGCACACAGGCGCTAAGGCTGAACAGGCTTCCCAGCAGGAGCAGTCGAAAAACGTTAGTATTCATAGGTCGTCCCTTT  >  NC_000913/432674‑432864                                                                                                  |                                                                                              aTTTTGCATTCAGCCTGTTTTGCTGAGTGATTTTGACTGTTTCTTGATTGAGTTGCGTCATCTCCTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGGCg                                                                                             <  2:258581/100‑1 (MQ=255)    ttGCATTCAGCCTGTTTTGCTGAGTGATTTTGACTGTTTCTTGATTGAGTTGCGTCATCTCCTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGGCGcac                                                                                          <  1:528213/100‑1 (MQ=255)       cATTCAGCCTGTTTTGCTGAGTGATTTTGACTGTTTCTTGATTGAGTTGCGTCATCTCCTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGGCGcacttc                                                                                       >  2:707469/1‑100 (MQ=255)          tCAGCCTGTTTTGCTGAGTGATTTTGACTGTTTCTTGATTGAGTTGCGTCATCTCCTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGGCGcacttcact                                                                                    <  2:911162/100‑1 (MQ=255)            aGCCTGTTTTGCTGAGTGATTTTGACTGTTTCTTGATTGAGTTGCGTCATCTCCTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGGCGCACTTCACttt                                                                                  <  1:541687/100‑1 (MQ=255)            aGCCTGTTTTGCTGAGTGATTTTGACTGTTTCTTGATTGAGTTGCGTCATCTCCTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGGCGCACTTCACttt                                                                                  <  1:112095/100‑1 (MQ=255)             gCCTGTTTTGCTGAGTGATTTTGACTGTTTCTTGATTGAGTTGCGTCATCTCCTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGGCGCACTTCACTTtg                                                                                 >  2:1125560/1‑100 (MQ=255)             gCCTGTTTTGCTGAGTGATTTTGACTGTTTCTTGATTGAGTTGCGTCATCTCCTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGGCGCACTTCACTTtg                                                                                 >  1:476188/1‑100 (MQ=255)              ccTGTTTTGCTGAGTGATTTTGACTGTTTCTTGATTGAGTTGCGTCATCTCCTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGGCGCACTTCACTTtgc                                                                                <  1:699834/100‑1 (MQ=255)              ccTGTTTTGCTGAGTGATTTTGACTGTTTCTTGATTGAGTTGCGTCATCTCCTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGGCGCACTTCACTTtgc                                                                                <  1:150259/100‑1 (MQ=255)               cTGTTTTGTTGAGTGATTTTGACTGTTTCTTGATTGAGTTGCGTCATCTCCTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGGCGCACTTCACTTtgct                                                                               <  1:306454/100‑1 (MQ=255)               cTGTTTTGCTGAGTGATTTTGACTGTTTCTTGATTGAGTTGCGTCATCTCCTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGGCGCACTTCACTTtgct                                                                               <  1:1228849/100‑1 (MQ=255)               cTGTTTTGCTGAGTGATTTTGACTGTTTCTTGATTGAGTTGCGTCATCTCCTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGGCGCACTTCACTTtgct                                                                               <  2:1303008/100‑1 (MQ=255)                 gTTTTGCTGAGTGATTTTGACTGTTTCTTGATTGAGTTGCGTCATCTCCTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGGCGCACTTCACTTtgctgc                                                                             <  1:7454/100‑1 (MQ=255)                  ttttGCTGAGTGATTTTGACTGTTTCTTGATTGAGTTGCGTCATCTCCTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGGCGCACTTCACTTTGCTGca                                                                            <  2:761762/100‑1 (MQ=255)                  ttttGCTGAGTGATTTTGACTGATTCTTGATTGAGTTCCGTCATCTCCTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGGCGCACTTCACTTTGCTGca                                                                            <  2:878413/100‑1 (MQ=255)                      gCTGAGTGATTTTGACTGTTTCTTGATTGAGTTGCGTCATCTCCTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGGCGCACTTCACTTTGCTGcacaca                                                                        <  2:289553/100‑1 (MQ=255)                             gATTTTGACTGTTTCTTGATTGAGTTGCGTCATCTCCTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGGCGCACTTCACTTTGCTGCACACAGGCGCTa                                                                 <  1:467627/100‑1 (MQ=255)                              aTTTTGACTGTTTCTTGATTGAGTTGCGTCATCTCCTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGGCGCACTTCACTTTGCTGCACACAGGCGCTaa                                                                >  1:914760/1‑100 (MQ=255)                              aTTTTGACTGTTTCTTGATTGAGTTGCGTCATCTCCTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGGCGCACTTCACTTTGCTGCACACAGGCGCTaa                                                                >  1:914580/1‑100 (MQ=255)                                tttGACTGTTTCTTGATTGAGTTGCGTCATCTCCTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGTCGCACTTCACTTTGCTGCACACAGGCGCTAAgg                                                              <  2:410494/100‑1 (MQ=255)                                tttGACTGTTTCTTGATTGAGTTGCGTCATCTCCTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGGCGCACTTCACTTTGCTGCACACAGGCGCTAAgg                                                              <  2:962025/100‑1 (MQ=255)                                 ttGACTGTTTCTTGATTGAGTTGCGTCATCTCCTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGGCGCACTTCACTTTGCTGCACACAGGCGCTAAGGc                                                             <  1:876938/100‑1 (MQ=255)                                      tGTTTCTTGATTGAGTTGCGTCATCTCCTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGGCGCACTTCACTTTGCTGCACACAGGCGCTAAGGCTGAAc                                                        <  2:410383/100‑1 (MQ=255)                                         ttCTTGATTGAGTTGCGTCATCTCCTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGGCGCACTTCACTTTGCTGCACACAGGCGCTAAGGCTGAACAgg                                                     >  1:977056/1‑100 (MQ=255)                                          tCTTGATTGAGTTGCGTCATCTCCTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGGCGCACTTCACTTTGCTGCACACAGGCGCTAAGGCTGAACAGGc                                                    >  1:393976/1‑100 (MQ=255)                                          tCTTGATTGAGTTGCGTCATCTCCTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGGCGCACTTCACTTTGCTGCACACAGGCGCTAAGGCTGAACAGGc                                                    >  2:228908/1‑100 (MQ=255)                                              gattgaGTTGCGTCATCTCCTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGGCGCACTTCACTTTGCTGCACACAGGCGCTAAGGCTGAACAGGCTTcc                                                >  2:972189/1‑100 (MQ=255)                                                  gaGTTGCGTCATCTCCTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGGCGCACTTCACTTTGCTGCACACAGGCGCTAAGGCTGAACAGGCTTCCcagc                                            <  2:409742/100‑1 (MQ=255)                                                   aGTTGCGTCATCTCCTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGGCGCACTTCACTTTGCTGCACACAGGCGCTAAGGCTGAACAGGCTTCCcagca                                           >  2:248272/1‑100 (MQ=255)                                                    gTTGCGTCATCTCCTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGGCGCACTTCACTTTGCTGCACACAGGCGCTAAGGCTGAACAGGCTTCCcagcag                                          >  2:1263152/1‑100 (MQ=255)                                                    gTTGCGTCATCTCCTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGGCGCACTTCACTTTGCTGCACACAGGCGCTAAGGCTGAACAGGCTTCCcagcag                                          >  2:999747/1‑100 (MQ=255)                                                     ttGCGTCATCTCCTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGGCGCACTTCACTTTGCTGCACACAGGCGCTAAGGCTGAACAGGCTTCCCAGCAgg                                         <  2:137545/100‑1 (MQ=255)                                                     ttGCGTCATCTCCTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGGCGCACTTCACTTTGCTGCACACAGGCGCTAAGGCTGAACAGGCTTCCCAGCAgg                                         <  1:1163713/100‑1 (MQ=255)                                                      tGCGTCATCTCCTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGGCGCACTTCACTTTGCTGCACACAGGCGCTAAGGCTGAACAGGCTTCCCAGCAGGa                                        >  2:59900/1‑100 (MQ=255)                                                      tGCGTCATCTCCTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGGCGCACTTCACTTTGCTGCACACAGGCGCTAAGGCTGAACAGGCTTCCCAGCAGGa                                        >  2:59977/1‑100 (MQ=255)                                                         gTCATCTCCTGGTTCAGCGTGCTGACGCTGTGTTTCAGCTGGCGCACTTCACTTTGCTGCACACAGGCGCTAAGGCTGAACAGGCTTCCCAGCAGGAGCa                                     <  2:880384/100‑1 (MQ=255)                                                             tctcCTGGTTCAGCGTGCTTACGCTGTGTTTCATCTGGCGCTCTTCACTTTGCTGCACACAGGCGCCAAGGCTGAACAGGCTTCCCAGCAGGAGCAGTCg                                 >  2:1270773/1‑100 (MQ=255)                                                             tctcCTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGGCGCACTTCACTTTGCTGCACACAGGCGCTAAGGCTGAACAGGCTTCCCAGCAGGAGCAGTCg                                 >  2:506747/1‑100 (MQ=255)                                                             tctcCTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGGCGCACTTCACTTTGCTGCACACAGGCGCTAAGGCTGAACAGGCTTCCCAGCAGGAGCAGTCg                                 <  2:1175185/100‑1 (MQ=255)                                                              ctcCTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGGCGCACTTCACTTTGCTGCACACAGGCGCTAAGGCTGAACAGGCTTCCCAGCAGGAGCAGTCGa                                >  2:606767/1‑100 (MQ=255)                                                                 cTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGGCGCACTTCACTTTGCTGCACACAGGCGCTAAGGCTGAACAGGCTTCCCAGCAGGAGCAGTCGaaaa                             <  1:1247313/100‑1 (MQ=255)                                                                   ggTTCAGCGTGCTGACGCTGTGTTTCATCTGGCGCACTTCACTTTGCTGCACACAGGCGCTAAGGCTGAACAGGCTTCCCAGCAGGAGCAGTCGAAAAAc                           >  1:335407/1‑100 (MQ=255)                                                                   ggTTCAGCGTGCTGACGCTGTGTTTCATCTGGCGCACTTCACTTTGCTGCACACAGGCGCTAAGGCTGAACAGGCTTCCCAGCAGGAGCAGTCGAAAAAc                           <  1:275637/100‑1 (MQ=255)                                                                     ttCAGCGTGCTGACGCTGTGTTTCATCTGGCGCACTTCACTTTGCTGCACACAGGCGCTAAGGCTGAACAGGCTTCCCAGCAGGAGCAGTCGAAAAACGt                         >  1:17315/1‑100 (MQ=255)                                                                         gCGTGCTGACGCTGTGTTTCATCTGGCGCACTTCACTTTGCTGCACACAGGCGCTAAGGCTGAACAGGCTTCCCAGCAGGAGCAGTCGAAAAACGTTAGt                     >  1:417692/1‑100 (MQ=255)                                                                            tGCTGACGCTGTGTTTCATCTGGCGCACTTCACTTTGCTGCACACAGGCGCTAAGGCTGAACAGGCTTCCCAGCAGGAGCAGTCGAAAAACGTTAGTAtt                  >  1:203798/1‑100 (MQ=255)                                                                            tGCTGACGCTGTGTTTCATCTGGCGCACTTCACTTTGCTGCACACAGGCGCTAAGGCTGAACAGGCTTCCCAGCAGGAGCAGTCGAAAAACGTTAGTAtt                  <  1:113999/100‑1 (MQ=255)                                                                             gCTGACGCTGTGTTTCATCTGGCGCACTTCACTTTGCTGCACACAGGCGCTAAGGCTGAACAGGCTTCCCAGCAGGAGCAGTCGAAAAACGTTAGTATTc                 >  2:161567/1‑100 (MQ=255)                                                                                  cGCTGTGTTTCATCTGGCGCACTTCACTTTGCTGCACACAGGCGCTAAGGCTGAACAGGCTTCCCAGCAGGAGCAGTCGAAAAACGTTAGTATTCATAgg            <  2:417688/100‑1 (MQ=255)                                                                                     tgtgTTTCATCTGGCGCACTTCACTTTGCTGCACACAGGCGCTAAGGCTGAACAGGCTTCCCAGCAGGAGCAGTCGAAAAACGTTAGTATTCATAGgtcg         >  1:359366/1‑100 (MQ=255)                                                                                        gtTTCATCTGGCGCACTTCACTTTGCTGCACACAGGCGCTAAGGCTGAACAGGCTTCCCAGCAGGAGCAGTCGAAAAACGTTAGTATTCATAGGTCGTcc      >  1:447338/1‑100 (MQ=255)                                                                                           tCATCTGGCGCACTTCACTTTGCGGCACACAGTCTCTAAGCCTGAACAGGCTGCCCAGCAGGATCAGTCGAAAAACGTTAGTATGCATAGTTCGTCCCtt   >  2:760533/1‑100 (MQ=255)                                                                                            cATCTGGCGCACTTCACTTTGCTGCACACAGGCGCTAAGGCTGAACAGGCTTCCCAGCAGGAGCAGTCGAAAAACGTTAGTATTCATAGGTCGTCCCttt  <  1:1228490/100‑1 (MQ=255)                                                                                                  |                                                                                              ATTTTGCATTCAGCCTGTTTTGCTGAGTGATTTTGACTGTTTCTTGATTGAGTTGCGTCATCTCCTGGTTCAGCGTGCTGACGCTGTGTTTCATCTGTCGCACTTCACTTTGCTGCACACAGGCGCTAAGGCTGAACAGGCTTCCCAGCAGGAGCAGTCGAAAAACGTTAGTATTCATAGGTCGTCCCTTT  >  NC_000913/432674‑432864

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: ‑