Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,776,362	T→G	I235M (ATT→ATG)	ECB_03516 →	Int

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,776,362	0	T	G	96.0%	314.6 / 3.3	76	I235M (ATT→ATG)	ECB_03516	Int
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base G (37/36);  total (39/37)
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.38e-01

CTTCCCCACTTCCTCAACGATCTGGCGGGTTATACCGGAAGTATAATCACTAAGACAGCTACTCAGATCATTATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGCTGGGAAGATATCGATTTTGAGACAAAAT  >  REL606/3776291‑3776432                                                                        |                                                                       cTTCCCCACTTCCTCAACGATCTGGCGGGTTATACCGGAAGTATAATCACTAAGACAGCTACTCAGATCatgatg                                                                     <  1:80468/75‑1 (MQ=255)   tCCCCACTTCCTCAACGATCTGGCGGGTTATACCGGAAGTATAATCACTAAGACAGCTACTCAGATCATGATGtt                                                                   <  2:824469/75‑1 (MQ=255)   tCCCCACTTCCTCAACGATCTGGCGGGTTATACCGGAAGTATAATCACTAAGACAGCTACTCAGATCATGATGtt                                                                   >  2:1376761/1‑75 (MQ=255)   tCCCCACTTCCTCAACGATCTGGCGGGTTATACCGGAAGTATAATCACTAAGACAGCTACTCAGATCATGATGtt                                                                   <  2:1155037/75‑1 (MQ=255)    ccccACTTCCTCAACGATCTGGCGGGTTATACCGGAAGTATAATCACTAAGACAGCTACTCAGATCATGATGTTa                                                                  <  1:1791893/75‑1 (MQ=255)     cccACTTCCTCAACGATCTGGCGGGTTATACCGGAAGTATAATCACTAAGACAGCTACTCAGATCATGATGTTaa                                                                 >  2:1194702/1‑75 (MQ=255)     cccACTTCCTCAACGATCTGGCGGGTTATACCGGAAGTATAATCACTAAGACAGCTACTCAGATCATGATGTTaa                                                                 >  1:776427/1‑75 (MQ=255)       cACTTCCTCAACGATCTGGCGGGTTATACCGGAAGTATAATCACTAAGACAGCTACTCAGATCATGATGTTAAcc                                                               >  2:580099/1‑75 (MQ=255)         cTTCCTCAACGATCTGGCGGGTTATACCGGAAGTATAATCACTAAGACAGCTACTCAGATCATGATGTTAACCgg                                                             <  2:99713/75‑1 (MQ=255)           tCCTCAACGATCTGGCGGGTTATACCGGAAGTATAATCACTAAGACAGCTACTCAGATCATGATGTTAACCGGCg                                                           >  1:1193724/1‑75 (MQ=255)           tCCTCAACGATCTGGCGGGTTATACCGGAAGTATAATCACTAAGACAGCTACTCAGATCATGATGTTAACCGGCg                                                           >  2:779137/1‑75 (MQ=255)           tCCTCAACGATCTGGCGGGTTATACCGGAAGTATAATCACTAAGACAGCTACTCAGATCATGATGTTAACCGGCg                                                           <  1:154914/75‑1 (MQ=255)             cTCAACGATCTGGCGGGTTATACCGGAAGTATAATCACTAAGACAGCTACTCAGATCATGATGTTAACCGGCGTa                                                         >  2:1581768/1‑75 (MQ=255)              tCAACGATCTGGCGGGTTATACCGGAAGTATAATCACTAAGACAGCTACTCAGATCATGATGTTAACCGGCGTAc                                                        >  2:1835791/1‑75 (MQ=255)                aaCGATCTGGCGGGTTATACCGGAAGTATAATCACTAAGACAGCTACTCAGATCATTATGTTAACAGGCGTACGa                                                      >  2:51546/1‑75 (MQ=255)                aaCGATCTGGCGGGTTATACCGGAAGTATAATCACTAAGACAGCTACTCAGATCATGATGTTAACCGGCGTACGa                                                      >  1:772013/1‑75 (MQ=255)                  cGATCTGGCGGGTTATACCGGAAGTATAATCACTAAGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAac                                                    <  2:371600/75‑1 (MQ=255)                     tCTGGCGGGTTATACCGGAAGTATAATCACTAAGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAAcaca                                                 <  1:316859/75‑1 (MQ=255)                      cTGGCGGGTTATACCGGAAGTATAATCACTAAGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACaa                                                <  1:531437/75‑1 (MQ=255)                      cTGGCGGGTTATACCGGAAGTATAATCACTAAGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACaa                                                >  1:577846/1‑75 (MQ=255)                      cTGGCGGGTTATACCGGAAGTATAATCACTAAGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACaa                                                <  2:1573786/75‑1 (MQ=255)                        ggCGGGTTATACCGGAAGTATAATCACTAAGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGa                                              <  1:1561058/75‑1 (MQ=255)                          cGGATTATACCGGAAGTATAATCACTAAGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAAc                                            <  2:1124813/75‑1 (MQ=255)                            ggTTATACCGGAAGTATAATCACTAAGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTg                                          >  2:179638/1‑75 (MQ=255)                            ggTTATACCGGAAGTATAATCACTAAGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTg                                          >  2:1117481/1‑75 (MQ=255)                               tataCCGGAAGTATAATCACTAAGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGt                                       >  2:1604149/1‑75 (MQ=255)                                  aCCGGAAGTATAATCACTAAGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGtttt                                    >  1:1200989/1‑75 (MQ=255)                                   ccGGAAGTATAATCACTAAGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTg                                   >  2:579600/1‑75 (MQ=255)                                   ccGGAAGTATAATCACTAAGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTg                                   >  2:864693/1‑75 (MQ=255)                                   ccGGAAGTATAATCACTAAGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTg                                   <  1:235383/75‑1 (MQ=255)                                    cGGAAGTATAATCACTAAGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGc                                  <  2:23670/75‑1 (MQ=255)                                    cGGAAGTATAATCACTAAGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGc                                  <  2:741988/75‑1 (MQ=255)                                    cGGAAGTATAATCACTAAGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGc                                  <  2:767345/75‑1 (MQ=255)                                     ggAAGTATAATCACTAAGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCa                                 >  2:524240/1‑75 (MQ=255)                                      gAAGTATAATCACTAAGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCAc                                >  1:1620765/1‑75 (MQ=255)                                        aGTATAATCACTAAGACAGCTACTCAGATCATGATGTTAACCGGCGTACGATCACAAGAACTGCGTTTTGCACGc                              <  2:1424507/75‑1 (MQ=255)                                        aGTATAATCACTAAGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGc                              <  2:936852/75‑1 (MQ=255)                                        aGTATAATCACTAAGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGc                              >  2:1592444/1‑75 (MQ=255)                                        aGTATAATCACTAAGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGc                              <  1:32046/75‑1 (MQ=255)                                          tataATCACTAAGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGCTg                            >  2:1961835/1‑75 (MQ=255)                                           ataATCACTAAGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGCTgg                           >  2:1571382/1‑75 (MQ=255)                                            taATCACTAAGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGCTggg                          >  1:136445/1‑75 (MQ=255)                                            taATCACCAAGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGCTggg                          >  1:1685068/1‑75 (MQ=255)                                              aTCACTAAGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGCTGGGaa                        <  2:603426/75‑1 (MQ=255)                                              aTCACTAAGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGCTGGGaa                        <  1:193163/75‑1 (MQ=255)                                              aTCACTAAGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGCTGGGaa                        >  2:2120477/1‑75 (MQ=255)                                              aTCACTAAGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGCTGGGaa                        <  1:6451/75‑1 (MQ=255)                                                cACTAAGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGCTGGGAAGa                      <  1:499325/75‑1 (MQ=255)                                                  cTAAGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGCTGGGAAGata                    <  1:1885221/75‑1 (MQ=255)                                                   tAAGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGCTGGGAAGatat                   <  1:329688/75‑1 (MQ=255)                                                   tAAGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGCTGGGAAGatat                   <  2:1926530/75‑1 (MQ=255)                                                   tAAGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGCTGGGAAGatat                   <  1:1392006/75‑1 (MQ=255)                                                   tAAGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGCTGGGAAGatat                   <  1:879744/75‑1 (MQ=255)                                                   tAAGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGCTGGGAAGatat                   <  2:1148897/75‑1 (MQ=255)                                                    aaGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGCTGGGAAGATATc                  >  1:825257/1‑75 (MQ=255)                                                    aaGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGCTGGGAAGATATc                  >  2:2101026/1‑75 (MQ=255)                                                    aaGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGCTGGGAAGATATc                  >  1:1353565/1‑75 (MQ=255)                                                     aGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGCTGGGAAGATATCg                 >  1:1140031/1‑75 (MQ=255)                                                     aGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGCTGGGAAGATATCg                 >  1:96062/1‑75 (MQ=255)                                                     aGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGCTGGGAAGATATCg                 >  1:24099/1‑75 (MQ=255)                                                     aGACAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGCTGGGAAGATATCg                 >  1:1636829/1‑75 (MQ=255)                                                      gACAGCTACTCAGATCATTATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGCTGGGAAGATATCGa                <  1:1410658/75‑1 (MQ=255)                                                        cAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGCTGGGAAGATATCGAtt              <  2:1610565/75‑1 (MQ=255)                                                        cAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGCTGGGAAGAGATCGact              >  2:966741/1‑73 (MQ=255)                                                        cAGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGCGGGGGAGAGATCGAtt              >  2:966747/1‑75 (MQ=255)                                                         aGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGCTGGGAAGATATCGAttt             <  1:1373353/75‑1 (MQ=255)                                                         aGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGCTGGGAAGATATCGAttt             <  1:362324/75‑1 (MQ=255)                                                         aGCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGCTGGGAAGATATCGAttt             >  1:351476/1‑75 (MQ=255)                                                          gCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGCTGGGAAGATATCGAtttt            <  1:99093/75‑1 (MQ=255)                                                          gCTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGCTGGGAAGATATCGAtttt            >  1:1660086/1‑75 (MQ=255)                                                           cTACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGCTGGGAAGATATCGGTTTTg           >  1:312513/1‑75 (MQ=255)                                                            tACTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGCTGGGAAGATATCGATTTTga          <  1:295795/75‑1 (MQ=255)                                                             aCTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGCTGGGAAGATATCGATTTTgag         <  2:1450026/75‑1 (MQ=255)                                                             aCTCAGATCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGCTGGGAAGATATCGATTTTgag         <  2:94209/75‑1 (MQ=255)                                                                cAGATCATTATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGCTGGGAAGATATCGATTTTGAGACa      >  2:1267633/1‑75 (MQ=255)                                                                   aTCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGCTGGGAAGATATCGATTTTGAGACaaaa   >  2:577001/1‑75 (MQ=255)                                                                    tCATGATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGCTGGGAAGATATCGATTTTGAGACAAAAt  <  1:339297/75‑1 (MQ=255)                                                                        |                                                                       CTTCCCCACTTCCTCAACGATCTGGCGGGTTATACCGGAAGTATAATCACTAAGACAGCTACTCAGATCATTATGTTAACCGGCGTACGAACACAAGAACTGCGTTTTGCACGCTGGGAAGATATCGATTTTGAGACAAAAT  >  REL606/3776291‑3776432

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