Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,041,634	A→C	T108P (ACC→CCC)	ygfZ →	folate‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,041,634	0	A	C	98.2%	161.8 / NA	55	T108P (ACC→CCC)	ygfZ	folate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (22/31);  total (24/31)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

TGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGGT  >  NC_000913/3041499‑3041772                                                                                                                                        |                                                                                                                                           tGCGGCGCCAAATGTAAAAGGTGGAGCAATTTACGTCTGTTCCGCGGCGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATcgc                                                                                                                                     <  2:661885/143‑1 (MQ=255)   cgcggCCAAAGGTAACATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCgcgc                                                                                                                                   <  2:2014/139‑1 (MQ=255)       ggctgtggTAAAATGTTGAGTAAGTGACGTCTGTTCCGCGACGGCGATGGTTTTGCATGGATTGAACGGTGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAga                                                                                                                               <  2:868692/137‑1 (MQ=255)        ccAAAGGTAAAATGTGGAGCAATTTACATCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGACAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAgac                                                                                                                              >  2:532538/1‑143 (MQ=255)          aaaGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAgacga                                                                                                                            <  2:296404/143‑1 (MQ=255)          aaaGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAgacga                                                                                                                            >  1:296404/1‑143 (MQ=255)                   aTGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGGTTGAACGGCGCCGCGTGCGTGAACCACAGCTGACTGCACAGAAATAATATGCGGTATTCCCTAAAGTGCCCATCCAGCCAAACGACGAGCgtgt                                                                                                                   >  2:626307/1‑143 (MQ=255)                         ctccaTTTCCGGCTGTTGCGAGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTgctgct                                                                                                             <  1:136587/139‑1 (MQ=255)                         aGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTgctgct                                                                                                             >  1:611703/1‑143 (MQ=255)                         aGCAATTTACGTCTGTGCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTgctgct                                                                                                             <  2:406205/143‑1 (MQ=255)                          gCAATGTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCtt                                                                                                            <  1:205566/143‑1 (MQ=255)                           cAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTg                                                                                                           <  1:641295/143‑1 (MQ=255)                                    gcTGTTCCGCGGCGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACGGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCg                                                                                                  <  1:565894/142‑1 (MQ=255)                                       gtgccGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGTGGTATTTTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGtt                                                                                               <  1:673692/140‑1 (MQ=255)                                       gTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGtt                                                                                               >  1:552417/1‑143 (MQ=255)                                        ttCCGCTACGGCTATGGCTTTGTATGGAGTGAACGGCGCTGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGttt                                                                                              <  1:983194/143‑1 (MQ=255)                                             cgACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAgg                                                                                         >  2:274428/1‑143 (MQ=255)                                                    gATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATCTGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGgcgcgcgc                                                                                  <  1:762594/143‑1 (MQ=255)                                                           ttGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAACGTGGCCATCGCGCAAGGCGACGAGCGCGGGCTGCTTGTTGTCGCCGGTTTTCCGGCTCGCTCCGCGGTg                                                                           >  1:40233/1‑143 (MQ=255)                                                           ttGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTg                                                                           <  1:218097/143‑1 (MQ=255)                                                               aTGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTATAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGCTTTTCAGGCGCGCGCCGCCCTGGcca                                                                       >  1:181084/1‑141 (MQ=255)                                                                 ggATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAGAt                                                                     >  1:738651/1‑143 (MQ=255)                                                                   gttGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAAtct                                                                   <  2:149087/142‑1 (MQ=255)                                                                   gttGAACGGCGCAGCGTGCGGGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAAtct                                                                   <  2:150582/142‑1 (MQ=255)                                                                   agtGAGCGGCGCAGGGGGCGTGAACCGCAGCTGACTGAACGGAAAAAATATGCGGTATTCTCTAAAGTGCCCGTCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAAtct                                                                   <  2:12712/141‑1 (MQ=255)                                                                   aTTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAAtct                                                                   <  2:949907/143‑1 (MQ=255)                                                                    ttGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACCGAAAAAATATGCGGTATTCTATAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGCGTTCAGGCGCGCGCCGCGCTGGCcaacccc                                                                  >  1:456393/1‑139 (MQ=255)                                                                          ggcggcgcgtgcgtGAGCCGCAGCTGGCTGCACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGc                                                            <  2:880020/137‑1 (MQ=255)                                                                                   tgtgtgAACCGCAGCTGACTGAACTGAAAAGATATGCGGTATTCTCTAAAGTGCCTTTCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCt                                                   <  1:1026099/140‑1 (MQ=255)                                                                                    gcgtgCACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCtt                                                  <  2:841904/143‑1 (MQ=255)                                                                                       tgAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCga                                               >  2:912523/1‑143 (MQ=255)                                                                                           ccGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTGgaaaga                                           >  1:986070/1‑143 (MQ=255)                                                                                               aGCAGACTGAACTGAAAAAATATGCGGTATTCGCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGaaaaa                                       <  1:484232/143‑1 (MQ=255)                                                                                                 ctgtctgaACTGAAAAAATATGCGGTATTCTCGAAAGGGCCCAGCGCGCCAGACGACGAGCGTGTGCTGCTGGGTGGTGCCGGTTTTCAGGCGCGCGCCGCGCTGGTAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACa                                     <  1:691474/143‑1 (MQ=255)                                                                                                 ctgactgaACTGAAAAAATATGTGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACa                                     <  1:94527/143‑1 (MQ=255)                                                                                                      tgaactgaaAAAATATGCGGTATTCTCTAAAGTGACCATCGCACGAGACGACGAGCGTGTCCAGCTTAGAGTAGACGGATTTCAGGCGAGCGAAGCACAGCCAAAACACTGTAGCGAAATCACATCGAACAAAAAACAGATAg                                >  2:118920/1‑143 (MQ=255)                                                                                                       gaactgaaAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGt                               <  1:469092/143‑1 (MQ=255)                                                                                                         actgaaAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTTGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCCAATCTCTTTAGCGACCTGCCTTCGGAAGAAAAACAGGGAGTcc                             >  1:478430/1‑142 (MQ=255)                                                                                                           tgaaAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACACGTAGTCaaa                           >  1:477789/1‑143 (MQ=255)                                                                                                           tgaaAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAAAAATTCGAACCTACTGCCTTCGAAAGAAAAACAGGTAGGCaaa                           >  2:220150/1‑143 (MQ=255)                                                                                                           tgaaAAAATATGCGGTATTCTCTAAAGTGCCAATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTGGCCGGATTTGAGGCGCGCGCCGCGCTGGCAAATCTCTTTCGCGAACTGCCTTGGAAAGGAAGACAGGTAGtaaca                           >  2:314684/1‑139 (MQ=255)                                                                                                             aaaaatAATGCGTTATTCGTTAACGGGCCCATCGCGTCAGACGCGCAGCGTGTGTTGCGTGGTGTTGCCGGTAGTCAGGCGCGCGCCGCGCTGGCAAATCGCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAaaga                         <  2:109508/143‑1 (MQ=255)                                                                                                             aaaaaaTGTGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGCCGCGCGCCGCGCTGGCAAATATCTTTGGCGAACGGCCGTCGAAAAAAAAACAGCTAGTACaaga                         >  2:954931/1‑143 (MQ=255)                                                                                                             aaaaaaTATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGCGCGTGTGCTGCTTGGTGTTGCCGGTTTTAAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGAAGTCAacga                         >  2:520488/1‑143 (MQ=255)                                                                                                             aaaaaaTATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTCGCGTACTTGCTGCGAAAGAAAAACAGGTAGTCAaaga                         >  1:287120/1‑143 (MQ=255)                                                                                                              aaaaaGAGGCGGTCTTCTCTAAAGTGCCCAGCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTCTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAaagaa                        <  2:787876/143‑1 (MQ=255)                                                                                                                    atGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGcgcg                  <  2:827386/143‑1 (MQ=255)                                                                                                                     tGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAACACCAGGTCGTCAAAGAAGgcgcca                 >  2:1060074/1‑141 (MQ=255)                                                                                                                     cgCGGTAGTGTCTCAAGTGCCCACCGCGCCAGACGACGAGCGTGTGCTGCTTAGTGTTGCCGGTTTTCAGGCGCGCGCCGAGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGa                 <  2:466710/142‑1 (MQ=255)                                                                                                                      gCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCGCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGAc                <  2:717380/143‑1 (MQ=255)                                                                                                                            ttCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGG‑TTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACGGGTAGTCCAAGAAGGCGCGACCACGTTg         >  1:668/1‑143 (MQ=255)                                                                                                                              ctcgAAAGTGCCCATCGCGCGAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTACCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGc        <  2:711561/143‑1 (MQ=255)                                                                                                                               tctAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTggt       >  1:429309/1‑141 (MQ=255)                                                                                                                                 tAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTAt     <  1:207186/143‑1 (MQ=255)                                                                                                                                    aGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGAAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGGCAAAGAAGGCGCGACCACGTTGCTATGGt  >  1:974762/1‑143 (MQ=255)                                                                                                                                        |                                                                                                                                           TGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGGT  >  NC_000913/3041499‑3041772

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