Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AP009048	197,395	T→C	5.2%	S284P (TCT→CCT)	yaeL →	zinc metallopeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AP009048	197,395	0	T	C	5.2%	211.8 / 2.6	77	S284P (TCT→CCT)	yaeL	zinc metallopeptidase
Reads supporting (aligned to +/- strand):  ref base T (36/37);  new base C (2/2);  total (38/39)
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.42e-01

TGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCC  >  AP009048/197325‑197462                                                                       |                                                                    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:897626/1‑70 (MQ=255)                                 ccTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCggg                                      >  1:3307857/1‑70 (MQ=255)                                 ccTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCggg                                      >  1:1993208/1‑70 (MQ=255)                                 ccTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCggg                                      >  1:237401/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:6244903/44‑1 (MQ=255)                                                     aGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTaa                                            <  1:1576187/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:5553728/71‑1 (MQ=255)                                                        cAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGt              <  1:4922697/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:6114528/53‑1 (MQ=38)                                                           gggACTCCCCTGCCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTaa                             <  1:1824813/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:6411573/1‑71 (MQ=255)                                                                    ttGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGcc  >  1:638727/1‑71 (MQ=255)                                                                    ttGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGcc  >  1:629298/1‑71 (MQ=255)                                                                     tGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTa         >  1:4125635/1‑63 (MQ=255)                                                                     tGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGcc  >  1:6731482/1‑70 (MQ=255)                                                                     tGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGcc  >  1:5821652/1‑70 (MQ=255)                                                                      gTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGtttt                <  1:5350199/55‑1 (MQ=255)                                                                       |                                                                    TGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATTGAGCC  >  AP009048/197325‑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: ‑