Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AP009048	197,392	T→C	5.0%	L283L (TTG→CTG)	yaeL →	zinc metallopeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AP009048	197,392	0	T	C	5.0%	224.6 / 2.9	79	L283L (TTG→CTG)	yaeL	zinc metallopeptidase
Reads supporting (aligned to +/- strand):  ref base T (36/39);  new base C (2/2);  total (38/41)
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.36e-01

CTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCC  >  AP009048/197322‑197462                                                                       |                                                                       cTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCt                                                                        >  1:5470426/1‑71 (MQ=255)   ttGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTg                                                                      >  1:7466262/1‑71 (MQ=255)   ttGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTg                                                                      >  1:7457029/1‑71 (MQ=255)   ttGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTg                                                                      >  1:6567932/1‑71 (MQ=255)   ttGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTg                                                                      >  1:4188778/1‑71 (MQ=255)    tgtgATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTg                                                                      <  1:945900/70‑1 (MQ=255)    tgtgATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGt                                                                     >  1:3317047/1‑71 (MQ=255)      tgATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCt                                                                   <  1:1385580/71‑1 (MQ=255)      tgATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCt                                                                   <  1:3834946/71‑1 (MQ=255)        aTGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCttt                                                                 <  1:4209517/71‑1 (MQ=255)         tGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGt                                                                     >  1:1515477/1‑66 (MQ=255)         tGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCttt                                                                 >  1:3452303/1‑70 (MQ=255)          gCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGa                                                               >  1:5793100/1‑71 (MQ=255)            ttGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACa                                                             >  1:7486774/1‑71 (MQ=255)              gTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACa                                                             <  1:4798441/69‑1 (MQ=255)               tCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCtt                                                                  <  1:2054973/63‑1 (MQ=255)                 cGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAAt                                                        >  1:7789048/1‑71 (MQ=255)                 cGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAAt                                                        <  1:3041545/71‑1 (MQ=255)                    gATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATccc                                                     >  1:3979165/1‑71 (MQ=255)                     aTAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCg                                                    <  1:3207197/71‑1 (MQ=255)                     aTAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGGCCCTTGTCTTTGACATTAATCCCg                                                    <  1:3427132/71‑1 (MQ=255)                         cccGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGga                                                  >  1:2190408/1‑69 (MQ=255)                            gggTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTT‑TCTTTGACATTAATCCCGGAGAGt                                              >  1:1452639/1‑69 (MQ=255)                             ggTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGga                                                  <  1:4899522/65‑1 (MQ=255)                             ggTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTaa                                            <  1:5833893/71‑1 (MQ=255)                             ggTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTaa                                            <  1:3343685/71‑1 (MQ=255)                             ggTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTa                                             >  1:2436857/1‑70 (MQ=255)                              gTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTaaa                                           <  1:1821578/71‑1 (MQ=255)                               tAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAAc                                          >  1:3602359/1‑71 (MQ=255)                                aaaTCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGga                                                  <  1:4636894/62‑1 (MQ=255)                                 aaTCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGa                                                               >  1:1416368/1‑48 (MQ=255)                                 aaTCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGAc                                                              <  1:523108/49‑1 (MQ=255)                                 aaTCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGAc                                                              <  1:532508/49‑1 (MQ=255)                                    ccTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTg                                                                >  1:1695216/1‑44 (MQ=255)                                    ccTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTaa                                                         >  1:204824/1‑51 (MQ=255)                                    ccTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCggg                                      >  1:237401/1‑70 (MQ=255)                                    ccTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCggg                                      >  1:1993208/1‑70 (MQ=255)                                    ccTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCggg                                      >  1:3307857/1‑70 (MQ=255)                                    ccTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCggg                                      >  1:897626/1‑70 (MQ=255)                                    ccTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCggg                                      >  1:3053792/1‑70 (MQ=255)                                     cTTAGCGTTAGAAATCGAAAGGCAGTGGAGTCCCTTGTCTTTGACa                                                             >  1:7258646/1‑46 (MQ=255)                                     cTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACa                                                             >  1:4291081/1‑46 (MQ=255)                                     cTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTa                                    <  1:1676505/71‑1 (MQ=255)                                       tAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGt                                                                     <  1:3673409/36‑1 (MQ=255)                                       tAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGt                                                                     <  1:5600713/36‑1 (MQ=255)                                       tAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTa                                             >  1:5032656/1‑60 (MQ=255)                                        aGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGga                                                  >  1:5177872/1‑54 (MQ=255)                                         gCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGa                                                               >  1:678515/1‑40 (MQ=255)                                           gTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTa                              <  1:7501580/71‑1 (MQ=255)                                            ttAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTa                              >  1:6033077/1‑70 (MQ=255)                                                aaaTCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGc                          >  1:1378508/1‑70 (MQ=255)                                                aaaTCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCg                         <  1:3119412/71‑1 (MQ=255)                                                  aTCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATgg                                <  1:3044962/62‑1 (MQ=255)                                                     gAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCTATTgg                    <  1:1684330/71‑1 (MQ=255)                                                      aaaGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGt                   <  1:5635083/71‑1 (MQ=255)                                                      aaaGACAAGGGACTCCCCTGCCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTg                     >  1:4534769/1‑69 (MQ=25)                                                      aaaGACAAGGGACTCCCCTGCCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTg                     >  1:7653797/1‑69 (MQ=25)                                                       aaGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTTAAGCGATTGGtt                  >  1:6626766/1‑71 (MQ=255)                                                       aaGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGtt                  <  1:3489967/71‑1 (MQ=255)                                                        aGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTaa                                            <  1:1576187/44‑1 (MQ=255)                                                        aGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTaa                                            <  1:6244903/44‑1 (MQ=255)                                                        aGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGttt                 <  1:2736806/71‑1 (MQ=255)                                                        aGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGttt                 >  1:4424518/1‑71 (MQ=255)                                                        aGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGt                   <  1:3372853/69‑1 (MQ=255)                                                          gCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTg               <  1:3276904/71‑1 (MQ=255)                                                           cAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGt                               >  1:7321514/1‑54 (MQ=255)                                                           cAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGt              <  1:4922697/71‑1 (MQ=255)                                                           cAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGt              <  1:5553728/71‑1 (MQ=255)                                                            aGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGt              <  1:6787180/70‑1 (MQ=255)                                                              gggAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCgg           <  1:285673/71‑1 (MQ=255)                                                              gggAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCgg           >  1:4511659/1‑71 (MQ=255)                                                              gggACTCCCCTGCCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTaa                             <  1:1824813/53‑1 (MQ=38)                                                              gggACTCCCCTGCCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTaa                             <  1:6114528/53‑1 (MQ=38)                                                                  gTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTAtt       <  1:890886/71‑1 (MQ=255)                                                                     ccTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGt                                     <  1:6425517/38‑1 (MQ=255)                                                                     ccTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAg                           >  1:3857667/1‑48 (MQ=255)                                                                     ccTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCgg           <  1:2590117/64‑1 (MQ=255)                                                                     ccTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGa     <  1:2473833/70‑1 (MQ=255)                                                                     ccTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGa     >  1:2879534/1‑70 (MQ=255)                                                                      cTTGTCTTTGACATTAATCCCGGAGAGTAAACCGTGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGa     <  1:7346222/69‑1 (MQ=255)                                                                       ttGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGcc  >  1:638727/1‑71 (MQ=255)                                                                       ttGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGcc  >  1:629298/1‑71 (MQ=255)                                                                       ttGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGcc  >  1:6411573/1‑71 (MQ=255)                                                                       |                                                                       CTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCC  >  AP009048/197322‑197462

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