Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	798,273	(A)7→8	coding (1003/1146 nt)	ECB_00739 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	798,266	1	.	A	91.4%	257.7 / 16.9	70	coding (1010/1146 nt)	ECB_00739	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (3/3);  new base A (29/35);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TTACCTATATTTATCATTGGCTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGG‑AAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCTGCCTTATAATAGAAGAACAAAGCAAG  >  REL606/798191‑798339                                                                             |                                                                          ttACCTATATTTATCATTGGCTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGATAATAACCAAgg                                                                            >  2:961195/1‑76 (MQ=255)  tACCTATATTTATCATTGGCTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGG‑a                                                                          >  1:985246/1‑76 (MQ=255)   aCCTATATTTATCATTGGCTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGG‑aa                                                                         <  2:851079/76‑1 (MQ=255)   aCCTATATTTATCATTGGCTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGG‑aa                                                                         <  1:293657/76‑1 (MQ=255)    ccTATATTTATCATTGGCTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGG‑aaa                                                                        <  1:1141893/76‑1 (MQ=255)      tataTTTATCATTGGCTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGG‑aaaaa                                                                      <  1:569565/76‑1 (MQ=255)       atatTTATCATTGGCTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGG‑aaaaaa                                                                     <  2:339118/76‑1 (MQ=255)       atatTTATCATTGGCTGTTCAATGAGGATATAGATTGGATATGATGCAGTTATTGAGATAATAACCAAGG‑aaaaaa                                                                     <  2:339150/76‑1 (MQ=255)        tatTTATCTTTGGCTTTTAAATGAGTATATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGG‑aaaaaaa                                                                    <  1:456783/76‑1 (MQ=255)        tatTTATCATTGGCTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGG‑aaaaaaa                                                                    <  1:456774/76‑1 (MQ=255)         atTTATCATTGGCTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGG‑aaaaaaaa                                                                   <  2:655822/76‑2 (MQ=255)          tttATCATTGGCTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGG‑aaaaaaaac                                                                  >  2:154848/1‑74 (MQ=255)          tttATCATTGGCTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGG‑aaaaaaaac                                                                  >  1:1350062/1‑74 (MQ=255)           ttATCATTGGCTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGG‑aaaaaaaaca                                                                 >  1:804666/1‑73 (MQ=255)              tCATTGGCTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAg                                                               <  1:609256/76‑1 (MQ=255)               cATTGGCTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAGa                                                              <  2:1256417/76‑1 (MQ=255)                  tGGCTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAGaaaa                                                           >  2:579454/1‑76 (MQ=255)                   ggCTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAAt                                                          >  1:1532496/1‑76 (MQ=255)                   ggCTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAAt                                                          >  1:1536524/1‑76 (MQ=255)                     cTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATcc                                                        >  2:997690/1‑76 (MQ=255)                     cTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATcc                                                        >  2:412193/1‑76 (MQ=255)                        ttCAATGAGAATATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATccacc                                                     >  1:494942/1‑76 (MQ=255)                        ttCAATGAGAATATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATccacc                                                     >  1:506014/1‑76 (MQ=255)                         tCAATGAGAATATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATccacca                                                    >  1:206105/1‑76 (MQ=255)                         tCAATGAGAATATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATccacca                                                    <  1:1607871/76‑1 (MQ=255)                          cAATGAGAATATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATCCACCAg                                                   >  2:404068/1‑76 (MQ=255)                           aaTGAGAATATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGt                                                  >  2:812311/1‑76 (MQ=255)                            aTGAGAATATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGtt                                                 >  2:11977/1‑76 (MQ=255)                            aTGAGAATATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGtt                                                 <  1:784088/76‑1 (MQ=255)                            aTGAGAATATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGtt                                                 >  2:11986/1‑76 (MQ=255)                             tGAGAATATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGTTa                                                >  1:416253/1‑76 (MQ=255)                              gagaATATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGTTAt                                               <  2:1643517/76‑1 (MQ=255)                              gagaATATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGTTAt                                               >  1:1027501/1‑76 (MQ=255)                              gagaATATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGTTAt                                               <  2:471179/76‑1 (MQ=255)                               agaATATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGTTATg                                              >  2:408943/1‑76 (MQ=255)                                 aaTATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGTTATGcc                                            <  2:394699/76‑1 (MQ=255)                                 aaTATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGTTATGcc                                            >  1:1483407/1‑76 (MQ=255)                                 aaTATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGG‑AAAAAAACACAGAAAATCCACCAGTTATGCCa                                           >  2:56962/1‑76 (MQ=255)                                 aaTATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGG‑AAAAAAACACAGAAAATCCACCAGTTATGCCa                                           >  2:1216290/1‑76 (MQ=255)                                  atatatatTGGATATGATGCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGTTATGCCa                                           <  2:805933/76‑1 (MQ=255)                                    atatatTGGATATGATGCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGTTATGCCAAt                                         <  2:1665896/76‑1 (MQ=255)                                      atatTGGATATGATGCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGTTATGCCAAttt                                       <  1:821861/76‑1 (MQ=255)                                      atatTGGATATGATGCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGTTATGCCAAttt                                       <  1:178314/76‑1 (MQ=255)                                      atatTGGATATGATGCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGTTATGCCAAttt                                       <  2:1251021/76‑1 (MQ=255)                                       tatTGGATATGATGCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGTTATGCCAAtttt                                      >  1:1028375/1‑75 (MQ=255)                                       tatTGGATATGATGCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTc                                      <  1:485219/76‑1 (MQ=255)                                       tatTGGATATGATGCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTc                                      <  1:1044062/76‑1 (MQ=255)                                        atTGGATATGATGCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCa                                     >  1:519508/1‑76 (MQ=255)                                         ttGGATATGATGCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAg                                    <  2:1470488/76‑1 (MQ=255)                                             atatGATGCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAg                                <  1:57117/76‑1 (MQ=255)                                             atatGATGCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAg                                <  2:1439686/76‑1 (MQ=255)                                                tgatgCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTc                             >  2:1280357/1‑76 (MQ=255)                                                 gatgCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCt                            <  1:1334207/76‑1 (MQ=255)                                                 gatgCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCt                            <  1:1334225/76‑1 (MQ=255)                                                 gatgCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCt                            <  2:744813/76‑1 (MQ=255)                                                 gatgCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCt                            <  2:1263240/76‑1 (MQ=255)                                                  atgCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCTg                           >  2:485884/1‑76 (MQ=255)                                                    gCAGTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCTGcc                         >  2:844640/1‑76 (MQ=255)                                                     cAGTTATTGAGATAATAACCAAGG‑AAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCTGCCtt                       <  1:9571/76‑1 (MQ=255)                                                      aGTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCTGCCtt                       >  2:650844/1‑76 (MQ=255)                                                      aGTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCTGCCtt                       <  1:573928/76‑1 (MQ=255)                                                       gTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCTGCCTta                      <  1:1655819/76‑1 (MQ=255)                                                       gTTATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCTGCCTta                      >  2:45927/1‑76 (MQ=255)                                                        ttATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCTGCCTtat                     <  2:1725443/76‑1 (MQ=255)                                                        ttATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCTGCCTtat                     <  2:1352195/76‑1 (MQ=255)                                                         tATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCTGCCTtata                    >  2:1315843/1‑76 (MQ=255)                                                         tATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCTGCCTtata                    >  2:105407/1‑76 (MQ=255)                                                         tATTGAGATAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCTGCCTtata                    >  2:1315866/1‑76 (MQ=255)                                                          aTTGAGATAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCTGCCTTataa                   >  2:97079/1‑76 (MQ=255)                                                             gagaTAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCTGCCTTATAATag                <  2:1449627/76‑1 (MQ=255)                                                              agaTAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCTGCCTTATAATaga               <  1:1325392/76‑1 (MQ=255)                                                              agaTAATAACCAAGGAAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCTGCCTTATAATaga               >  2:353214/1‑76 (MQ=255)                                                              agaTAATAACCAAGG‑AAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCTGCCTTATAATagaa              <  1:1545426/76‑1 (MQ=255)                                                                ataataACCAAGGAAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCTGCCTTATAATagaag             <  1:1420766/76‑1 (MQ=255)                                                                ataataACCAAGGAAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCTGCCTTATAATagaag             <  2:94229/76‑1 (MQ=255)                                                                 taataaCCAAGGAAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCTGCCTTATAATagaaga            <  1:124900/76‑1 (MQ=255)                                                                   ataaCCAAGGAAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCTGCCTTATAATAGAAGAAc          <  1:577389/76‑1 (MQ=255)                                                                    taaCCAAGG‑AAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCTGCCTTATAATAGAAGAACaa        <  2:769264/76‑1 (MQ=255)                                                                     aaCCAAGGAAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCTGCCTTATAATAGAAGAACaa        <  1:249390/76‑1 (MQ=255)                                                                       ccAAGGAAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCTGCCTTATAATAGAAGAACAAAg      <  1:503632/76‑1 (MQ=255)                                                                        cAAGGAAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCACGGTCTGCCTTATAATAGAAGAACAAAGc     <  1:1079623/76‑1 (MQ=255)                                                                        cAAGGAAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCTGCCTTATAATAGAAGAACAAAGc     <  2:682126/76‑1 (MQ=255)                                                                        cAAGGAAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCTGCCTTATAATAGAAGAACAAAGc     >  1:1528869/1‑76 (MQ=255)                                                                          aGG‑AAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCTGCCTTATAATAGAAGAACAAAGCAAg  >  2:27924/1‑76 (MQ=255)                                                                             |                                                                          TTACCTATATTTATCATTGGCTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGG‑AAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCTGCCTTATAATAGAAGAACAAAGCAAG  >  REL606/798191‑798339

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