Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	340,151	T→A	A14A (GCT→GCA)	ybaS →	predicted glutaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	340,151	0	T	A	97.1%	216.2 / ‑1.4	70	A14A (GCT→GCA)	ybaS	predicted glutaminase
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base A (30/37);  total (30/40)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

GCATCTATAAAAAGGAGTTAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTACATTCCCTTTCTGGCGAATGT  >  minE/340083‑340216                                                                     |                                                                  gCATCTATAAAAAGGAGTTAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCATa                                                                 <  1:343233/71‑1 (MQ=255)      tataAAAAGGAGTTAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCATACAcc                                                             >  1:443006/1‑70 (MQ=255)        taAAAAGGAGTTAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCATACACCCaa                                                          <  1:369798/71‑1 (MQ=255)         aaaaaGGAGTTAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCATACACCCaa                                                          <  1:1683611/70‑1 (MQ=255)          aaaaGGAGTTAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCATACACCCAAtt                                                        <  1:2341361/71‑1 (MQ=255)          aaaaGGAGTTAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCATACACCCAAtt                                                        >  1:2184395/1‑71 (MQ=255)          aaaaGGAGTTAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCATACACCCAAtt                                                        >  1:1868754/1‑71 (MQ=255)          aaaaGGAGTTAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCATACACCCAAtt                                                        >  1:157001/1‑71 (MQ=255)           aaaGGAGTTAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCATa                                                                 >  1:441423/1‑61 (MQ=255)           aaaGGAGTTAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCATACACCCAAttt                                                       <  1:40438/71‑1 (MQ=255)           aaaGGAGTTAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCATACACCCAAttt                                                       <  1:1412503/71‑1 (MQ=255)            aaGGAGTTAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCATACACCCAAttt                                                       >  1:336791/1‑70 (MQ=255)            aaGGAGTTAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCATACACCCAAttt                                                       >  1:1768325/1‑70 (MQ=255)            aaGGAGTTAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCATACACCCAAtt                                                        >  1:1550576/1‑69 (MQ=255)             aGGAGTTAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCATACACCCAATTtca                                                     <  1:1015793/71‑1 (MQ=255)             aGGAGTTAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCATACACCCAATTtca                                                     <  1:1013027/71‑1 (MQ=255)             aGGAGTTAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCATACACCCAATTtca                                                     <  1:1132104/71‑1 (MQ=255)                aGTTAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCATACACCCAATTtcactc                                                  <  1:1623408/71‑1 (MQ=255)                aGTTAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCATACACCCAATTtcact                                                   >  1:1915172/1‑70 (MQ=255)                  ttAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCATa                                                                 <  1:1540/54‑1 (MQ=255)                   tAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCATACACCCAAt                                                         <  1:2192445/61‑1 (MQ=255)                     aCAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCATa                                                                 <  1:740254/51‑1 (MQ=255)                     aCAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCATACACCCAATTTCACTCACTTa                                             <  1:1267551/71‑1 (MQ=255)                        aaaGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCATACACCCAATTTCACTCACTTAAcg                                          <  1:119074/71‑1 (MQ=255)                         aaGATGTTAGATGCAAACAAATTAGAGCAGGCAGTGGATCAGGCATACACCCAATTTCACTCACTTAAcgg                                         >  1:993672/1‑71 (MQ=255)                         aaGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCATACACCCAATTTCACTCACTTAAcgg                                         <  1:376260/71‑1 (MQ=255)                                tAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCATACACCCAATTTCACTCACTTaa                                            >  1:2205008/1‑61 (MQ=255)                                tAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCATACACCCAATTTCACTCACTTAACGGCGGACa                                   >  1:2247967/1‑70 (MQ=255)                                tAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCATACACCCAATTTCACTCACTTAACGGCGGACa                                   >  1:144812/1‑70 (MQ=255)                                 aGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACaa                                  <  1:2312332/70‑1 (MQ=255)                                 aGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCATACACCCAATTTCACTCACTTAACGGCGGACaa                                  <  1:32697/70‑1 (MQ=255)                                 aGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCATACACCCAATTTCACTCACTTAACGGCGGACaa                                  <  1:2383499/70‑1 (MQ=255)                                 aGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCATACACCCAATTTCACTCACTTAACGGCGGACaa                                  >  1:1890908/1‑70 (MQ=255)                                   aTGCAAACAAATTACAGCAGGCAGTGGATCAGGCATACACCCAATTTCACTCACTTAACGGCGGACAAAAt                               <  1:1738086/71‑1 (MQ=255)                                    tGCAAACAAATTACAGCAGGCAGTGGATCAGGCATACACCCAATTTCACTCACTTAACGGCGGACAAAATg                              >  1:1459394/1‑71 (MQ=255)                                     gCAAACAAATTACAGCAGGCAGTGGATCAGGCATACACCCAATTtcactcac                                                >  1:1948576/1‑52 (MQ=255)                                     gCAAACAAATTACAGCAGGCAGTGGATCAGGCATACACCCAATTtcactcac                                                >  1:2275047/1‑52 (MQ=255)                                     gCAAACAAATTACAGCAGGCAGTGGATCAGGCATACACCCAATTTCACTCACTTAACGGCGGACAAAAt                               >  1:1099029/1‑69 (MQ=255)                                      caaacaaaTTACAGCAGGCAGTGGATCAGGCATACACCCAATTTCACTCACTTAACGGCGGACaa                                  <  1:294972/65‑1 (MQ=255)                                      caaacaaaTTACAGCAGGCAGTGGATCAGGCATACACCCAATTTCACTCACTTAACGGCGGACAAAATGcc                            <  1:1920852/71‑1 (MQ=255)                                      caaacaaaTTACAGCAGGCAGTGGATCAGGCATACACCCAATTTCACTCACTTAACGGCGGACAAAATGcc                            <  1:1670577/71‑1 (MQ=255)                                       aaacaaaTTACAGCAGGCAGTGGATCAGGCATACACCCAATTTCACTCACTTAAcggcg                                       <  1:25212/59‑1 (MQ=255)                                       aaacaaaTTACAGCAGGCAGTGGATCAGGCATACACCCAATTTCACTCACTTAACGGCGGACAAAATGcc                            >  1:218568/1‑70 (MQ=255)                                       aaacaaaTTACAGCAGGCAGTGGATCAGGCATACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCg                           <  1:2194315/71‑1 (MQ=255)                                       aaacaaaTTACAGCAGGCAGTGGATCAGGCATACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCg                           <  1:648880/71‑1 (MQ=255)                                        aacaaaTTACAGCAGGCAGTGGATCAGGCATACACCCAATTTCACTCACTTAACGGCGGa                                     >  1:236790/1‑60 (MQ=255)                                        aacaaaTTACAGCAGGCAGTGGATCAGGCATACACCCAATTTCACTCACTTAACGGCGGa                                     >  1:602595/1‑60 (MQ=255)                                        aacaaaTTACAGCAGGCAGTGGATCAGGCATACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGa                          <  1:987522/71‑1 (MQ=255)                                         acaaaTTACAGCAGGCAGTGGATCAGGCATACACCCAATTTCACTCACTTAACGGCGGACaaaa                                >  1:1265354/1‑64 (MQ=255)                                          caaaTTACAGCAGGCAGTGGATCAGGCATACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGAt                         >  1:344891/1‑70 (MQ=255)                                          caaaATACAGCAGGCAGTGGATCAGGCATACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGAt                         >  1:928817/1‑70 (MQ=255)                                           aaaTTACAGCAGGCAGTGGATCAGGCATACACCCAATTtcactca                                                 >  1:1920484/1‑45 (MQ=255)                                           aaaTTACAGCAGGCAGTGGATCAGGCATACACCCAATTtcactca                                                 >  1:1430975/1‑45 (MQ=255)                                           aaaTTACAGCAGGCAGTGGATCAGGCATACACCCAATTTCACTCACTTAAc                                           <  1:866376/51‑1 (MQ=255)                                             aTTACAGCAGGCAGTGGATCAGGCATACACCCAATTtcactc                                                  <  1:2235131/42‑1 (MQ=255)                                             aTTACAGCAGGCAGTGGATCAGGCATACACCCAATTtcactc                                                  <  1:1613855/42‑1 (MQ=255)                                             aTTACAGCAGGCAGTGGATCAGGCATACACCCAATTtcactc                                                  <  1:1410291/42‑1 (MQ=255)                                              ttACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTACAt                    <  1:1969092/71‑1 (MQ=255)                                              ttACAGCAGGCAGTGGATCAGGCATACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTACAt                    <  1:984985/71‑1 (MQ=255)                                               tACAGCTGGGAGTGGATCAGGCATACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCg                           <  1:709587/63‑1 (MQ=255)                                                        cagTGGATCAGGCCTACACCCAATTTCACTCACTTAACGGCGGACaa                                  <  1:1433817/47‑1 (MQ=255)                                                         agTGGATCAGGCATACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTACATTCCCTTTCTgg         <  1:1291949/71‑1 (MQ=255)                                                          gTGGATCAGGCATACACCCAATTTCACTCACTTAACGGCGGACa                                   <  1:1873426/44‑1 (MQ=255)                                                           tGGATCAGGCATACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTACATTCCCTTTCTGGCg       <  1:522844/71‑1 (MQ=255)                                                             gATCAGGCATACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTACATTCCCTTTc            >  1:1346060/1‑64 (MQ=255)                                                             gATCAGGCATACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTACATTCCCTTTCTGGCGaa     >  1:1907697/1‑71 (MQ=255)                                                              aTCAGGCATACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTACATTCCCTTTCTGGCGaa     >  1:2247545/1‑70 (MQ=255)                                                              aTCAGGCATACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTACATTCCCTTTCTGGCGa      <  1:314954/69‑1 (MQ=255)                                                                cAGGCATACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGAtt                        >  1:203461/1‑49 (MQ=255)                                                                 aGGCATACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTACATTCCCTTTCTGGCGAATGt  <  1:1959436/70‑1 (MQ=255)                                                                     |                                                                  GCATCTATAAAAAGGAGTTAACAAAAGATGTTAGATGCAAACAAATTACAGCAGGCAGTGGATCAGGCTTACACCCAATTTCACTCACTTAACGGCGGACAAAATGCCGATTACATTCCCTTTCTGGCGAATGT  >  minE/340083‑340216

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