Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,041,398	T→G	L29R (CTT→CGT)	ygfZ →	folate‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,041,398	0	T	G	98.3%	179.1 / NA	58	L29R (CTT→CGT)	ygfZ	folate‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (26/31);  minor base A (1/0);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTCTTTTACCATTAGTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTG  >  NC_000913/3041261‑3041537                                                                                                                                          |                                                                                                                                            gCTCTTTTACCATTAGTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCCTTTACCGCTGCCGCTGATGACGCTCGATGCCTGGGCGCGTGcca                                                                                                                                        >  2:776425/1‑143 (MQ=255) gCTCGTTTACCATTAGTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGcca                                                                                                                                        <  1:988800/143‑1 (MQ=255)        tACCAGTAGTCATTAATACGCCGCTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTAc                                                                                                                                 <  2:996881/143‑1 (MQ=255)             ttAGTCATTAATACGCCGTTAAGCATCTCAGGCCTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTCCGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGgcg                                                                                                                            <  2:688542/143‑1 (MQ=255)              tAGTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGgcgc                                                                                                                           <  1:1053908/143‑1 (MQ=255)              tAGTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGgcgc                                                                                                                           <  1:483593/143‑1 (MQ=255)                gTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGCGCGCGTGCCACCCCTGCTGGCGCgg                                                                                                                         >  1:1076747/1‑143 (MQ=255)                gTCATTAATACGCCGTTAAGCAACTCAGGACCCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGCCGCGCGCCACCATTACTGGCGCgg                                                                                                                         >  2:785373/1‑143 (MQ=255)                 tCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGa                                                                                                                        <  2:276362/143‑1 (MQ=255)                 tCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGa                                                                                                                        <  1:390442/143‑1 (MQ=255)                 tCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTATTACACCTTTTCCTCCCCGTAAGCCTAAGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCATGCCATCGTTACTGTCGCGGa                                                                                                                        >  1:709906/1‑143 (MQ=255)                   aTTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACa                                                                                                                      <  2:1023886/143‑1 (MQ=255)                    ttAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTCGGGCTTCTGTCCGTTTACCGCTGACGCTGATGACGCTTGCTGACTGGGCCCGTGCCACCACTACTGGCGCGGAccg                                                                                                                     >  2:691842/1‑141 (MQ=255)                    ttAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCCGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAg                                                                                                                     <  1:596471/143‑1 (MQ=255)                          cGCCGTTACGCACCACAGGCCTCTACTATGGCTTTTACACCTTTTCCCCCCCGTCAGCCTCCGGCTTCTGCCCGTTTACCGCTGACGCTGATGTCGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGaaaa                                                                                                               <  2:633773/143‑1 (MQ=255)                          cGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGaaaa                                                                                                               <  1:728117/143‑1 (MQ=255)                            ccGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAat                                                                                                             >  1:435821/1‑143 (MQ=255)                              ggtAACCAACTCAGGCCTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAatat                                                                                                           <  2:973623/141‑1 (MQ=255)                              gTTAAGCAACTCAGTTCTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAatat                                                                                                           <  1:991920/143‑1 (MQ=255)                                    cAACTCTGGACTTTATTATGGCTTGTACACCTTGTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAg                                                                                                     <  2:522462/143‑1 (MQ=255)                                    cAACTCAGCACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAg                                                                                                     <  1:552049/143‑1 (MQ=255)                                     aaCTCAGGGCTCTATTGTGGCTTCTACACCGTTGCCGCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAgg                                                                                                    <  1:398522/143‑1 (MQ=255)                                            gACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGt                                                                                             <  1:767221/143‑1 (MQ=255)                                              cgctATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGa                                                                                           <  1:1047684/141‑1 (MQ=255)                                               tctATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTGCTGGCGCGGACAGCGAAAACTATATGCAGGGTCAGGTGag                                                                                          >  1:557746/1‑142 (MQ=255)                                                  attatGGCTTTTACTCCTTTCGCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGAcagc                                                                                       <  2:180118/143‑1 (MQ=255)                                                    tatGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGAcagcag                                                                                     <  2:435821/143‑1 (MQ=255)                                                     atGGCTTTTACACCTTGTCCACCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGTCAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGa                                                                                    <  1:743880/143‑1 (MQ=255)                                                         cTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGCGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCGGGTGACAGCAGATggc                                                                                >  2:173821/1‑141 (MQ=255)                                                                   ttttCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATAGATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATgg                                                                      >  2:19290/1‑143 (MQ=255)                                                                     ttCCTCCCCGTCAGGCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCa                                                                    <  1:366656/143‑1 (MQ=255)                                                                      tcctccCCGTCAGCCTACGGCTTCTGCCCGTCTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCag                                                                   <  1:281868/143‑1 (MQ=255)                                                                        ctccCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCagaa                                                                 <  1:744957/143‑1 (MQ=255)                                                                          ccccGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGCTGACTGGGCGCGTGCCACCCGTACTGGCGCGGACAGCGAAACAGAGACGCAGGGTCCGGGGACAGCAGATGGCAGCCCGATGGCaggcgc                                                               >  2:312760/1‑139 (MQ=255)                                                                           cccGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAACATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGAt                                                              >  1:119333/1‑143 (MQ=255)                                                                            ccGGCAGCCTACGGCTGCTGCCCGTCTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCAGCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATc                                                             <  2:95949/143‑1 (MQ=255)                                                                             cGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCa                                                            <  1:23386/143‑1 (MQ=255)                                                                                    cTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTg                                                     >  1:650009/1‑143 (MQ=255)                                                                                    cTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCCCGTGCCACCATTACTGGCGCGGACAGCGAAAAAAATATGCCGGCTCACGTGACAGCCCACGTCAGCCAGATGCCACAAGACCAGCACCtc                                                     >  2:1106985/1‑142 (MQ=255)                                                                                        ggCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCCGAAGATCAGCACCTGCTcg                                                 >  1:604243/1‑143 (MQ=255)                                                                                               gCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAACTATATGCCGCGGCAGGTGACCGCAGATGGCAGCCAGATGGCAAGAGATCAGCACCTGCCCGCCGcccc                                          >  2:545888/1‑142 (MQ=255)                                                                                                cccGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCCGCACCTGCCCGcccccccg                                         >  2:831406/1‑141 (MQ=255)                                                                                                 ccGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCAtt                                        <  1:592396/143‑1 (MQ=255)                                                                                                  cGTTTGCCGCTGAGGCTGGTGATGCTTGAGGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAGGTATTTGCAGGGTCCGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTg                                       <  1:312760/143‑1 (MQ=255)                                                                                                  cGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTg                                       >  2:144439/1‑143 (MQ=255)                                                                                                         cgctgacgctgaTGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATCTATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTACTCGCCGCCCCTTGCgccgcc                                >  2:650743/1‑143 (MQ=255)                                                                                                         cgctgacgctgaTGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCCGGGTCAGGTGACAGCAGACGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCAATGCGACGcc                                >  2:533315/1‑143 (MQ=255)                                                                                                                   gatgcCGTTTGGTGACCGGGCGCGTGCTGCCATTACTGGGGCGCACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTTGCCGCCCATTGCGACGCCAAAGGTaaaa                      <  2:887089/143‑1 (MQ=255)                                                                                                                       actcTTGATGACTGGGCGCGTGCCACCATTACCGGCGCGGACAGCGAAAAATATATGCAGCGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCCGCTCGCCGACCATTGCGACGCCAAAGGTAAAAtgtg                  <  1:137576/140‑1 (MQ=255)                                                                                                                         gCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAGAGGTCCGATGTGGa                >  1:322842/1‑143 (MQ=255)                                                                                                                         gCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGa                >  1:263572/1‑143 (MQ=255)                                                                                                                           ttGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACCGCGAAAAATATATGCATGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAACGTGGAGc              >  1:30591/1‑143 (MQ=255)                                                                                                                               tgaCTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGCGGAGCaagt          >  1:411129/1‑141 (MQ=255)                                                                                                                                gaCTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAACGGTAAAATGTGGAGCAAttt         >  1:975972/1‑143 (MQ=255)                                                                                                                                gaCTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAAGGTGGAGCAAttt         >  1:906974/1‑143 (MQ=255)                                                                                                                                   tGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGTTGTAGAGCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATATGGAGCAATTTACg      >  1:835927/1‑143 (MQ=255)                                                                                                                                    gggCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTACAATGTGGAGCAATTTACGt     >  2:1097498/1‑143 (MQ=255)                                                                                                                                      gcgcGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCATGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCt   <  2:70055/143‑1 (MQ=255)                                                                                                                                       cgcgTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCCTAGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTg  >  1:1081999/1‑143 (MQ=255)                                                                                                                                          |                                                                                                                                            GCTCTTTTACCATTAGTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGATGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTG  >  NC_000913/3041261‑3041537

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