Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,993,572	G→A	24.5%	V289M (GTG→ATG)	nupG →	nucleoside transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,993,572	0	G	A	24.5%	101.9 / 30.4	53	V289M (GTG→ATG)	nupG	nucleoside transporter
Reads supporting (aligned to +/- strand):  ref base G (18/22);  new base A (9/4);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 2.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

CTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGG  >  minE/1993503‑1993641                                                                      |                                                                      cTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATtgt                                                                      >  1:1053670/1‑71 (MQ=255) cTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATtgt                                                                      >  1:135429/1‑71 (MQ=255)   gTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTgg                                                                    >  1:888908/1‑71 (MQ=255)        ttCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTATGGCGTg                                                                >  1:1992019/1‑70 (MQ=255)         tCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTATGGCGTGGa                                                              >  1:1285117/1‑71 (MQ=255)         tCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTATGGCGTGGa                                                              >  1:1685703/1‑71 (MQ=255)          cTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTATGGCg                                                                  >  1:863574/1‑66 (MQ=255)          cTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTATGGCGTGGa                                                              <  1:595192/70‑1 (MQ=255)          cTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTATGGCGTGGa                                                              <  1:2160796/70‑1 (MQ=255)               cATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTg                                                         <  1:123800/70‑1 (MQ=255)               cATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTg                                                         <  1:2196343/70‑1 (MQ=255)                 tCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGt                                                      >  1:1855538/1‑71 (MQ=255)                     gTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGa                                                              <  1:2014886/59‑1 (MQ=255)                       tcttctTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTgc                                                 >  1:980779/1‑70 (MQ=255)                           cttAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCt                                              <  1:813331/69‑1 (MQ=255)                           cttAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGt                                            >  1:2016985/1‑71 (MQ=255)                           cttAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCATGCGTTTTGCGCt                                              <  1:182077/69‑1 (MQ=255)                           cttAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCt                                              <  1:219214/69‑1 (MQ=255)                             tAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTTTTGTGGCg                                                                  <  1:1589937/47‑1 (MQ=255)                             tAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGttt                                          >  1:885331/1‑71 (MQ=255)                              aaGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGt                                            <  1:920758/68‑1 (MQ=255)                               aGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTg                                         <  1:1947564/70‑1 (MQ=255)                                gCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCt                                       <  1:1587855/71‑1 (MQ=255)                                   gCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTAc                                    <  1:614656/71‑1 (MQ=255)                                   gCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTAc                                    <  1:1477584/71‑1 (MQ=255)                                      aCGGTATTAAGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGc                                 >  1:55835/1‑71 (MQ=255)                                      aCGGTATTAAGAACGTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGc                                 >  1:1483545/1‑71 (MQ=255)                                       cGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGc                                 <  1:97382/70‑1 (MQ=255)                                           aTTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTAc                                    <  1:365416/63‑1 (MQ=255)                                            ttAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGc                                        <  1:1751601/58‑1 (MQ=255)                                            ttAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGAccc                            <  1:1646911/70‑1 (MQ=255)                                            ttAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGAccc                            <  1:1399079/70‑1 (MQ=255)                                              aaGAACTTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTAc                                    >  1:285040/1‑60 (MQ=255)                                              aaGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTg                                             >  1:721024/1‑51 (MQ=255)                                              aaGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTAc                                    >  1:511178/1‑60 (MQ=255)                                              aaGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTAc                                    >  1:98346/1‑60 (MQ=255)                                              aaGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTAc                                    >  1:501167/1‑60 (MQ=255)                                               aGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGAc                         <  1:1511725/70‑1 (MQ=255)                                                    gTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGtt                   >  1:1102717/1‑71 (MQ=255)                                                    gTAATGATGATCAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGtt                   >  1:1334015/1‑71 (MQ=255)                                                     tAATGATGATCAGTATTGTGGCGTGGATCCTGCGttt                                                    >  1:1767181/1‑37 (MQ=255)                                                     tAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTc                  <  1:2152439/71‑1 (MQ=255)                                                       atgatgATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCgg                <  1:180751/71‑1 (MQ=255)                                                       atgatgATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCgg                <  1:1245958/71‑1 (MQ=255)                                                            gatCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGgtactg           >  1:685203/1‑71 (MQ=255)                                                             atCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGgtactg           <  1:1922535/70‑1 (MQ=255)                                                               cAGTATTATGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGAACCGACTCCGTTc                  >  1:1653750/1‑61 (MQ=255)                                                                aGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGgtactgtact       >  1:954082/1‑71 (MQ=255)                                                                  tATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCgg                <  1:580458/60‑1 (MQ=255)                                                                   aTTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCg                 >  1:655349/1‑58 (MQ=255)                                                                    ttGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGtt                   <  1:906456/55‑1 (MQ=255)                                                                    ttGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCgg                >  1:785270/1‑58 (MQ=255)                                                                    ttGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCgg                >  1:1968654/1‑58 (MQ=255)                                                                    ttGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTActgctg   >  1:167178/1‑71 (MQ=255)                                                                    ttGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTActgctg   >  1:1369352/1‑71 (MQ=255)                                                                      gtgGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTgg  >  1:2114310/1‑70 (MQ=255)                                                                      |                                                                      CTGTTCATTCTGACCATCCCGTTCTTCTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGG  >  minE/1993503‑1993641

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