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%	151.1 / NA	57	T108P (ACC→CCC)	ygfZ	folate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (27/29);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGG  >  NC_000913/3041495‑3041771                                                                                                                                            |                                                                                                                                          ccAGTGCGACGCCAAAGGTAAAATGGGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCa                                                                                                                                        <  1:974444/143‑1 (MQ=255)  cATTGCGGCGCCAAAGGTAAAATGTGGAGCAATGTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCAt                                                                                                                                       <  1:172062/143‑1 (MQ=255)       cGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCgcgc                                                                                                                                  >  1:296012/1‑143 (MQ=255)        gACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTGGCAAGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGGATTCTCGAAAGTGACCATCGCGcc                                                                                                                                 >  2:740750/1‑143 (MQ=255)           gCCAAAGGTAAAATGCGGAGCAATTTACGTCTGTTCCGCGACTGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTTAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAga                                                                                                                              <  1:274349/143‑1 (MQ=255)            ccAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAgac                                                                                                                             <  1:816858/143‑1 (MQ=255)             cAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAgacg                                                                                                                            <  1:763568/143‑1 (MQ=255)              aaaGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTGAAGTGCCCATCGCGCCAgaggg                                                                                                                           >  2:312980/1‑140 (MQ=255)                 ggTAAAATGGGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGCACTGAAAAAATATGCGGTTTTCTCTAAAGTGCCCATCGCGCCAGAcgacga                                                                                                                        <  2:524735/143‑1 (MQ=255)                    aaaaTGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCg                                                                                                                     >  1:630407/1‑143 (MQ=255)                           gcaGCAATTTACGTCTGTTCCGCGACTGCGCTGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGtgctg                                                                                                              <  2:789580/141‑1 (MQ=255)                             aGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTgctgct                                                                                                            >  1:129143/1‑143 (MQ=255)                              gCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGTAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCtt                                                                                                           >  2:491622/1‑143 (MQ=255)                               caatgcaCGTCGGTCCCGCGACGGCGATGGCTTGGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTg                                                                                                          <  1:447521/137‑1 (MQ=255)                                  tttACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTCAAGTGCCCCTCGCGCCAGACGAGGGCCGTGGGCTGCTTggcg                                                                                                       >  2:724136/1‑141 (MQ=255)                                       gTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTACTCTCTAAAGTGCCCATAGCGCCAGACGACGAGCGCGTGCTGCTTGGGGTTGcc                                                                                                  >  2:753377/1‑143 (MQ=255)                                        tCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGCAAAAATATACGGTATTCTCTAAAGTCCCCATCCCGCCAGACGACGAGCGTTTGCTGCTTGGTGTTCCCg                                                                                                 >  2:765731/1‑143 (MQ=255)                                           gTTCCCCGGCGGGGATGGCTTGGCATGGAAATCACGTCTCAGTGTGCGTGAACCGCAGCTGACGGAACTGAAAAGATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGtt                                                                                              <  2:437816/143‑1 (MQ=255)                                              ccGCGACGGCGCTGGTCTTGGATGGATTGAACGGCGTAGGGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTc                                                                                           <  2:113214/143‑1 (MQ=255)                                              ccGCGACGGCGATGGCTTTGCATGGATTGAATGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGCGTTGCCGGTTTTc                                                                                           >  1:814378/1‑143 (MQ=255)                                                              tttGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCt                                                                           >  1:494276/1‑143 (MQ=255)                                                                tGCCTGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTgg                                                                         <  2:695943/143‑1 (MQ=255)                                                                   aTGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCaa                                                                      >  1:487857/1‑143 (MQ=255)                                                                       aTTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAAtct                                                                  >  1:406401/1‑143 (MQ=255)                                                                          gAACGGCGCAGCGTGCGTGAACCGCAGCTGACAGAACTGAAAAAATATGCGGTATTCTCTACAGTGCCCATCGCGCCAGACGAAGCGCGCATGCTGCTTGGCGTTGCCGGTTTTCAGGCGCGCGCCGCGCAGGCGCATCTCtt                                                               >  2:349491/1‑143 (MQ=255)                                                                           aaCGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCCGGCGCGCGCCCCGCCGGCAAATccctct                                                              >  2:87436/1‑141 (MQ=255)                                                                           aaCGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAACTCTCttt                                                              >  2:519258/1‑143 (MQ=255)                                                                           aaCGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCttt                                                              >  2:142125/1‑143 (MQ=255)                                                                           aaCGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATGCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGCTGACGGTTTTCACGCGCGCGCCGCGCTGGCATATctctcc                                                              >  2:699598/1‑141 (MQ=255)                                                                             cGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGAAACTTTCAGTTGTTCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAg                                                            <  1:991746/143‑1 (MQ=255)                                                                                cgcAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCAATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGGATTCAGGCGCGGGCCGCGCTGGCAAATCTCTTTAGCGa                                                         >  2:861839/1‑143 (MQ=255)                                                                                 gcAGCGTGCTTCAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGaa                                                        <  1:830194/143‑1 (MQ=255)                                                                                   aGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAAGATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACt                                                      <  1:962791/143‑1 (MQ=255)                                                                                          gtgAACCGCAGCTGACTGAACTGACAAAATATGCGGTATTCTCTACGGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTTGTGTTGCCGGTTTTCAGGCGCTCGTCGCGCTGGCAACTCTCTTTAGCGAACTGCCTTCg                                               <  2:767888/143‑1 (MQ=255)                                                                                          gtgAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCg                                               <  2:659707/143‑1 (MQ=255)                                                                                               ccGCAGCTGACTGGACTGAAAAACTATGCGGTGTTCTCTAAAGTTCCCATCGGGGCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCgaaaga                                          <  2:706769/143‑1 (MQ=255)                                                                                               ccGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGGGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTCGGGATCTGCGTGCgagcga                                          >  1:440908/1‑139 (MQ=255)                                                                                                ggCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGGGCCCATCGCGCCAGACGACGAGCGTGCGCTGCTTGGTGTTGCCGGTCTTCAGGCGCGCGCCGCGCCCGGAAATCGCTGTAGCGAACTGGCATCggaagaa                                         >  2:360838/2‑143 (MQ=255)                                                                                                cGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGTTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCgaaagaa                                         <  2:296012/143‑1 (MQ=255)                                                                                                  cAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGaaaa                                       >  1:560620/1‑143 (MQ=255)                                                                                                  cAGCTGACTGAACTGAAAAAATATGCGGTATTATCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCgaaagaaga                                       >  1:57221/1‑141 (MQ=255)                                                                                                   agttgactgaACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGaaaaa                                      <  1:589622/140‑1 (MQ=255)                                                                                                      tgactgaACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAg                                   <  1:867982/143‑1 (MQ=255)                                                                                                             actgaaAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTGGTCa                            >  1:970468/1‑143 (MQ=255)                                                                                                               tgaaAAAATATCGGGTATTCTCTAAGGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGAGCGAGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCaaa                          <  2:812245/143‑1 (MQ=255)                                                                                                                  gaaaaTATGCGGTATGCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAaagaa                       <  1:312980/142‑1 (MQ=255)                                                                                                                  aaaaaTATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGACAGAAAGACAGGTAGTCAaagaa                       >  1:900645/1‑143 (MQ=255)                                                                                                                    aaaGATGCGGTATGCTCTAAAGTGCCCATCGCGCCAGACGACGGGCGTGTGCTGCTTGGTGTTGCGGGTTTCCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAgg                     <  2:776915/143‑1 (MQ=255)                                                                                                                         tGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAACGAACACCCGGTCGCCAAAGAACGCGCGa                >  1:40128/1‑143 (MQ=255)                                                                                                                          gCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGAc               <  1:60129/143‑1 (MQ=255)                                                                                                                             ctatccgctAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTGCAGGCGCGCGCCGCGCTGGCAAATCTCTTCAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCAc            <  2:576777/136‑1 (MQ=255)                                                                                                                                 tcgctAAGGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCCACGCGCTGGCAAATCTCCTTAGCGAACGGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTg        <  2:603135/140‑1 (MQ=255)                                                                                                                                   tctAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTCTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCCAATCTCTTTAGCGCACTGCCTGCGAAAGAAAAACAGGCAGTCCATGAAGGCGCGACCACTTCGcc      >  2:279550/1‑142 (MQ=255)                                                                                                                                     tAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGCAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTAt    >  2:78354/1‑143 (MQ=255)                                                                                                                                      caaGTGCCCATCGCGTCAGACGACGAGCGTGGGCTGCTTGGTGTGGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATg   <  1:457027/142‑1 (MQ=255)                                                                                                                                      aaaGTGCCCATCGCGCCCGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCCTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATg   <  2:401265/143‑1 (MQ=255)                                                                                                                                       aaGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATgg  <  1:399515/143‑1 (MQ=255)                                                                                                                                            |                                                                                                                                          CCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGG  >  NC_000913/3041495‑3041771

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