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	97.9%	131.8 / ‑4.4	47	T108P (ACC→CCC)	ygfZ	folate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (15/31);  total (15/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

CCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTAT  >  NC_000913/3041494‑3041769                                                                                                                                             |                                                                                                                                        cccATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTGTTCTCTAAAGTGccc                                                                                                                                       >  1:54916/1‑143 (MQ=255)        cGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCgcgc                                                                                                                                >  1:1033164/1‑143 (MQ=255)           cGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGACGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAg                                                                                                                             <  1:739208/143‑1 (MQ=255)             ccAAAGGTAAAGTGTGGAGCAATTGACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCTTCAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAgac                                                                                                                           <  1:427298/143‑1 (MQ=255)               aaaGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAgacga                                                                                                                         >  2:707496/1‑143 (MQ=255)                    tAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGAGCCAGACGACGAGc                                                                                                                    >  1:1008878/1‑143 (MQ=255)                        gtgtgGAGCATTTGACGTCTCTTCCGAGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTCAACCGCAGCTGACTGAACCGAACAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCgtgt                                                                                                                <  2:1081806/142‑1 (MQ=255)                         tgtgGAGCAATTTACGTCTGTGCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGtgtg                                                                                                               <  2:1037668/143‑1 (MQ=255)                         tgtgGAGCAATTTACGTCTGTCCCGCGACGGCGATGGCGTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGtgtg                                                                                                               <  2:738879/143‑1 (MQ=255)                           tgGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAANATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGCCGACGAGCGTGtgct                                                                                                             >  1:553522/1‑143 (MQ=255)                            ggAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGtgctg                                                                                                            <  1:1062594/143‑1 (MQ=255)                                  aTTTACGTCGGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGgt                                                                                                      <  2:950402/143‑1 (MQ=255)                                   tttACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGGAAAAATATGCGGTATTCTCTAAGGTGCGCGTTGCGCGAGACCACCGGCGTGGGCTGCTTGgtg                                                                                                     >  1:207368/1‑143 (MQ=255)                                       cGTCGGTTCCGCGACGGCGATCGCTTTGCATGGACTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGc                                                                                                 <  1:921327/143‑1 (MQ=255)                                         tCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTCGGTGTTGCCg                                                                                               >  1:57528/1‑143 (MQ=255)                                               ccGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTc                                                                                         <  1:668440/143‑1 (MQ=255)                                                   caCGGCGATGGCATTGCATGGATTGAACGGCGCAGCGTGCGCGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGgc                                                                                     <  1:284776/142‑1 (MQ=255)                                                        cGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGGGTGCTGCTTGGTGTTGCCGGTTTTCAGgcgcgcg                                                                                >  1:884649/1‑143 (MQ=255)                                                             gCTGTGCATGGATGGAACGGCGCAGCGGGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTAGTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGcgcgccgcg                                                                           <  2:217294/143‑1 (MQ=255)                                                                    aTGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTAAGGCGCGCGCCGCGCTGGCag                                                                    >  1:496899/1‑142 (MQ=255)                                                                           gAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCtt                                                             <  1:838400/143‑1 (MQ=255)                                                                             aCGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTa                                                           >  1:661332/1‑143 (MQ=255)                                                                              cGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAg                                                          <  1:583465/143‑1 (MQ=255)                                                                              cGGCGCAGCGTGCGGCAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAg                                                          <  1:861677/143‑1 (MQ=255)                                                                                  gcAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGac                                                      <  1:898980/143‑2 (MQ=255)                                                                                      cgtgcgtgAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGCCGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGc                                                  <  1:592824/143‑1 (MQ=255)                                                                                       gtgcgtgAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTgag                                                 >  2:656674/1‑141 (MQ=255)                                                                                        tgcgtgAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCt                                                <  2:706017/143‑1 (MQ=255)                                                                                             gAACCGCAGCTGACTGAATTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCgaa                                           <  1:438756/143‑1 (MQ=255)                                                                                              aaCCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCgaaa                                          <  2:1033164/143‑1 (MQ=255)                                                                                              aaCCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCgaaa                                          <  1:51737/143‑1 (MQ=255)                                                                                              aaCCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCGCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCgaaa                                          <  2:1077827/143‑1 (MQ=255)                                                                                              aaCCGCAGCTCACTCAACTGAAAAACTCTGCGGTATTCTCTCAAGCGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTCGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCgaaa                                          <  1:960292/143‑1 (MQ=255)                                                                                                     gcgagctgaacggaaAAAAGATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAAc                                   <  2:328921/138‑1 (MQ=255)                                                                                                      ctgactgaACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACa                                  <  1:102390/143‑1 (MQ=255)                                                                                                          cggaactgcgAAAACATGCGTTATGCTATAAATGGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCTGGTTTTCAGGCGCGCGTCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTa                              <  2:1111865/141‑1 (MQ=255)                                                                                                            gaactgaaAAGATGGGCGGTATGCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCNGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGt                            <  2:553522/143‑1 (MQ=255)                                                                                                             aactgaaAAAATATGCGGTATGCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGAATTCGAAAGAAAAACAGGTAGTc                           <  2:165784/143‑1 (MQ=255)                                                                                                               ctgaaAAAACATGCGGTATTCTCTAACGCGCTCATCTCCCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCaa                         <  2:233443/143‑1 (MQ=255)                                                                                                                tgaaAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGCTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACCGCCTTCGAAAGAAAAAACGGGAGTaaaa                        >  2:695461/1‑143 (MQ=255)                                                                                                                 gaaAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGACAGAAAAACAGGTAGTaaaac                       >  2:300816/1‑142 (MQ=255)                                                                                                                  aaaaaaTATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTGGTCAaagc                      >  1:969219/1‑142 (MQ=255)                                                                                                                     aaaTATGCGGTATTCTCTAAAGTGCCCGTCGCGCCAGACGACGCGCCTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAgg                   <  1:827852/143‑1 (MQ=255)                                                                                                                        gctGCGGTAGACTCTAAAGTGCCCGTCGCGCCAGACGACGAGCGTGTGCTGGTTGGTGGTGCCGGTGTACAGGCGCGCGCCGCGCCGGAAAAACTCTTGAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGgcgc                <  2:256906/141‑1 (MQ=255)                                                                                                                         atGCGGTATTCTCTAAAGTTCTCATCGCGCCAGACGACGAACGTGTGATGCTTGGTGTTGCCGGTTTTAAGGCGCGCACCGCAAAGGCAAATATCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGcgcg               <  2:527421/143‑1 (MQ=255)                                                                                                                                     ctAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTa   >  1:346576/1‑143 (MQ=255)                                                                                                                                     ctAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCGCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTa   <  2:57528/143‑1 (MQ=255)                                                                                                                                      tAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTTGCGAACTGCCTTCGAAAGAAAAACAGGTAATCAAAGAAGGCGCGAACACTATGCGAt  >  1:401801/1‑143 (MQ=255)                                                                                                                                             |                                                                                                                                        CCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTAT  >  NC_000913/3041494‑3041769

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