Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,791,747	A→G	R129R (CGT→CGC)	btuE ←	vitamin B12 transport

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,791,747	0	A	G	98.7%	235.5 / ‑4.7	76	R129R (CGT→CGC)	btuE	vitamin B12 transport
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (55/20);  total (55/21)
Fisher's exact test for biased strand distribution p-value = 2.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

GATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATT  >  NC_000913/1791657‑1791836                                                                                           |                                                                                          gATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCg                                                                                         >  1:1425258/1‑93 (MQ=255)   tGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCggg                                                                                       >  2:4051689/1‑93 (MQ=255)   tGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGCTGCACGGCCTTTGCTGACCATGCggg                                                                                       >  1:4129785/1‑93 (MQ=255)    gACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGc                                                                                      >  2:1515117/1‑93 (MQ=255)    gACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGc                                                                                      >  1:4097628/1‑93 (MQ=255)     aCTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCa                                                                                     >  2:1490205/1‑93 (MQ=255)     aCTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCGTTGCTGACCATGCGGGCa                                                                                     >  1:836774/1‑93 (MQ=255)       ttttCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATa                                                                                   >  2:3711189/1‑93 (MQ=255)        tttCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAg                                                                                  <  2:288326/93‑1 (MQ=255)         ttCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGa                                                                                 <  1:2500126/93‑1 (MQ=255)              tCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCg                                                                            >  1:3690651/1‑93 (MQ=255)               cccTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGc                                                                           <  1:1144691/93‑1 (MQ=255)                 cTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGctc                                                                         >  2:1370158/1‑93 (MQ=255)                 cTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGctc                                                                         >  1:3911065/1‑93 (MQ=255)                 cTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGctc                                                                         >  2:3003979/1‑93 (MQ=255)                  tGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATCGAATCCGctct                                                                        >  2:3618106/1‑93 (MQ=255)                  tGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGctct                                                                        >  1:575751/1‑93 (MQ=255)                     cAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCtt                                                                     >  1:3967770/1‑93 (MQ=255)                     cAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCtt                                                                     >  2:4635072/1‑93 (MQ=255)                      aaCCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTc                                                                    >  1:4379490/1‑93 (MQ=255)                         cAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCgg                                                                 >  2:4085157/1‑93 (MQ=255)                         cAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCgg                                                                 >  2:954991/1‑93 (MQ=255)                           ggAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGgcg                                                               >  1:4199758/1‑93 (MQ=255)                            gAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGgcgc                                                              <  2:784217/93‑1 (MQ=255)                             aaTTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGcgcg                                                             >  2:641741/1‑93 (MQ=255)                              aTTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGa                                                            >  1:4461863/1‑93 (MQ=255)                              aTTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGa                                                            <  1:2020357/93‑1 (MQ=255)                               tttttCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGAc                                                           >  1:2293763/1‑93 (MQ=255)                                ttttCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACt                                                          <  2:1617048/93‑1 (MQ=255)                                 tttCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTg                                                         <  2:1073815/93‑1 (MQ=255)                                 tttCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTg                                                         >  2:3236695/1‑93 (MQ=255)                                  ttCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGc                                                        <  2:1206386/93‑1 (MQ=255)                                  ttCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGc                                                        <  2:525400/93‑1 (MQ=255)                                  ttCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGc                                                        <  1:734293/93‑1 (MQ=255)                                    cAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCgg                                                      >  1:2940318/1‑93 (MQ=255)                                     aaaaTTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGt                                                     <  1:1624935/93‑1 (MQ=255)                                     aaaaTTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGt                                                     >  2:740887/1‑93 (MQ=255)                                     aaaaTTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGt                                                     >  2:699892/1‑93 (MQ=255)                                      aaaTTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTc                                                    >  2:1127044/1‑93 (MQ=255)                                      aaaTTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTc                                                    >  2:1940027/1‑93 (MQ=255)                                       aaTTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCg                                                   <  1:2993248/93‑1 (MQ=255)                                        aTTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCgg                                                  >  1:2465170/1‑93 (MQ=255)                                            cATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGgcgcg                                              >  1:463128/1‑92 (MQ=255)                                            cATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCt                                              <  2:4476887/93‑1 (MQ=255)                                             aTAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGgcgccg                                             >  1:2985715/1‑91 (MQ=255)                                              tAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGc                                            >  1:1228997/1‑93 (MQ=255)                                                aaaTATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCgg                                          >  2:780160/1‑93 (MQ=255)                                                  atatCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGca                                        >  1:3718415/1‑93 (MQ=255)                                                  atatCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGca                                        >  2:2672878/1‑93 (MQ=255)                                                   tatCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGcaa                                       >  2:2692442/1‑93 (MQ=255)                                                   tatCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGcaa                                       >  1:2480481/1‑93 (MQ=255)                                                    atcatcCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGcaat                                      <  2:3201561/93‑1 (MQ=255)                                                     tcatcCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGcaatc                                     >  2:950424/1‑93 (MQ=255)                                                       atcCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGcaatcaa                                   >  2:2008324/1‑93 (MQ=255)                                                       atcCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGcaatcaa                                   >  2:2316927/1‑93 (MQ=255)                                                        tcCGGGTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGcaatcaat                                  <  1:3389257/93‑1 (MQ=255)                                                            ggTACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAAttttt                              <  1:690816/93‑1 (MQ=255)                                                              tACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGa                            >  1:4497619/1‑93 (MQ=255)                                                              tACAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAACCAATTTTTGa                            >  2:3099548/1‑93 (MQ=255)                                                               aCAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGAt                           >  2:1504189/1‑93 (MQ=255)                                                               aCAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGAt                           >  2:38533/1‑93 (MQ=255)                                                                cAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATa                          >  1:1986297/1‑93 (MQ=255)                                                                cAGCGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATa                          >  1:2820311/1‑93 (MQ=255)                                                                   cGGTGCACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAg                       <  1:941073/93‑1 (MQ=255)                                                                        cACGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGAt                  <  1:3341194/93‑1 (MQ=255)                                                                         aCGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATg                 >  1:37461/1‑93 (MQ=255)                                                                         aCGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATg                 <  2:1288025/93‑1 (MQ=255)                                                                         aCGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATg                 >  1:1934759/1‑93 (MQ=255)                                                                          cGGCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATgg                >  1:4103974/1‑93 (MQ=255)                                                                            gCCTTTGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCg              >  2:180898/1‑93 (MQ=255)                                                                              cTTTGCTGGCCATGCGGGCATAGAATCCGCTCTCTGCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTc            <  2:3043433/93‑1 (MQ=255)                                                                               tttGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTcc           >  2:1421924/1‑93 (MQ=255)                                                                                 tGCTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCtt         >  2:4245234/1‑93 (MQ=255)                                                                                   cTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCg       <  1:3934919/93‑1 (MQ=255)                                                                                   cTGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCg       >  2:3054928/1‑93 (MQ=255)                                                                                    tGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGc      >  1:2678286/1‑93 (MQ=255)                                                                                    tGACCATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGc      >  2:3134793/1‑93 (MQ=255)                                                                                        cATGCGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCAtt  <  2:514819/93‑1 (MQ=255)                                                                                           |                                                                                          GATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATT  >  NC_000913/1791657‑1791836

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