Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,152,760	G→T	N113K (AAC→AAA)	ppc ←	phosphoenolpyruvate carboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,152,760	0	G	T	98.6%	243.4 / ‑3.4	69	N113K (AAC→AAA)	ppc	phosphoenolpyruvate carboxylase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (32/36);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

ACGACGGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGGTGTTG  >  NC_000913/4152664‑4152853                                                                                                 |                                                                                              acgacgGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGtttt                                                                                            <  1:457245/100‑1 (MQ=255) acgacgGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGtttt                                                                                            >  1:1487909/1‑100 (MQ=255)  cgacgGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTtttt                                                                                           >  1:1485685/1‑100 (MQ=255)    acgGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTc                                                                                         <  1:70748/100‑1 (MQ=255)    acgGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTc                                                                                         <  2:1659050/100‑1 (MQ=255)    acgGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTc                                                                                         >  1:445161/1‑100 (MQ=255)     cgGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCa                                                                                        <  2:66734/100‑1 (MQ=255)      gggTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAg                                                                                       <  2:1130795/100‑1 (MQ=255)      gggTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAg                                                                                       <  2:708331/100‑1 (MQ=255)          aaTTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTa                                                                                   >  1:1039879/1‑100 (MQ=255)            tttCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACg                                                                                 >  2:119954/1‑100 (MQ=255)            tttCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACg                                                                                 >  2:634751/1‑100 (MQ=255)             ttCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGc                                                                                <  1:774672/100‑1 (MQ=255)                    gggTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGc                                                                         >  1:1711026/1‑100 (MQ=255)                    gggTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGc                                                                         <  2:329301/100‑1 (MQ=255)                        gAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGc                                                                     <  2:1909530/100‑1 (MQ=255)                        gAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGc                                                                     <  2:405102/100‑1 (MQ=255)                          gCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCg                                                                    >  1:178557/1‑99 (MQ=255)                          gCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGa                                                                   >  1:138652/1‑100 (MQ=255)                           ccGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGAt                                                                  <  2:195032/100‑1 (MQ=255)                            cGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATc                                                                 <  1:1345729/100‑1 (MQ=255)                             gTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCa                                                                <  2:1487917/100‑1 (MQ=255)                               gAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACt                                                              <  2:33816/100‑1 (MQ=255)                                   aCCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCg                                                          >  2:936441/1‑100 (MQ=255)                                   aCCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGAGGGCGATCACTTCCg                                                          >  2:1910410/1‑100 (MQ=255)                                    ccAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTTTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCg                                                          >  1:282936/1‑99 (MQ=255)                                    ccAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCgg                                                         >  2:1617114/1‑100 (MQ=255)                                     ctgTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCggg                                                        <  2:507747/98‑1 (MQ=255)                                       gTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGtt                                                      >  2:1505022/1‑100 (MQ=255)                                       gTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGtt                                                      <  2:1396781/100‑1 (MQ=255)                                        ttCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTg                                                     <  1:1738831/100‑1 (MQ=255)                                         tCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGc                                                    >  1:1162083/1‑100 (MQ=255)                                           cagcgacagcgaTTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTg                                                  <  1:279505/100‑1 (MQ=255)                                           cagcgacagcgaTTCCACTGCTTTTTTGATGGTATCTTCGCTCAGTTCCGGCTGTTTTTTCTGTTTTTGCGTGGTGCGGGCGATCACTTCCGGGTTGCt                                                   >  1:1172311/1‑99 (MQ=255)                                            agcgacagcgaTTCCACTGCTTTCTTGATGGTTTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCATATTTCGGGCGATCACTTCCGGCTTGCTgg                                                 >  1:1681633/1‑100 (MQ=255)                                             gcgacagcgaTTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGc                                                <  2:163598/100‑1 (MQ=255)                                              cgacagcgaTTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCa                                               <  2:585896/100‑1 (MQ=255)                                               gacagcgcTTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCa                                               <  2:575365/99‑1 (MQ=255)                                                 cagcgaTTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCt                                            >  2:262412/1‑100 (MQ=255)                                                 cagcgaTTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCt                                            >  2:461281/1‑100 (MQ=255)                                                   gcgaTTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGGGGGCGATCGCTTCCGGGTTGCTGGCAGCTTc                                          >  1:1311193/1‑100 (MQ=255)                                                    cgaTTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCg                                         <  1:79919/100‑1 (MQ=255)                                                      aTTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGcc                                       <  1:936440/100‑1 (MQ=255)                                                      aTTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGcc                                       <  1:1910397/100‑1 (MQ=255)                                                          cACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTc                                   >  2:1478937/1‑100 (MQ=255)                                                           aCTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCg                                  >  2:224652/1‑100 (MQ=255)                                                             tGCTTTTTTGATGGCGTCTGCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCTGGTTGCTGGCAGCTTCGCCTTTCGGc                                <  1:1323861/100‑1 (MQ=255)                                                              gCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGc                                >  1:1743168/1‑99 (MQ=255)                                                              gCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCg                               <  2:1534245/100‑1 (MQ=255)                                                               cTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGa                              >  2:309667/1‑100 (MQ=255)                                                               cTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGa                              <  2:623958/100‑1 (MQ=255)                                                                ttttttGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGaa                             <  2:1863084/100‑1 (MQ=255)                                                                  ttttGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAAt                           <  2:55309/100‑1 (MQ=255)                                                                  ttttGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAAt                           <  1:1505014/100‑1 (MQ=255)                                                                    ttGATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGc                         >  2:1466367/1‑100 (MQ=255)                                                                      gATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCtg                       <  1:1705977/100‑1 (MQ=255)                                                                      gATGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCtg                       <  2:1711037/100‑1 (MQ=255)                                                                       aTGGTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCtgt                      >  2:1768360/1‑100 (MQ=255)                                                                         ggTGTCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTgg                    <  2:1711504/100‑1 (MQ=255)                                                                          gtgtCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCGGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGt                   <  1:582383/100‑1 (MQ=255)                                                                           tgtCTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTa                  >  2:1496972/1‑100 (MQ=255)                                                                             tattCGCTCAGTTACGGCTGTTTTTTCAGTTCACGCAGGGTGAGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTAtt                >  2:1312862/3‑100 (MQ=255)                                                                              cTTCGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTg               >  1:1854833/1‑100 (MQ=255)                                                                                tCGCTCAGTTTCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCt             <  2:262785/100‑1 (MQ=255)                                                                                 cGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTc            <  1:1453519/100‑1 (MQ=255)                                                                                 cGCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTc            >  2:918924/1‑100 (MQ=255)                                                                                  gCTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTcg           >  2:1570218/1‑100 (MQ=255)                                                                                   cTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTcgg          >  2:494906/1‑100 (MQ=255)                                                                                   cTCAGTTCCGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTcgg          <  2:1487268/100‑1 (MQ=255)                                                                                          ccGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGGTGtt   >  2:1141184/1‑100 (MQ=255)                                                                                           cGGCTGTTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGGTGTTg  <  2:1018982/100‑1 (MQ=255)                                                                                                 |                                                                                              ACGACGGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGGTGTTG  >  NC_000913/4152664‑4152853

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