Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,791,743	C→G	10.5%	N157K (AAC→AAG)	gabD →	succinate‑semialdehyde dehydrogenase I, NADP‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,791,743	0	C	G	10.5%	191.8 / 11.1	57	N157K (AAC→AAG)	gabD	succinate‑semialdehyde dehydrogenase I, NADP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (27/24);  new base G (3/3);  total (30/27)
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.20e-01

CTCCTTTATTGAGTGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTA  >  NC_000913/2791611‑2791891                                                                                                                                     |                                                                                                                                                      ctcctTTATTGAGTGGGTTGCCGAAGAAGGCAAACGTATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGAATGTGATCAAGCAGCCGATTGGCGTCACCGCGCTTATCATGCCGTGGAACTTCC‑CGGCGGCgatga                                                                                                                                      <  3:590661/149‑1 (MQ=255)  tcctTTATTGAGTGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCgatgat                                                                                                                                     >  3:168418/1‑149 (MQ=255)    ctTTATTGAGTGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTa                                                                                                                                   >  1:141360/1‑149 (MQ=255)    ctTTATTGAGTGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTa                                                                                                                                   >  3:221618/1‑149 (MQ=255)       tATTGAGTGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTAccc                                                                                                                                <  4:221618/149‑1 (MQ=255)         ttGAGTGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTACCCGc                                                                                                                              >  1:158488/1‑149 (MQ=255)              tGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTACCCGCAAAGc                                                                                                                         <  4:244212/149‑1 (MQ=255)               ggTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTACCCGCAAAGcc                                                                                                                        <  3:337444/149‑1 (MQ=255)                            agGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTACCCGCAAAGcc                                                                                                                        <  4:560166/136‑1 (MQ=255)                            agGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTACCCGCAAAGcc                                                                                                                        >  3:560166/1‑136 (MQ=255)                            agGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTgg                                                                                                           <  1:127148/149‑1 (MQ=255)                             ggCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGgc                                                                                                          <  3:562459/149‑1 (MQ=255)                                aaaCGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCTCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGgcagc                                                                                                       <  3:22485/149‑1 (MQ=255)                                aaaCGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCG‑GGAAGTTCCAcggcg                                                                                                                                             <  4:323110/111‑1 (MQ=255)                                aaaCGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCG‑GGAAGTTCCAcggcg                                                                                                                                             >  3:323110/1‑111 (MQ=255)                                   cGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAgg                                                                                                    <  4:168418/149‑1 (MQ=255)                                               gACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGt                                                                                        >  2:13177/1‑149 (MQ=255)                                               gACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCG‑GGAAGTTCCAcggcgtgata                                                                                                                                        >  4:474296/1‑96 (MQ=255)                                               gACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCG‑GGAAGTTCCAcggcgtgata                                                                                                                                        <  3:474296/101‑6 (MQ=255)                                                    cATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCga                                                                                                                                         <  1:32758/95‑1 (MQ=255)                                                    cATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCga                                                                                                                                         >  2:32758/1‑95 (MQ=255)                                                    cATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGa                                                                                   >  1:115746/1‑149 (MQ=255)                                                     aTTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGaa                                                                                  <  1:13177/149‑1 (MQ=255)                                                     aTTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGaa                                                                                  <  3:91453/149‑1 (MQ=255)                                                          tGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCg                                                                             >  3:76186/1‑149 (MQ=255)                                                             tcatcaGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCa                                                                          >  3:168451/1‑149 (MQ=255)                                                               atcaGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGc                                                                                      >  3:229704/1‑135 (MQ=255)                                                               atcaGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGc                                                                                      <  4:229704/135‑1 (MQ=255)                                                               atcaGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGt                                                                        >  1:19419/1‑149 (MQ=255)                                                               atcaGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGt                                                                        >  3:357255/1‑149 (MQ=255)                                                                tcaGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGACAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTc                                                                       <  1:361380/149‑1 (MQ=255)                                                                  aGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGCTTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAg                                                                     <  4:357255/149‑1 (MQ=255)                                                                  aGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAg                                                                     <  2:19419/149‑1 (MQ=255)                                                                   ggCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGa                                                                    >  2:7372/1‑149 (MQ=255)                                                                       gATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCga                                                                                                                                         >  3:13238/1‑76 (MQ=255)                                                                       gATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCga                                                                                                                                         <  4:13238/76‑1 (MQ=255)                                                                       gATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCG‑GGAAGTTCCAcggcgt                                                                                                                                            <  3:222272/73‑2 (MQ=37)                                                                       gATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCG‑GGAAGTTCCAcggcgt                                                                                                                                            >  4:222272/1‑72 (MQ=37)                                                                          aaaCGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGtt                                                             >  2:393597/1‑149 (MQ=255)                                                                            aCGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTACCCGCAAAGCCGGTcc                                                                                                                   <  4:601142/93‑1 (MQ=255)                                                                            aCGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTACCCGCAAAGCCGGTcc                                                                                                                   >  3:601142/1‑93 (MQ=255)                                                                              gCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTAccc                                                                                                                                <  1:618564/78‑1 (MQ=255)                                                                              gCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTAccc                                                                                                                                >  2:618564/1‑78 (MQ=255)                                                                                cTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTgc                                                       >  4:373264/1‑149 (MQ=255)                                                                                 tGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTgcg                                                      >  4:478462/1‑149 (MQ=255)                                                                                        tATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTgcgctggcgc                                               >  1:38075/1‑149 (MQ=255)                                                                                         aTCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTgcgctggcgct                                              >  1:161518/1‑149 (MQ=255)                                                                                         aTCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTgcgctggcgct                                              >  2:30792/1‑149 (MQ=255)                                                                                          tCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTgcgctggcgctg                                             >  3:133071/1‑149 (MQ=255)                                                                                            aaGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTACCCGc                                                                                                                              >  4:485776/1‑66 (MQ=255)                                                                                            aaGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTACCCGc                                                                                                                              <  3:485776/66‑1 (MQ=255)                                                                                                        ggCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGAt                               <  2:38075/149‑1 (MQ=255)                                                                                                         gCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGCTGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATc                              >  2:509398/1‑149 (MQ=255)                                                                                                                cgcgGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCgg                       <  2:140136/149‑1 (MQ=255)                                                                                                                cgcgGCTATCACGCCGTGGAACATCC‑CGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGCTCCGCGCgg                       <  1:535555/149‑1 (MQ=255)                                                                                                                 gcgGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCggg                      <  3:478462/149‑1 (MQ=255)                                                                                                                                 ggAACTTCC‑CGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTa      >  2:401240/1‑149 (MQ=255)                                                                                                                                     cTTCC‑CGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTa  >  4:109802/1‑149 (MQ=255)                                                                                                                                     |                                                                                                                                                      CTCCTTTATTGAGTGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCC‑CGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTA  >  NC_000913/2791611‑2791891

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑