Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	216,231	A→T	5.1%	I25F (ATT→TTT)	yafT →	predicted aminopeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	216,231	0	A	T	5.1%	226.5 / 3.8	78	I25F (ATT→TTT)	yafT	predicted aminopeptidase
Reads supporting (aligned to +/- strand):  ref base A (43/31);  new base T (2/2);  total (45/33)
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.46e-01

ATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACT  >  minE/216159‑216298                                                                         |                                                                    aTGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTT‑‑‑GGATGTGCCTCTGAATCTTCTAt                                                                    >  1:1328736/1‑71 (MQ=255)   gAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGACTCTGAATCTTCTa                                                                     >  1:574633/1‑71 (MQ=255)        cAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATtgat                                                                >  1:1949272/1‑71 (MQ=255)        cAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATtgat                                                                >  1:223120/1‑71 (MQ=255)        cAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATtgat                                                                >  1:642288/1‑71 (MQ=255)         aaaaaaGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATtgat                                                                <  1:1978649/70‑1 (MQ=255)         aaaaaaGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATtgat                                                                >  1:1224211/1‑70 (MQ=255)          aaaaaGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATtgatga                                                              <  1:2846944/71‑1 (MQ=255)           aaaaGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATtgatga                                                              >  1:1578429/1‑70 (MQ=255)            aaaGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCGTCTTTtcatca                                                              >  1:871485/1‑67 (MQ=255)            aaaGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCATCTTTTCATCaa                                                             >  1:868091/1‑70 (MQ=255)             aaGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGaa                                                             >  1:2176897/1‑69 (MQ=255)               gCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAaag                                                          >  1:1114266/1‑70 (MQ=255)                   tgttgCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTTTTCATCAATagaaga                                                      <  1:2964797/70‑1 (MQ=38)                   tgttgCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTTTTCATCAATagaaga                                                      <  1:2329642/70‑1 (MQ=38)                   tgttgCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAagaaga                                                      <  1:2005521/70‑1 (MQ=255)                   tgttgCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAagaaga                                                      >  1:766097/1‑70 (MQ=255)                       gCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGaaaaa                                                  <  1:1425865/70‑1 (MQ=255)                        cATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATtgatg                                                               <  1:576789/56‑1 (MQ=255)                        cATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATtgatg                                                               <  1:1680341/56‑1 (MQ=255)                        cATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATtgatg                                                               <  1:2563706/56‑1 (MQ=255)                        cATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAg                                                <  1:433313/71‑1 (MQ=255)                        cATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAg                                                <  1:577769/71‑1 (MQ=255)                          tatGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCt                                              >  1:2376565/1‑71 (MQ=255)                           atGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTAtt                                                                   <  1:1843099/49‑1 (MQ=255)                            tgtgtgTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAg                                                >  1:2902450/1‑67 (MQ=255)                             gtgtgtTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAaagaag                                                       <  1:267521/59‑1 (MQ=255)                             gtgtgtTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTc                                             >  1:826765/1‑69 (MQ=255)                                tgtTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGt                                         <  1:2263035/70‑1 (MQ=255)                                  ttATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATtgat                                                                >  1:1375749/1‑45 (MQ=255)                                  ttATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCt                                              >  1:1423484/1‑63 (MQ=255)                                   tATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTcaca                                     >  1:1165688/1‑71 (MQ=255)                                    aTTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTa                                                                     >  1:2298184/1‑38 (MQ=255)                                    aTTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTcaca                                     >  1:949078/1‑70 (MQ=255)                                    aTTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTcaca                                     >  1:2128928/1‑70 (MQ=255)                                    aTTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTcaca                                     >  1:1570144/1‑70 (MQ=255)                                     ttCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTcacaca                                   <  1:153546/71‑1 (MQ=255)                                     ttCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTcacac                                    >  1:1250899/1‑70 (MQ=255)                                     ttCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTca                                       <  1:1821613/67‑1 (MQ=255)                                      tctcGCTGCTTGCAGGATGTGCCTCTGAATCTTCTa                                                                     >  1:183597/1‑36 (MQ=255)                                      tctcGCTGCTTGCAGGATGTGCCTCTGAATCTTCTa                                                                     >  1:1236018/1‑36 (MQ=255)                                      tctcGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAagaaga                                                      <  1:607515/51‑1 (MQ=255)                                       ctcGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATtg                                                                  >  1:2756750/1‑38 (MQ=255)                                        tcGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGaaaa                                                           >  1:720511/1‑44 (MQ=255)                                          gctgctTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGt                               <  1:1410746/70‑1 (MQ=255)                                           ctgctTGCAGGATGTGCCTCTGAATCTTCTATtgatg                                                               <  1:2294811/37‑1 (MQ=255)                                           ctgctTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGt                               <  1:1670682/69‑1 (MQ=255)                                             gctTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAat                            >  1:241886/1‑70 (MQ=255)                                              ctTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAatat                          <  1:2663473/71‑1 (MQ=255)                                              ctTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAatat                          <  1:368636/71‑1 (MQ=255)                                              ctTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAata                           <  1:2621872/70‑1 (MQ=255)                                                tGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATAtt                         >  1:1221497/1‑70 (MQ=255)                                                 gCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATtaa                       >  1:550369/1‑71 (MQ=255)                                                     gATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATtaataa                    >  1:1705292/1‑70 (MQ=255)                                                     gATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATtaataa                    >  1:1971700/1‑70 (MQ=255)                                                     gATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATaaa                   <  1:594244/71‑1 (MQ=255)                                                      aTGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCt                                              >  1:1562737/1‑43 (MQ=255)                                                       tgtgCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGt                                         >  1:519683/1‑47 (MQ=255)                                                         tgCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAaca               >  1:1042311/1‑71 (MQ=255)                                                         tgCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAaca               >  1:2781508/1‑71 (MQ=255)                                                         tgCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAaca               >  1:1215413/1‑71 (MQ=255)                                                          gCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACaaa                                 >  1:2893703/1‑52 (MQ=255)                                                          gCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAaca               >  1:174449/1‑70 (MQ=255)                                                           ccTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACt             >  1:1499291/1‑71 (MQ=255)                                                               tGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTca                                       >  1:748722/1‑41 (MQ=255)                                                               tGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTcc           <  1:52013/69‑1 (MQ=255)                                                                gAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGt                                         >  1:1652283/1‑38 (MQ=255)                                                                gAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTccc          >  1:1648841/1‑69 (MQ=255)                                                                 aaTCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCa        <  1:2253016/70‑1 (MQ=255)                                                                 aaTCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCa        <  1:1678235/70‑1 (MQ=255)                                                                 aaTCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAg       >  1:2713396/1‑71 (MQ=255)                                                                  aTCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATta                        <  1:1516782/53‑1 (MQ=255)                                                                  aTCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAg       <  1:2524399/70‑1 (MQ=255)                                                                  aTCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAg       >  1:1741174/1‑70 (MQ=255)                                                                  aTCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGc      <  1:700931/71‑1 (MQ=255)                                                                     ttctATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATta                        >  1:711365/1‑50 (MQ=255)                                                                       ctATTGATGAAAAGAAGAAAAAAGCTCAAGTcaca                                     >  1:2720681/1‑35 (MQ=255)                                                                       ctATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATtaat                      <  1:65403/50‑1 (MQ=255)                                                                       ctATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATtaat                      <  1:798363/50‑1 (MQ=255)                                                                       ctATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATtaat                      <  1:803694/50‑1 (MQ=255)                                                                       ctATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACt  >  1:363555/1‑70 (MQ=255)                                                                       ctATTGATGAAAAGAAGAAAAAAGCTAAAGTCACACAAAGTAATATTAATAAAAACACt             >  1:1471286/1‑59 (MQ=255)                                                                         |                                                                    ATGAATTCAAAAAAGCTTTGTTGCATATGTGTGTTATTCTCGCTGCTTGCAGGATGTGCCTCTGAATCTTCTATTGATGAAAAGAAGAAAAAAGCTCAAGTCACACAAAGTAATATTAATAAAAACACTCCCCAGCAACT  >  minE/216159‑216298

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