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.8%	137.6 / NA	46	T108P (ACC→CCC)	ygfZ	folate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (25/20);  total (26/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGG  >  NC_000913/3041494‑3041771                                                                                                                                             |                                                                                                                                          cccATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGccc                                                                                                                                         <  1:766049/143‑1 (MQ=255)         gACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGCGAACCGCAGCTGACTGAACTGAAAAAATACGCGGTATTCTCTAACGTCCCCATCGCGcc                                                                                                                                 >  2:1121938/1‑143 (MQ=255)               aaaGGTAAAATGTGGAGCAAGTTACGTCTGTTCCGCGACGGCGATGGCTTTGCAGGGATTGAACGGCGCAGCGTGCGTGAACCGAAGCTGACTGAACTGAAAAAGGGTGCGGGATGCTCTAAAGAGCCCAGAGCGCGAggcga                                                                                                                           >  1:431785/1‑143 (MQ=255)                   gTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGacgacgcg                                                                                                                       >  1:753457/1‑141 (MQ=255)                   gTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCCCTAAAGTGCCCATCGCGCCACACGACGAg                                                                                                                       >  2:553188/1‑143 (MQ=255)                    tAAAATGTGGAGCAATTTACGTCTGTTCCTCTACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGc                                                                                                                      <  2:48253/143‑1 (MQ=255)                        aTGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCGTTGCAGGGACTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCgtgt                                                                                                                  <  1:3147/143‑1 (MQ=255)                         agtgGAGCAATTTACATCTGTTCCGCGACGGCGATGGCTCTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGtgtg                                                                                                                 <  1:702927/142‑1 (MQ=255)                             gAGCAATTTACGTCTGTTCCGCGACGGCAATGGATATACATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGCACTGCAAAAATATGCGGGAGTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGtgcttc                                                                                                             >  2:555789/1‑141 (MQ=255)                               gCAACTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCtt                                                                                                           <  1:637348/143‑1 (MQ=255)                                gaagtgaCGTCTGTTCCGCGACGGCGGTGTCTTTGCATGGATTGCGCCGGGGCCCATTCCGGCACAGCAGTGGTTGGCACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTg                                                                                                          <  2:806212/137‑1 (MQ=255)                                cactttCCGTCTGTGCCGCGACGGCGATGGCTTGGCATGCATTCAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCAGCGCGCCAGACGACGAGCGTGTGCTGCTTg                                                                                                          <  1:473375/140‑1 (MQ=255)                                cAATTTACGTCTGTGCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTg                                                                                                          <  1:639215/143‑1 (MQ=255)                                 aaTTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTgg                                                                                                         >  1:217881/1‑143 (MQ=255)                                        gTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGcc                                                                                                  <  1:307163/143‑1 (MQ=255)                                        gTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGcc                                                                                                  >  1:791050/1‑143 (MQ=255)                                         tCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCg                                                                                                 <  1:1030406/143‑1 (MQ=255)                                             ttCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAGAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCCTGGGGTTGCCGGgtt                                                                                             >  2:565164/1‑143 (MQ=255)                                             ttCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGttt                                                                                             >  1:920709/1‑143 (MQ=255)                                                cgcgACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGCCTGAACTGAAAAAATATGCGGTCTTCTCTAAAGTGCCCATCGCGCCAGACGACGAACGCGTGCTGCGTGCTGTTGCCGGTTCTCa                                                                                          >  2:210084/1‑143 (MQ=255)                                                 gcgcCGGCGATGGCTTTGCATGGATTGAGCGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAGAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAg                                                                                         <  2:204858/143‑1 (MQ=255)                                                 gcgACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCAAGGACCTAAAGAACTATTGTGTGCAGCCTGGTGTTGCCGGTTTTCAg                                                                                         >  2:234825/1‑143 (MQ=255)                                                            ggCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGcgcgccgc                                                                              >  1:1089922/1‑143 (MQ=255)                                                                  gCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCGCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGGGTTGCCGATTTTCCGGCGCGCGCCGCGCTGGc                                                                        >  2:763992/1‑143 (MQ=255)                                                                       gATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGAGGTATGCTCTAAAGTGCCCATCGCGCCAGACGACGCGCGTGTGCTGCTTGGTGTTGCCGCTTTTCAGGGGCGCGCCGCGCTGGCagctc                                                                   >  1:614931/1‑139 (MQ=255)                                                                        aTTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCCCGCGCCGCGCTGGCAAAtct                                                                  >  1:422749/1‑143 (MQ=255)                                                                              cggcgcggcgGGCGTGAACCGCTGCCGGCTGAACTGAAAAAATGTGCGGTATTCTCTAAAGTGCCCATCGCGCCAGGCGACGAGCGGGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAg                                                            <  2:430887/143‑1 (MQ=255)                                                                              cGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAg                                                            >  2:430428/1‑143 (MQ=255)                                                                               gtcgtatcgtgcgtgATCCTCAGTTTTCTTTTCTGATTAAATTTGCGGTTTTCTCTAAAGTGCCCATCGAGCCAGACGACGAGCGTGGGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGc                                                           <  1:137297/136‑1 (MQ=255)                                                                                  gccgcgtgcgtGAACCGCAGCTGACTGAACTGACAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGaa                                                        <  1:230288/140‑1 (MQ=255)                                                                                      cgtgcgtgAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCa                                                                                          >  1:96998/1‑105 (MQ=255)                                                                                      cgtgcgtgAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCa                                                                                          <  2:96998/105‑1 (MQ=255)                                                                                         gcgtgAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCCTGTGCTGCTGGGTGTTGCCGGTTCTCAGGCGCGCGCCGCGCTCGGAAACGTCTTTGGCGGACCGccgt                                                 >  2:414181/1‑141 (MQ=255)                                                                                              aaCCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCCGGCAAATCTCTTTAGGGCACGGCCGTCgaaa                                            >  2:818469/1‑143 (MQ=255)                                                                                              aaCCGCAGCTGACTGAACTGAAAAAATATGCGGTATGCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCgaaa                                            <  1:386715/143‑1 (MQ=255)                                                                                                      ctgactgaACTGAAAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGGGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAACTCTCTTTAGCGGACTGCCTGCgcaagaaacgca                                    >  1:1086460/1‑139 (MQ=255)                                                                                                       gtgctgaactgaaAAATGATGCGGTATTCTCTAAAGAGCCCCTCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTGGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCGTCGAAAGAAAAACAg                                   <  2:276575/140‑1 (MQ=255)                                                                                                           gagatagaaAAAATACTCGGTATTCTCGAAAGGGCGAGTCTCGCGAGACGACGAGCGTGTGCTGATTGGTGTCACCGGTTGGAAGGCGCGCGCCGCGCTGGAAAATCTCTTTAGGGAACTGCCTTCGAAAGAAAAACAGGTAg                               <  2:419912/137‑1 (MQ=255)                                                                                                                tgaaAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGCTGCCGGTTTTCAGGCGCGCGACGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGACAAACAGGTAGTCaaa                          >  1:1016574/1‑143 (MQ=255)                                                                                                                 gaaAAAATATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAaag                         >  1:1029785/1‑143 (MQ=255)                                                                                                                   aaaaaTATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAaagac                       >  1:784092/1‑142 (MQ=255)                                                                                                                   aaaaaTATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAaagaa                       <  1:796867/143‑1 (MQ=255)                                                                                                                      aaTATGCGGTATTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGATGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAGTCTCTGTAGCGAACCGCCTTCGAAAGAAAAACAGGTAGTCAaagcacgc                    >  1:635582/1‑138 (MQ=255)                                                                                                                                aTTCTCTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGCTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGCCCAcgt          >  2:324082/1‑141 (MQ=255)                                                                                                                                  tctcTAAAGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTg        <  2:920709/143‑1 (MQ=255)                                                                                                                                    tctAGAGGGCCCAGCGCGCCAGACGACGAGCGTGGGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCt      <  2:695291/143‑1 (MQ=255)                                                                                                                                        aaGTGCCCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATgg  >  2:408350/1‑143 (MQ=255)                                                                                                                                             |                                                                                                                                          CCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGG  >  NC_000913/3041494‑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: ‑