Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,448,947	A→C	T477P (ACG→CCG)	feaB →	phenylacetaldehyde dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,448,947	0	A	C	97.1%	240.4 / ‑3.0	69	T477P (ACG→CCG)	feaB	phenylacetaldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base C (37/30);  total (38/31)
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.97e-01

GAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAA  >  NC_000913/1448851‑1449041                                                                                                 |                                                                                               gAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCgg                                                                                             <  2:1674842/100‑1 (MQ=255)   atatAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGc                                                                                           <  1:1014577/100‑1 (MQ=255)     ataGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCg                                                                                         <  2:1540011/100‑1 (MQ=255)     ataGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCg                                                                                         >  1:1823998/1‑100 (MQ=255)      taGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGt                                                                                        >  2:447373/1‑100 (MQ=255)       aGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTg                                                                                       <  2:1619883/100‑1 (MQ=255)       aGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGACCGTg                                                                                       >  2:1765668/1‑100 (MQ=255)         cGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGAt                                                                                     <  2:1545472/100‑1 (MQ=255)          gATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGAtt                                                                                    >  2:670899/1‑100 (MQ=255)           aTCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGAttt                                                                                   >  1:669539/1‑100 (MQ=255)           aTCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTATGGTGGGATGAAGCAGTCAGGACCGGGCCGTGAttt                                                                                   >  2:594588/1‑100 (MQ=255)            tCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGAtttt                                                                                  >  1:1278393/1‑100 (MQ=255)            tCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGAtttt                                                                                  >  1:20152/1‑100 (MQ=255)              gCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGAATGAAGCAGTCAGGACCGGGCCGTGATTTTg                                                                                 >  2:586146/1‑99 (MQ=255)                 tACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGccc                                                                             >  2:1341927/1‑100 (MQ=255)                     ggcaggGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGAc                                                                         >  2:1105011/1‑100 (MQ=255)                              ggTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGAc                                                                <  1:217552/100‑1 (MQ=255)                                tgtgGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACgg                                                              >  1:836392/1‑100 (MQ=255)                                  tgGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGtt                                                            >  2:1501636/1‑100 (MQ=255)                                   gggTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTg                                                           >  1:1342625/1‑100 (MQ=255)                                     gTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGgt                                                         >  2:1935270/1‑100 (MQ=255)                                      tAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGgtg                                                        >  1:448513/1‑100 (MQ=255)                                      tAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGgtg                                                        <  1:975339/100‑1 (MQ=255)                                       aaaCAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGgtgt                                                       <  2:1824009/100‑1 (MQ=255)                                        aaCAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGtgtg                                                      >  1:708540/1‑100 (MQ=255)                                         aCAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGa                                                     <  1:926569/100‑1 (MQ=255)                                          cAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTAGggtgga                                                    >  1:1647531/1‑98 (MQ=255)                                            gCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAAc                                                  >  2:1057762/1‑100 (MQ=255)                                              cATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTa                                                <  2:60260/100‑1 (MQ=255)                                              cATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTa                                                >  1:1652673/1‑100 (MQ=255)                                              cATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTa                                                <  2:1206685/100‑1 (MQ=255)                                                tACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAg                                              <  1:673402/100‑1 (MQ=255)                                                 aCCTTCATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGACACTAAg                                              >  2:730460/1‑99 (MQ=255)                                                  ccTTAATTGACGCTTACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTc                                            <  1:1326535/100‑1 (MQ=255)                                                  ccTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTc                                            >  2:1193500/1‑100 (MQ=255)                                                  ccTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTc                                            <  2:1242239/100‑1 (MQ=255)                                                  ccTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTc                                            <  2:1705918/100‑1 (MQ=255)                                                      aaTTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGgtg                                        <  2:1549326/100‑1 (MQ=255)                                                       aTTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGgtgt                                       >  2:1108205/1‑100 (MQ=255)                                                        ttGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGgtgtg                                      <  2:927962/100‑1 (MQ=255)                                                          gACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGgtgtgtg                                    >  1:1756445/1‑100 (MQ=255)                                                              cTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACg                                >  2:1799396/1‑100 (MQ=255)                                                              cTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACg                                <  2:708541/100‑1 (MQ=255)                                                               tAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACgg                               <  1:1464924/100‑1 (MQ=255)                                                                aaCTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGt                              <  2:983476/100‑1 (MQ=255)                                                                 aCTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTa                             <  2:578850/100‑1 (MQ=255)                                                                    tACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTa                          <  2:199165/100‑1 (MQ=255)                                                                      ccGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAAt                        >  2:346139/1‑100 (MQ=255)                                                                       cGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATc                       >  2:1514777/1‑100 (MQ=255)                                                                       cGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATc                       >  2:1042828/1‑100 (MQ=255)                                                                          ttggttggATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTgg                    <  2:160690/100‑1 (MQ=255)                                                                          ttGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTgg                    <  2:836393/100‑1 (MQ=255)                                                                          ttGGTGGGATGAAGCAGTCAGGACCGGGCCATGATTTTGGCCCCGGCTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTgg                    <  1:951901/100‑1 (MQ=255)                                                                             gtggGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTc                 <  2:626170/100‑1 (MQ=255)                                                                                ggATGAAGCAGTGAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCt              <  1:447371/100‑1 (MQ=255)                                                                                ggATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGc               >  1:1370715/1‑99 (MQ=255)                                                                                ggATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCt              >  2:1803642/1‑100 (MQ=255)                                                                                ggATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCt              <  1:1105010/100‑1 (MQ=255)                                                                                ggATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCt              >  2:1192187/1‑100 (MQ=255)                                                                                 gATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTc             >  1:1810150/1‑100 (MQ=255)                                                                                 gATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTc             >  1:148992/1‑100 (MQ=255)                                                                                   tGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCAt           <  2:326065/100‑1 (MQ=255)                                                                                     aaGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATaa         >  1:1395754/1‑100 (MQ=255)                                                                                     aaGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATaa         >  2:493097/1‑100 (MQ=255)                                                                                      aGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAg        <  1:1341919/100‑1 (MQ=255)                                                                                       gCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTCGTAATCTGGTTCGATCATAAGt       >  2:675032/1‑100 (MQ=255)                                                                                       gCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGt       >  2:1646688/1‑100 (MQ=255)                                                                                         aGTCAGGACCGGGGCGTGAGTTTGGCCCCGACTTGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGt       <  1:1620191/98‑1 (MQ=255)                                                                                            cAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTaaaaa  >  1:258310/1‑100 (MQ=255)                                                                                                 |                                                                                               GAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAA  >  NC_000913/1448851‑1449041

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