Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,093,040	G→A	A215T (GCT→ACT)	hisC →	histidinol‑phosphate aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,093,040	0	G	A	95.7%	132.2 / NA	47	A215T (GCT→ACT)	hisC	histidinol‑phosphate aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (23/22);  minor base C (0/2);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATTTTCGCACCCTGCTGGAGTTAACCCGCGGTAAGGCGATTGTGGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAGCTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCGGCCCAGGCGTTAAG  >  NC_000913/2092896‑2093170                                                                                                                                                 |                                                                                                                                   gATTTTCGCACCCTGCTGGAGTTAACCCGCGGTAAGGCGATTGTGGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCCACCTGGCGAAAACttt                                                                                                                                >  2:495100/1‑149 (MQ=255) gATTTTCGCACCCTGCTGGAGTTAACCCGCGGTAAGGCGATTGTGGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACttt                                                                                                                                <  2:114890/149‑1 (MQ=255)  ttgttCGCACCCTGCTGGAGTTAACCCGCGGTAAGGCGATTGTGGTTGCCGATGAAGCCTATATCGAGCTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACtttt                                                                                                                               <  1:236962/146‑1 (MQ=255)    tttCGCACCCTGCTGGAGTTAACCCGCGGTAAGGCGATTGTGGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCCATTTTACGCACACTGTCGAAAACTTTTGc                                                                                                                             >  1:514480/1‑149 (MQ=255)          gtccTGCTGGAGTTAAGCCGCGGTGAGGCGATTGTGGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTTGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGc                                                                                                                       <  1:23965/147‑1 (MQ=255)          aCCCTGCTGGAGTTAACCCGCGGTAAAGCGATTGTGGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCCGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCCCTGGc                                                                                                                       >  2:42042/1‑149 (MQ=255)            ccTGCTGGAGTTAACCCGCGGTAAGGCGATTGTGGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCgg                                                                                                                     >  1:302036/1‑149 (MQ=255)              tgcggGAGTTAACCCGCGGTAAGGCGATTGTGGTTGCCGATGCAGCCTAGATCGAGTTTTGCCCGCAGGCATCGCTGGCCGGCTGGCTGGCGGAATATCCGCACCTGGCGATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCgggg                                                                                                                   <  2:378395/149‑1 (MQ=255)                 tagAGTTAACCCGCGGTAAGGCGATTGTGGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCtt                                                                                                                >  2:103835/3‑149 (MQ=255)                 tagAGTTAACCCGCGGTAAGGCGATTGTGGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCtt                                                                                                                <  1:103835/147‑1 (MQ=255)                                   aGGGGATGGGGGTTCCGGATGAAGCCTATACCGCGGGTTGCCACCAGGCATCGCGGCGTGGCTGGCTGGCGGAATGTCGGCACCTGGCTATTTTACGCACGCTGTCGAAACCTTTTGCTCTGGCGGGGCTTCGTCGCGGATTTACGCTg                                                                                              <  2:630810/149‑1 (MQ=255)                                            tgGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACgaa                                                                                     <  1:140743/149‑1 (MQ=255)                                            tgGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATAACCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGACTCGTTTCGGATTTACGCTGGCAAACgaa                                                                                     >  1:138699/1‑149 (MQ=255)                                             ggTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACgaag                                                                                    <  1:220458/149‑1 (MQ=255)                                             ggTTGCCGATGAAGCCTATATCGAGCTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACgaag                                                                                    <  2:199130/149‑1 (MQ=255)                                               ttGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCGCCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACgaagaa                                                                                  >  1:456334/1‑149 (MQ=255)                                               ttGCCGAGCTAGCCTATATCGAGCTTTGCCCGCAGGCATCGCTGGCGGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACgaagaa                                                                                  <  2:544388/149‑1 (MQ=255)                                                   cGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAAGGAAGAGGtca                                                                              >  2:58431/1‑149 (MQ=255)                                                   cGATGAAGCCTATATCGAGTTTTGCCCGCAGGCAGCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGtca                                                                              <  1:333423/149‑1 (MQ=255)                                                     aTGAGGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATACGCACCTGGCTATTTTACACACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGCATTTACGCTGGCAAACCAAGAAGacacc                                                                            >  2:77966/1‑147 (MQ=255)                                                           ccTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACctg                                                                      >  2:194772/1‑149 (MQ=255)                                                              aaaccGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACctgctg                                                                   <  1:244871/145‑1 (MQ=255)                                                                atCGGGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCtgat                                                                 <  2:663591/149‑1 (MQ=255)                                                                  cGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCCGCtgatga                                                               >  1:325861/1‑149 (MQ=255)                                                                   gAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTACTTTACGCACACTGGCGCAAACTTTTGCCCTGGCGGGGCCTCGTTGCGGATTTACGCTGGCAAACCAAGAAGGCATCAACccgcttctgca                                                              >  2:568646/1‑143 (MQ=255)                                                                       tttGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTg                                                          <  2:458195/149‑1 (MQ=255)                                                                          gCCCGCCGGACTGGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTCTTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATc                                                       <  2:456334/149‑1 (MQ=255)                                                                              gCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGACTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGccc                                                   <  2:423942/149‑1 (MQ=255)                                                                                  gCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTGATCAACCTGCTGATGAAAGTGATCGCCCCCTa                                               >  2:459224/1‑149 (MQ=255)                                                                                       gctggctggctggctgGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGc                                          >  1:559110/1‑149 (MQ=255)                                                                                        ctggctggctggctggCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGCTGAAAGGGATGGCCCCCTACCCGct                                         >  1:183346/1‑149 (MQ=255)                                                                                         tggctggctggctggcGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGctc                                        >  1:224960/1‑149 (MQ=255)                                                                                         tggcgggcgggctggaCGCAGATCCGCACCTGGCTATTTCATCCACACTGTTGAAACCTTTTGCTCTGGCGGGGCTTCGTGGCGGCTTTTCGCTGGCAAACGAAGTAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTCCCCGctc                                        <  2:324921/149‑1 (MQ=255)                                                                                           gctggctggctggcGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCtctc                                      <  2:559110/149‑1 (MQ=255)                                                                                             tggctggctggtGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGa                                    <  2:188679/149‑1 (MQ=255)                                                                                               gctggcttgcGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCCTCGTTGCGGATTTACGCTGGCAAACGAACAAGTCATCACCCTGCTGATCAAAGTCATCCCCCCCTACCCCCCCTCGACg                                  >  2:238363/1‑149 (MQ=255)                                                                                                ctggctggcGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGc                                 <  1:637222/149‑1 (MQ=255)                                                                                                   gctggcGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCACTCTCGagcacgg                              >  2:202712/1‑143 (MQ=255)                                                                                                     tggcGGAATATCCGCACCTGGCTATTTTACGCCCACTGTCGAAAACTTTTGCTCTGGCGGAGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGtt                            >  2:581443/1‑149 (MQ=255)                                                                                                        cGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGcc                         <  2:514480/149‑1 (MQ=255)                                                                                                                   ccccTGGCTTTTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCGGc              <  1:625386/147‑1 (MQ=255)                                                                                                                     ccTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTACTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCGGccc            <  1:33771/149‑1 (MQ=255)                                                                                                                      cTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGCTTGCCGACATTGCGGCCCa           <  1:593712/149‑1 (MQ=255)                                                                                                                        ggCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCGGCCCAgg         >  1:43323/1‑149 (MQ=255)                                                                                                                           tATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGGAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCCCTCGACGCCGGTTGCCGACATTGCGGcccgcgggt      >  2:13410/1‑143 (MQ=255)                                                                                                                             ttttACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCGGCCCAGGCGTTa    >  2:607007/1‑149 (MQ=255)                                                                                                                               ttACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCGGCCCAGGCGTGAAg  >  1:541216/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                   GATTTTCGCACCCTGCTGGAGTTAACCCGCGGTAAGGCGATTGTGGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAGCTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCGGCCCAGGCGTTAAG  >  NC_000913/2092896‑2093170

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