Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	2,114,852	(T)8→7	coding (187/366 nt)	ECB_02002 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	2,114,845	0	T	.	98.0%	241.1 / ‑1.6	50	coding (194/366 nt)	ECB_02002	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (26/23);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

TTAAAATATCCATCAGTCAATGAGCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTATTTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTTATAAGATATTTCACTTTGGAATTGTCG  >  REL606/2114770‑2114916                                                                            |                                                                        ttAAAATATCCATCAGTCAATGAGCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCattat                                                                         <  1:1653700/76‑1 (MQ=255)  tAAAATATCCATCAGTCAATGAGCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCattatt                                                                        >  2:392442/1‑76 (MQ=255)    naaTATCCATCAGTCAATGAGCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTAtttt                                                                      >  1:1556346/2‑76 (MQ=255)     aaTATCCATCAGTCAATGCGCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCCTTAttttt                                                                     >  2:914935/1‑76 (MQ=255)      atatCCATCAGTCAATGAGCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTAtttttt                                                                    >  2:490783/1‑76 (MQ=255)        atCCATCAGTCAATGAGCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTAtttttttc                                                                  >  1:1851071/1‑75 (MQ=255)         tCCATCAGTCAATGAGCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTAtttttttct                                                                 >  2:1504683/1‑74 (MQ=255)          ccATCAGTCAATGAGCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTAtttttttctc                                                                >  1:1366684/1‑73 (MQ=255)          ccATCAGTCAATGAGCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTAtttttttctc                                                                <  1:1794994/76‑4 (MQ=255)            aTCAGTCAATGAGCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAg                                                             >  1:1393803/1‑76 (MQ=255)            aTCAGTCAATGAGCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAg                                                             >  2:1580385/1‑76 (MQ=255)             tCAGTCAATGAGCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGc                                                            <  2:1760692/76‑1 (MQ=255)                gTCAATGAGCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTga                                                         <  1:517483/76‑1 (MQ=255)                 tCAATGAGCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTgag                                                        <  1:867525/76‑1 (MQ=255)                 tCAATGAGCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTgag                                                        <  1:161881/76‑1 (MQ=255)                     tGAGCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCa                                                    >  1:407237/1‑76 (MQ=255)                         cTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATcc                                                <  2:325993/76‑1 (MQ=255)                              aaaTTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCa                                           >  1:1276012/1‑76 (MQ=255)                               aaTTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCaa                                          >  1:460659/1‑76 (MQ=255)                               aaTTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCaa                                          >  1:1547120/1‑76 (MQ=255)                                  tCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGa                                       >  2:945147/1‑76 (MQ=255)                                     ttACCTTTTTCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACAt                                    >  1:268845/1‑76 (MQ=255)                                     ttACCTTTTTCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACAt                                    >  1:1527729/1‑76 (MQ=255)                                      tACCTTTTTCAGATAAAAATAATTTAACACTCTCATTATTTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACAt                                    <  1:512748/76‑1 (MQ=255)                                      tACCTTTTTCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATa                                   >  2:1542349/1‑76 (MQ=255)                                       aCCTTTTTCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATaa                                  <  2:906290/76‑1 (MQ=255)                                       aCCTTTTTCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATaa                                  <  2:1725524/76‑1 (MQ=255)                                        ccTTTTTCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAAt                                 <  1:512249/76‑1 (MQ=255)                                           ttttCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAAtttt                              >  2:963604/1‑76 (MQ=255)                                           ttttCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAAtttt                              >  1:1329321/1‑76 (MQ=255)                                           ttttCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAAtttt                              >  1:280836/1‑76 (MQ=255)                                            tttCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAAttttt                             <  2:881399/76‑1 (MQ=255)                                            tttCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAAttttt                             <  2:1657542/76‑1 (MQ=255)                                            tttCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAAttttt                             <  2:1657520/76‑1 (MQ=255)                                            tttCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAAttttt                             >  1:1405340/1‑76 (MQ=255)                                             ttCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTta                            >  1:545565/1‑76 (MQ=255)                                               cAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACCTAATTTTtata                          >  2:1580944/1‑76 (MQ=255)                                               cAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTtata                          <  2:1510364/76‑1 (MQ=255)                                                aGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTTATaa                         >  1:749596/1‑76 (MQ=255)                                                    aaaaaTAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTTATAAGata                     <  2:1688941/76‑1 (MQ=255)                                                       aataatTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTTATAAGATAttt                  <  2:1762291/76‑1 (MQ=255)                                                       aataatTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTTATAAGATAttt                  <  2:1567806/76‑1 (MQ=255)                                                       aataatTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTTATAAGATAttt                  >  2:1246336/1‑76 (MQ=255)                                                        ataatTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTTATAAGATATTTc                 <  1:1037043/76‑1 (MQ=255)                                                         taatTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTTATAAGATATTTCa                <  2:1672975/76‑1 (MQ=255)                                                          aatTTAACAATCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTTATAAGATATTTCAc               <  1:668416/76‑1 (MQ=255)                                                           atTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTTATAAGATATTTCACt              <  2:1366064/76‑1 (MQ=255)                                                           atTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTTATAAGATATTTCACt              <  2:1366056/76‑1 (MQ=255)                                                             ttAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTTATAAGATATTTCACttt            >  1:955297/1‑76 (MQ=255)                                                             ttAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTTATAAGATATTTCACttt            >  1:63688/1‑76 (MQ=255)                                                              tAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTTATAAGATATTTCACTTTg           >  1:199101/1‑76 (MQ=255)                                                               aaCACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTTATAAGATATTTCACTTTgg          >  2:653233/1‑76 (MQ=255)                                                               aaCACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTTATAAGATATTTCACTTTgg          <  1:374628/76‑1 (MQ=255)                                                                acacTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTTATAAGATATTTCACTTTGGa         >  2:1071687/1‑76 (MQ=255)                                                                acacTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTTATAAGATATTTCACTTTGGa         <  1:854096/76‑1 (MQ=255)                                                                 cacTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTTATAAGATATTTCACTTTGGaa        >  2:1361802/1‑76 (MQ=255)                                                                   ctctCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTTATAAGATATTTCACTTTGGANtt      <  2:1322213/76‑1 (MQ=255)                                                                   ctctCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTTATAAGATATTTCACTTTGGAAtt      >  1:1180288/1‑76 (MQ=255)                                                                       cATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTTATAAGATATTTCACTTTGGAATTGTCg  >  2:760157/1‑76 (MQ=255)                                                                            |                                                                        TTAAAATATCCATCAGTCAATGAGCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTATTTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTTATAAGATATTTCACTTTGGAATTGTCG  >  REL606/2114770‑2114916

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