Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,041,632	T→A	V107E (GTG→GAG)	ygfZ →	folate‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,041,632	0	T	A	98.1%	140.0 / NA	53	V107E (GTG→GAG)	ygfZ	folate‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base A (24/28);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTG‑CGTGAACCGCAGCTGACTGA‑ACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGG  >  NC_000913/3041496‑3041771                                                                                                                                           |                                                                                                                                            cATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTG‑CGTGAACCGCAGCTGACTGA‑ACTGAAAAAATATGCGGTATTCTCTAAAGAGACCAt                                                                                                                                       <  1:973819/143‑1 (MQ=255)  aTTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTG‑CGTGAACAGCAGCTGCCGGG‑ACGGAAAAAATATGCGGTATTCTCTAAAGAGACCATc                                                                                                                                      <  1:737034/143‑1 (MQ=255)    ttcgcggCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGCCGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTG‑CGTGAACCTCAGCTGACTGA‑ACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATcgc                                                                                                                                    <  2:1029971/137‑1 (MQ=255)           ccAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCCTGGATTGAACGGCGCAGCGTG‑CGTGAACCGCAGCTGACTGA‑ACTGAAACACTATGCGGTCTTCTCTAACGAGCCCATCGCGcccccc                                                                                                                             >  1:322936/1‑139 (MQ=255)            cAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTG‑CGTGAACAAAAGCACA‑TGATGCTTAAAAAATATGCGGTGTTCTCTAAAGAGACCATCGCGCCAgacg                                                                                                                            >  1:781345/1‑143 (MQ=255)             aaaGGTAAAATGTGGAGCAATTTAAGTCTGTACCGCGACGGCGATGGCTTTGTGGGGCAAAAACACAACAGAGTA‑CGGTATTCGCAGCTAACTGA‑ACTGAAAAAATATTCGGTATTCCCTAAAGAAAACATCGCGCAAgacga                                                                                                                           >  2:1046881/1‑143 (MQ=255)              aaGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTG‑CGTGAACCGCAGCTGACTGA‑ACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAgacgac                                                                                                                          >  1:723711/1‑143 (MQ=255)                ggTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTG‑CGTGAACCGCAGCTGACTGA‑ACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGAcgacga                                                                                                                        <  2:595760/143‑1 (MQ=255)                 gTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTG‑CGTGAACCGCAGCTGACTGA‑ACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAg                                                                                                                       >  1:882118/1‑143 (MQ=255)                  tAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTG‑CGTGAACCGCAGCTGACTGA‑ACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGc                                                                                                                      >  1:98409/1‑143 (MQ=255)                    aaaTGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTG‑CGTGAACCGCAGCTGACTGA‑ACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCgt                                                                                                                    <  1:517554/143‑1 (MQ=255)                         tgGAGCAATTTATGTCTGTTCCGCGCCGGCGATGGCTTTGCGTGGCTTGAACGGCGCCGAGTG‑CGTGATATGCGGCTGACTGG‑CGTGCGAAAATATGCGGTATTCTCTATAGAGACCATCGCGCCAGACGACGAGCGTGtgct                                                                                                               >  1:182586/1‑143 (MQ=255)                           gAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTG‑CGTGAACCGCAGCTGACTGA‑ACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGtgctgc                                                                                                             >  1:503737/1‑143 (MQ=255)                             gCAATTGACGTCTGTTCCGCGACGGCGATGGCTCTGCGTGGATTGAACGGCGCACCGCGTCCCGATGC‑CAGCTGCCTGA‑AGTGAAGAGATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCtt                                                                                                           <  1:1067609/143‑1 (MQ=255)                                aTTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTG‑CGTGAACCGCAGCTGACTGA‑ACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGgt                                                                                                        <  1:438890/143‑1 (MQ=255)                                aTTTACGTCTCTTCCGGGACGGCGGTGGCTTTGCATGGATTGAACGGCGCAGCGTG‑CGTGAACCGCAGCTGACTGA‑ACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGgt                                                                                                        <  2:357078/143‑1 (MQ=255)                                 tttACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTG‑CGTGAACCGCAGCTGACTGA‑ACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACCAGCGTGTGCTGCTTGgtg                                                                                                       >  1:673083/1‑143 (MQ=255)                                     cGTCTGTTCCGCGACGGCTATGGCTTTGCATGGATTGAACGGCGCAGCGTG‑CGTGAACCGCAGCTGACTGA‑ACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGc                                                                                                   <  2:398881/143‑1 (MQ=255)                                          tctCCGTGACGGCGCGTGCTTTGAATCGCTTGAATGGCGCGGCTTT‑CGTAGAGCGCAGCTGACTGA‑ATTTAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGtt                                                                                              <  2:843061/141‑1 (MQ=255)                                           gtCCGCGACGGCGATGGCTGTGCATGGAGTGAACGGCGCAGCGTG‑CGTGAACCGCAGCTGACTGA‑ACGGAAAAAATATGCGGTATGCACTAAAGAGACCATCGCGCCAGGGGACGAGCGTGAGATAATTGGTGTTGCgggact                                                                                             >  2:859774/2‑140 (MQ=255)                                            tCCGCGACGGCGAGGGCTTTGCATGGATTGAACGGCGCAGCGTG‑CGTGAACCGCAGCTGACTGA‑ACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGtttt                                                                                            <  1:571836/143‑1 (MQ=255)                                             ccGCGACGGCGATGGCTTTGCATGGACTAGCGGGCGCAGCGGG‑CGTGAACCGCAGCTGACTGC‑ACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTc                                                                                           <  1:543166/143‑1 (MQ=255)                                                 gACGGCGTTGGCTTTGCATGGATTGAGCGGCGCAGCGTG‑CGTGTACCGCAGGTGACTGA‑ACTGTCAAAATATGCGGTATTCGCTAAAGTGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGgc                                                                                       <  2:581886/143‑1 (MQ=255)                                                   cGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTG‑CGTGAACCGCAGCTGACTGA‑ACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGgcgc                                                                                     <  2:104043/143‑1 (MQ=255)                                                            ctctccATCGATTCAACGGCGCAGCGTG‑CGTCAACCGCAGCTGACTGA‑ACTGAAAAAATATGCGGTATCCTCTAAAGAGACCATCACGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGcgcgccgcgc                                                                            <  2:343629/138‑1 (MQ=255)                                                             tgtGCATGGGTGGCCCGGCGCAGCGTG‑CGTGAACCGCAGCTGACTGA‑ACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGGGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCt                                                                           <  2:923831/141‑1 (MQ=255)                                                             gctGCAGGGATTGAACGGCGCAGCGTG‑CGTGAACCGCAGCTGACTGA‑ACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCt                                                                           <  2:550344/141‑1 (MQ=255)                                                              ttGCATGGATTGCACGGCGCAGCGTG‑CGTGAACCGCAGCTGACTGA‑ACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTg                                                                          <  1:285363/143‑1 (MQ=255)                                                                       ttGAACGGCGCAGCGTG‑CGTGCACCGCAGCTGACTGA‑AATGAAAACATATGCGGTATTCTCGGCAGAGCCCATCGCGCCAGACGACGAGCGTGTGCTGATGGGTGTTGCCAGCATTCAGACGCGAGCCGCGCGGGCAACtcac                                                                 >  2:553459/1‑141 (MQ=255)                                                                          aaCGGCGCAGCGTG‑CGTGAACCGCAGCTGACTGA‑ACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCttt                                                              >  2:357641/1‑143 (MQ=255)                                                                           gcggcgcggcgTG‑CGTGAACCTCAGCTGACTGA‑ACTGAAAAAATATGCGGTATTCGCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTa                                                             <  2:928615/142‑1 (MQ=255)                                                                           aCGGCGCAGCGTG‑CGTGAACCGCAGCTGACTGA‑ACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCCTTa                                                             >  1:417859/1‑143 (MQ=255)                                                                           aCGGCGCAGCGTG‑CGCGAACCGCAGCTGACTGA‑ACTGAAACAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGCGCGTGTGCTGCTTGGTGTTGCCGGCATTCAGGCGCGCGCCGCGCGCGCACACCTCttcc                                                             >  2:46319/1‑141 (MQ=255)                                                                               cgcAGCGTG‑CATGAACCGCAGCTGACTGA‑ACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACTACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGa                                                         <  2:428004/143‑1 (MQ=255)                                                                                             aaCCGCAGCCGACTGA‑ACTGAAACAATCAGCGGCGGTCTCTACAGAGCCCAGCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCgaaa                                            <  2:805485/143‑1 (MQ=255)                                                                                                cGCAGCTGTCTGA‑ACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCgaaagaa                                         <  1:19779/143‑1 (MQ=255)                                                                                                 gCAGCTGACTGA‑ACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTAGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTATTTAGCGAACTGCCTTCgaaagaaa                                        <  2:537262/143‑1 (MQ=255)                                                                                                 gCAGCTGACTGA‑ACTGAAAAAATATGCGGTATTCGCTAAAGAGACCATCGCGCCAGGCGACGAGCGTGTGCAGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCGGGCAAATCTCTGTAGCGAACTGCCTTGggaagaaa                                        >  2:848831/1‑143 (MQ=255)                                                                                                         ctga‑actggaAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTa                                <  1:788444/143‑1 (MQ=255)                                                                                                          tga‑actgaaAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGACAGAAACACCGGTAg                               >  2:703945/1‑143 (MQ=255)                                                                                                          tga‑actgaaAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCATGTGCTGCTTGGTGTTGCCGATTTTCAGGCGCCCGCCGCGCTGGCACATCTCTTTAGCGCACTCCCTGCGCGAGAAACACATGTAg                               >  1:854124/1‑143 (MQ=255)                                                                                                           ga‑acagaaAAAATATGCGGTATGCTATAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGATTGGTGTTGGCGGTTTTCAGGCGCACACCTCGCTGGGAAATCTCTTTTGCGAACTGCCTTCGAAAGAAAAACAGGTAGt                              <  2:182586/143‑1 (MQ=255)                                                                                                                        tatGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTCGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGgcgc                  >  2:514421/1‑143 (MQ=255)                                                                                                                         gggCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGcgcg                 <  2:184725/141‑1 (MQ=255)                                                                                                                             ggTATTCTCTAAAGAGACCATCGCGCAAGACGACGAGCGTGTGCTGCTTGGTGTTGACGGTTTTCAGGCGCGCGCAGCGCTGGCAAATCTCTTTAGCGAAGGGCCCTCGACACAAACACAGGTCGGCCAAGAAGGCGCGCCCa             >  2:599893/1‑143 (MQ=255)                                                                                                                               tATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCGTTGGCGAACTGCCTTCGCAAGCAGAACAGGTAGTCAACGCAGGCGCGCCCGCt           >  2:857910/1‑143 (MQ=255)                                                                                                                                     ctAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTGCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAAGTGCCTTCGGAAGAAAAACAGGTAGTAAAAGAATGCGCGACCACTGTGCTa     >  1:349376/1‑143 (MQ=255)                                                                                                                                     cgaaaGAGACCCTCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTa     <  2:600608/141‑1 (MQ=255)                                                                                                                                       aaaGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATg   <  1:973488/143‑1 (MQ=255)                                                                                                                                       aaaGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACCTTGCTATg   >  2:55000/1‑143 (MQ=255)                                                                                                                                        aaGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATgg  >  1:235245/1‑143 (MQ=255)                                                                                                                                        aaGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCGTTAGCGAACTGCCCTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCAGTTTGCGATgg  >  2:946638/1‑143 (MQ=255)                                                                                                                                        aaGAGACCAGCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTCCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATgg  <  1:1027596/143‑1 (MQ=255)                                                                                                                                           |                                                                                                                                            CATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTG‑CGTGAACCGCAGCTGACTGA‑ACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAAAGAAAAACAGGTAGTCAAAGAAGGCGCGACCACTTTGCTATGG  >  NC_000913/3041496‑3041771

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