Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,262,921	C→T	M305I (ATG→ATA)	actP ←	acetate permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	4,262,921	0	C	T	95.7%	165.0 / ‑0.2	46	M305I (ATG→ATA)	actP	acetate permease
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (26/18);  total (26/20)
Fisher's exact test for biased strand distribution p-value = 1.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.18e-01

TATATTCCGGATTCGCACCAACCAGCATGATCGCGCCGAAGCCGATAATAAAGGTCAGAATATAGAAGTAGCCCATAAACCCGGTGGCGTAGAACACGCTCTTACGTGCTTCGCGGGCATCGCTGACTGTAAAGAAGCGCATCA  >  REL606/4262848‑4262991                                                                          |                                                                       tataTTCCGGATTCGCACCAACCAGCATGATCGCGCCGAAGCCGATAATAAAGGTCAGAATATAGAAGTAGCCTAt                                                                      >  1:1646637/1‑76 (MQ=255)   tatTCCGGATTCGCACCAACCAGCATGATCGCGCCGAAGCCGATAATAAAGGTCAGAATATAGAAGTAGCCTATaa                                                                    <  1:1122246/76‑1 (MQ=255)   tatTCCGGATTCGCACCAACCAGCATGATCGCGCCGAAGCCGATAATAAAGGTCAGAATATAGAAGTAGCCCATaa                                                                    <  1:709319/76‑1 (MQ=255)   tatCCCGGATTCGCACCAACCAGCATGACCGCGCCGAAGCCGATAATAAAGGTCAGAATATAGAAGTAGCCTATaa                                                                    <  2:1355736/76‑1 (MQ=255)    atTCCGGATCCGCCCCAACCAGCATGATCGCGCCGAAGCCGATAATAAAGGTCAGAATATAGAAGTAGCCTATaaa                                                                   <  2:204087/76‑1 (MQ=255)       ccGGATTCGCACCAACCAGCATGATCGCGCCGAAGCCGATAATAAAGGTCAGAATATAGAAGTAGCCTATAAAccc                                                                >  1:1118963/1‑76 (MQ=255)            ttCGCACCAACCAGCATGATCGCGCCGAAGCCGATAATAAAGGTCAGAATATAGAAGTAGCCTATAAACCCGGTgg                                                           >  1:1392647/1‑76 (MQ=255)                ccccaaccaGCATGATCGCGCCGAAGCCGATAATAAAGGTCAGAATATAGAAGTAGCCTATAAACCCGGTGGCGTa                                                       <  1:1829112/74‑1 (MQ=255)                    aaccaGCATGATCGCGCCGAAGCCGATAATAAAGGTCAGAATATAGAAGTAGCCTATAAACCCGGTGGCGTAGAac                                                   >  1:367031/1‑76 (MQ=255)                    aaccaGCATGACCGCGCCGAAGCCGATAATAAAGGTCAGAATATAGAAGTAGCCTATAAACCCGGTGGCGTAGAac                                                   <  1:996803/76‑1 (MQ=255)                      ccaGCATGATCGCGCCGAAGCCGATAATAAAGGTCAGAATATAGAAGTAGCCTATAAACCCGGTGGCGTAGAacac                                                 >  1:1030472/1‑76 (MQ=255)                         gCATGATCGCGCCGAAGCCGATAATAAAGGTCAGAATATAGAAGTAGCCTATAAACCCGGTGGCGTAGAACACGct                                              >  1:1368687/1‑76 (MQ=255)                          cATGATCGCGCCGAAGCCGATAATAAAGGTCAGAATATAGAAGTAGCCTATAAACCCGGTGGCGTAGAACACGctc                                             >  1:1648429/1‑76 (MQ=255)                               tCGCGCCGAAGCCGATAATAAAGGTCAGAATATAGAAGTAGCCTATAAACCCGGTGGCGTAGAACACGCTCTTACg                                        <  2:132959/76‑1 (MQ=255)                                 gcgcCGAAGCCGATAATAAAGGTCAGAATATAGAAGTAGCCTATAAACCCGGTGGCGTAGAACACGctct                                            >  2:565715/1‑70 (MQ=255)                                 ccgcCAAAGCCGATAATAAAGGTCAGAATATAGAAGTAGCCTATAAACCCGGTGGCGTAGAACACGctct                                            <  1:565715/69‑1 (MQ=255)                                   gcCGAAGCCGATAATAAAGGTCAGAATATAGAAGTAGCCTATAAACCCGGTGGCGTAGAACACGCTCTTACGTGCt                                    >  1:1836735/1‑76 (MQ=255)                                    ccaaaGCCGATAATAAAGGTCAGAATATAGAAGTAGCCTATAAACCCGGTGGCGTAGAACACGCTCTTACGTGCtt                                   <  2:1618476/73‑1 (MQ=255)                                     cGAAGCCGATAATAAAGGTCAGCATATAGAAGTAGCCTATAAACCCGGTGGCGTAGAACACGCTCTTACGTGCTTc                                  >  2:868665/1‑76 (MQ=255)                                     cGAAGCCGATAATAAAGGTCAGAATATAGAAGTAGCCTATAAACCCGGTGGCGTAGAACACGCTCTTACGTGCTTc                                  >  2:98789/1‑76 (MQ=255)                                     cGAAGCCGATAATAAAGGTCAGAATATAGAAGTAGCCTATAAACCCGGTGGCGTAGAACACGCTCTTACGTGCTTc                                  >  2:1802261/1‑76 (MQ=255)                                      gAAGCCGATAATAAAGGTCAGAATATAGAAGTAGCCTATAAACCCGGTGGCGTAGAACACGCTCTTACGTGCTTcg                                 >  2:1558197/1‑76 (MQ=255)                                          ccGATAATAAAGGTCAGAATATAGAAGTAGCCTATAAACCCGGTGGCGTAGAACACGCTCTTACGTGCTTCGCggg                             >  1:797010/1‑76 (MQ=255)                                          ccGATAATAAAGGTCAGAATATAGAAGTAGCCTATAAACCCGGTGGCGTAGAACACGCTCTTACGTGCTTCGCggg                             >  2:102209/1‑76 (MQ=255)                                          ccGATAATAAAGGTCAGAATATAGAAGTAGCCTATAAACCCGGTGGCGTAGAACACGCTCTTACGTGCTTCGCggg                             >  2:645412/1‑76 (MQ=255)                                           cGATAATAAAGGTCAGAATATAGAAGTAGCCTATAAACCCGGTGGCGTAGAACACGCTCTTACGTGCTTCGCGGGc                            >  1:440171/1‑76 (MQ=255)                                            gATAATAAAGGTCAGAATATAGAAGTAGCCTATAAACCCGGTGGCGTAGAACACGCTCTTACGTGCTTCGCGGGCa                           <  1:1153004/76‑1 (MQ=255)                                               aagaaaGGTCAGAATATAGAAGTAGCCTATAAACCCGGTGGCGTAGAACACGCTCTTACGTGCTTCGCGGGCATCg                        <  2:894397/73‑1 (MQ=255)                                                  aaaGGTCAGAATATAGAAGTAGCCTATAAACCCGGTGGCGTAGAACACGCTCTTACGTGCTTCGCGTGCATCGCTg                     <  2:1223549/76‑1 (MQ=255)                                                  aaaGGTCAGAATATAGAAGTAGCCTATAAACCCGGTGGCGTAGAACACGCTCTTACGTGCTTCGCGGGCATCGCTg                     >  2:158048/1‑76 (MQ=255)                                                  aaaGGTCAGAATATAGAAGTAGCCTATAAACCCGGTGGCGTAGAACACGCTCTTACGTGCTTCGCGGGCATCGCTg                     >  1:633931/1‑76 (MQ=255)                                                   aaGGTCAGAATATAGAAGTAGCCTATAAACCCGGTGGCGTAGAACACGCTCTTACGTGCTTCGCGGGCATCGCTGa                    >  2:212216/1‑76 (MQ=255)                                                   aaGGTCAGAATATAGAAGTAGCCCATAAACCCGGTGGCGTAGAACACGCTCTTACGTGCTTCGCGGGCATCGCTGa                    <  2:1777445/76‑1 (MQ=255)                                                   aaGGGCAGAATATAGAAGTAGCCTATAAACCCGGTGGCGTAGAACACGCTCTTACGTGCTTCGCGGGCATCGCTGa                    <  2:782457/76‑1 (MQ=255)                                                    aGGTCAGAATATAGAAGTAGCCTATAAACCCGGTGGCGTAGAACACGCTCTTACGTGCTTCGCGGGCATCGCTGAc                   <  1:318064/76‑1 (MQ=255)                                                      gTCAGAATATAGAAGTAGCCTATAAACCCGGTGGCGTAGAACACGCTCTTACGTGCTTCGCGGGCATCGCTGACTg                 <  2:1538316/76‑1 (MQ=255)                                                        cAGAATATAGAAGTAGCCTATAAACCCGGTGGCGTAGAACACGCTCTTACGTGCTTCGCGGGCATCGCTGACTGTa               <  1:718005/76‑1 (MQ=255)                                                        cAGAATATAGAAGTAGCCTATAAACCCGGTGGCGTAGAACACGCTCTTACGTGCTTCGCGGGCATCGCTGACTGTa               >  2:1320150/1‑76 (MQ=255)                                                         aGAATATAGAAGTAGCCTATAAACCCGGTGGCGTAGAACACGCTCTTACGTGCTTCGCGGGCATCGCTGACTGTaa              <  2:1477671/76‑1 (MQ=255)                                                         aGAATATAGAAGTAGCCTATAAACCCGGTGGCGTAGAACACGCTCTTACGTGCTTCGCGGGCATCGCTGACTGTaa              <  1:1726238/76‑1 (MQ=255)                                                           aaTATAGAAGTAGCCTATAAACCCGGTGGCGTAGAACACGCTCTTACGTGCTTCGCGGGCATCGCTGACTGTCaag            >  2:1601156/1‑76 (MQ=255)                                                           aaTATAGAAGTAGCCTATAAACCCGGTGGCGTAGAACACGCTCTTACGTGCTTCGCGGGCATCGCTGACTGTAaag            >  1:1693194/1‑76 (MQ=255)                                                                ataaGTAGCCTATAAACCCGGTGGCGTAGAACACGCTCTTACGTGCTTCGCGGGCATCGCTGACTGTAAAGAAgcg       >  2:1318934/3‑76 (MQ=255)                                                                aGAAGTAGCCTATAAACCCGGTGGCGTAGAACACGCTCTTACGTGCTTCGCGGGCATCGCTGACTGTAAAGAAgcg       >  1:1687400/1‑76 (MQ=255)                                                                   aGTAGCCTATAAACCCGGTGGCGTAGAACACGCTCTTACGTGCTTCGCGGGCATCGCTGACTGTAAAGAAGCGCAt    <  1:1717372/76‑1 (MQ=255)                                                                     tAGCCTATAAACCCGGTGGCGTAGAACACGCTCTTACGTGCTTCGCGGGCATCGCTGACTGTAAAGAAGCGCATCa  >  1:1049512/1‑76 (MQ=255)                                                                          |                                                                       TATATTCCGGATTCGCACCAACCAGCATGATCGCGCCGAAGCCGATAATAAAGGTCAGAATATAGAAGTAGCCCATAAACCCGGTGGCGTAGAACACGCTCTTACGTGCTTCGCGGGCATCGCTGACTGTAAAGAAGCGCATCA  >  REL606/4262848‑4262991

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