Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	242,203	(G)8→7	coding (193/372 nt)	ECB_00212 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	242,196	0	G	.	93.8%	205.3 / 12.5	50	coding (186/372 nt)	ECB_00212	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (3/0);  new base . (24/21);  total (27/23)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

ATCAAAACCCAACTTACATGGAGTAATTAAAAAAGTTGATTTGAACACAAATGAAGCGGAAATT‑AAAATAAAAGAGGGGGGGGCAGTGAGTTTTATTTTTAATGAGGAAACATTGTATTTTACTGTTTTTGGTGGGGATTTTATTAT  >  REL606/242121‑242267                                                                             |                                                                        aTCAAAACCCAACTTACATGGAGTAATTAAAAAAGTTGATTTGAACACAAATGAAGCGGAAATT‑AAAATAAAagag                                                                         >  2:1028214/1‑76 (MQ=255)        cccAACTTACATGGAGTAATTAAAAAAGTTGATTTGAACACAAATGAAGCGGAAATT‑AAAATAAAAGAgggggggc                                                                  <  1:995125/76‑2 (MQ=255)          cAACTTACATGGAGTAATTAAAAAAGTTGATTTGAACACAAATGAAGCGGAAATT‑AAAATAAAAGAgggggggcag                                                                >  2:429571/1‑73 (MQ=255)            aCTTACATGGAATAATTAAAAAAGTTGATTTGAACACAAATGAAGCGGAAATT‑AAAATAAAAGA‑GGGGGGGCAGTg                                                             <  2:558202/76‑1 (MQ=255)                aCATGGAGTAATTAAAAAAGTTGATTTGAACACAAATGAAGCGGAAATT‑AAAATAAAAGA‑GGGGGGGCAGTGAGtt                                                         <  2:1085108/76‑1 (MQ=255)                 cATGGAGTAATTAAAAAAGTTGATTTGAACACAAATGAAGCGGAAATT‑AAAATAAAAGA‑GGGGGGGCAGTGAGttt                                                        >  2:465134/1‑76 (MQ=255)                  aTGGAGTAATTAAAAAAGTTGATTTGAACACAAATGAAGCGGAAATT‑AAAATAAAAGA‑GGGGGGGCAGTGAGtttt                                                       <  2:184414/76‑1 (MQ=255)                  aTGGAGTAATTAAAAAAGTTGATTTGAACACAAATGAAGCGGAAATT‑AAAATAAAAGA‑GGGGGGGCAGTGAGtttt                                                       <  2:197125/76‑1 (MQ=255)                    ggAGTAATTAAAAAAGTTGATTTGAACACAAATGAAGCGGAAATT‑AAAATAAAAGA‑GGGGGGGCAGTGAGTTTTAt                                                     <  2:349806/76‑1 (MQ=255)                      aGTAATTAAAAAAGTTGATTTGAACACAAATGAAGCGGAAATT‑AAAATAAAAGA‑GGGGGGGCAGTGAGTTTTAttt                                                   <  2:1670059/76‑1 (MQ=255)                      aGTAATTAAAAAAGTTGATTTGAACACAAATGAAGCGGAAATT‑AAAATAAAAGA‑GGGGGGGCAGTGAGTTTTAttt                                                   >  1:196172/1‑76 (MQ=255)                       gTAATTAAAAAAGTTGATTTGAACACAAATGAAGCGGAAATT‑AAAATAAAAGA‑GGGGGGGCAGTGAGTTTTAtttt                                                  <  2:649033/76‑1 (MQ=255)                        tAATTAAAAAAGTTGATTTGAACACAAATGAAGCGGAAATT‑AAAATAAAAGA‑GGGGGGGCAGTGAGTTTTAttttt                                                 >  1:999485/1‑76 (MQ=255)                         aaTTAAAAAAGTTGATTTGAACACAAATGAAGCGGAAATT‑AAAATAAAAGA‑GGGGGGGCAGTGAGTTTTATTTTTa                                                >  1:357667/1‑76 (MQ=255)                            tAAAAAAGTTGATTTGAACACAAATGAAGCGGAAATT‑AAAATAAAAGA‑GGGGGGGCAGTGAGTTTTATTTTTAATg                                             <  2:31281/76‑1 (MQ=255)                              aaaaaGTTTTTTTTAAAAAAAAAGAAAAGGAAAAAAAAAAAAAAAGGAGGGGGGGCAGTGAGTTTTATTTTTAATg                                             <  1:838897/76‑1 (MQ=255)                              aaaaaGTTGATTTGAACACAAATGAAGCGGAAATT‑AAAATAAAAGA‑GGGGGGGCAGTGAGTTTTATTTTTAATGAg                                           >  2:1198435/1‑76 (MQ=255)                              aaaaaGTTGATTTGAACAAAAATGAAGCGGAAATT‑AAAATAAAAGA‑GGGGGGGCAGTGAGTTTTATTTTTAATGAg                                           <  1:1476552/76‑1 (MQ=255)                               aaaaGTTGATTTGAACACAAATGAAGCGGAAATT‑AAAATAAAAGA‑GGGGGGGCAGTGAGTTTTATTTTTAATGAgg                                          >  2:1050445/1‑76 (MQ=255)                               aaaaGTTGATTT‑AAAAAAAAAGAAGGGGAAAAT‑AAAAAAAAAAGAGGGGGGGCAGTGAGTTTTATTTTTAATGAgg                                          <  1:1104049/76‑1 (MQ=255)                                  aGTTGATTTGAACACAAATGAAGCGGAAATT‑AAAATAAAAGAGGGGGGGGCAGTGAGTTTTATTTTTAATGAGGaa                                        >  2:39906/1‑76 (MQ=255)                                  aGTTGATTTGAACACAAATGAAGCGGAAATT‑AAAATAAAAGA‑GGGGGGGCAGTGAGTTTTATTTTTAATGAGGaaa                                       >  2:372387/1‑76 (MQ=255)                                  aGTTGATTTGAACACAAATGAAGCGGAAATT‑AAAATAAAAGA‑GGGGGGGCAGTGAGTTTTATTTTTAATGAGGaaa                                       <  1:1536007/76‑1 (MQ=255)                                    ttGATTTGAACACAAATGAAGCGGAAATT‑AAAATAAAAGA‑GGGGGGGCAGTGAGTTTTATTTTTAATGAGGAAACa                                     <  2:4034/76‑1 (MQ=255)                                    ttGATTTGAACACAAATGAAGCGGAAATT‑AAAATAAAAGA‑GGGGGGGCAGTGAGTTTTATTTTTAATGAGGAAACa                                     >  2:814833/1‑76 (MQ=255)                                     tGATTTGAACACAAATGAAGCGGAAATT‑AAAATAAAAGA‑GGGGGGGCAGTGAGTTTTATTTTTAATGAGGAAACAt                                    >  1:782381/1‑76 (MQ=255)                                     tGATTTGAACACAAATGAAGCGGAAATT‑AAAATAAAAGA‑GGGGGGGCAGTGAGTTTTATTTTTAATGAGGAAACAt                                    >  1:1176507/1‑76 (MQ=255)                                           gAACACAAATGAAGCGGAAATT‑AAAATAAAAGA‑GGGGGGGCAGTGAGTTTTATTTTTAATGAGGAAACATTGTAtt                              >  1:1196303/1‑76 (MQ=255)                                               acaAATGAAGCGGAAATT‑AAAATAAAAGA‑GGGGGGGCAGTGAGTTTTATCTTTANTGAGGAAACATTGTATTTTAc                          >  2:280163/1‑76 (MQ=255)                                                caAATGAAGCGGAAATT‑AAAATAAAAGA‑GGGGGGGCAGTGAGTTTTATTTTTAATGAGGAAACATTGTATTTTACt                         >  1:1084933/1‑76 (MQ=255)                                                 aaaTGAAGCGGAAATT‑AAAATAAAAGA‑GGGGGGGCAGTGAGTTTTATTTTTAATGAGGAAACATTGTATTTTACTg                        >  1:203945/1‑76 (MQ=255)                                                   aTGAAGCGGAAATT‑AAAATAAAAGA‑GGGGGGGCAGTGAGTTTTATTTTTAATGAGGAAACATTGTATTTTACTGtt                      <  1:1336503/76‑1 (MQ=255)                                                    tGAAGCGGAAATT‑AAAAAAAAAGA‑GGGGGGGCAGTGAGTTTTATTTTTAATGAGGAAACATTGTATTTTACTGttt                     <  2:999516/76‑1 (MQ=255)                                                     gAAGCGGAAATT‑AAAATAAAAGA‑GGGGGGGCAGTGAGTTTTATTTTTAATGAGGAAACATTGTATTTTACTGtttt                    <  1:191830/76‑1 (MQ=255)                                                     gAAGCGGAAATT‑AAAATAAAAGA‑GGGGGGGCAGTGAGTTTTATTTTTAATGAGGAAACATTGTATTTTACTGtttt                    >  2:459498/1‑76 (MQ=255)                                                     gAAGCGGAAATT‑AAAATAAAAGA‑GGGGGGGCAGTGAGTTTTATTTTTAATGAGGAAACATTGTATTTTACTGtttt                    >  2:421188/1‑76 (MQ=255)                                                     gAAGCGGAAATT‑AAAATAAAAGA‑GGGGGGGCAGTGAGTTTTATTTTTAATGAGGAAACATTGTATTTTACTGtttt                    >  2:1187387/1‑76 (MQ=255)                                                     gAAGCGGAAATT‑AAAATAAAAGA‑GGGGGGGCAGTGAGTTTTATTTTTAATGAGGAAACATTGTATTTTACTGtttt                    >  2:1134129/1‑76 (MQ=255)                                                     gAAGCGGAAATT‑AAAATAAAAGA‑GGGGGGGCAGTGAGTTTTATTTTCAATGAGGAAACATTGTATTTTACTGtttt                    >  2:291228/1‑76 (MQ=255)                                                      aaGCGGAAATT‑AAAATAAAAGA‑GGGGGGGCAGTGAGTTTTATTTTTAATGAGGAAACATTGTATTTTACTGttttt                   >  2:1574348/1‑76 (MQ=255)                                                          ggAAATT‑AAAATAAAAGA‑GGGGGGGCAGTGAGTTTTATTTTTAATGAGGAAACATTGTATTTTACTGTTTTtggtg               <  2:427891/76‑1 (MQ=255)                                                          ggAAATT‑AAAATAAAAGA‑GGGGGGGCAGTGAGTTTTATTTTTAATGAGGAAACATTGTATTTTACTGTTTTtggtg               >  1:1708643/1‑76 (MQ=255)                                                              aat‑aaaaaaaaaGGAGGGGGGGCAGTGAGTTTTATTTTTAATGAGGAAACATTGTATTTTACTGTTTTTGGTgggg            <  1:1080443/74‑1 (MQ=255)                                                              aTT‑AAAATAAAAGAGGGGGGGGCAGTGAGTTTTATTTTTAATGAGGAAACATTGTATTTTACTGTTTTTGGTgggg            >  2:1433432/1‑76 (MQ=255)                                                              aTT‑AAAATAAAAGAGGGGGGGGCAGTGAGTTTTATTTTTAATGAGGAAACATTGTATTTTACTGTTTTTGGTgggg            >  1:393228/1‑76 (MQ=255)                                                                t‑aaaataaaaGA‑GGGGGGGCAGTGAGTTTTATTTTTAATGAGGAAACATTGTATTTTACTGTTTTTGGTGGGGAtt         >  1:1182992/1‑76 (MQ=255)                                                                   aaataaaaGA‑GGGGGGGCAGTGAGTTTTATTTTTAATGAGGAAACATTGTATTTTACTGTTTTTGGTGGGGAtttt       <  1:1586280/76‑1 (MQ=255)                                                                   aaataaaaGA‑GGGGGGGCAGTGAGTTTTATTTTTAATGAGGAAACATTGTATTTTACTGTTTTTGGTGGGGAtttt       <  2:357667/76‑1 (MQ=255)                                                                   aaataaaaGA‑GGGGGGGCAGTGAGTTTTATTTTTAATGAGGAAACATTGTATTTTACTGTTTTTGGTGGGGAtttt       <  2:1302072/76‑1 (MQ=255)                                                                      taaaaGA‑GGGGGGGCAGTGAGTTTTATTTTTAATGAGGAAACATTGTATTTTACTGTTTTTGGTGGGGATTttatt    >  2:371587/1‑76 (MQ=255)                                                                      taaaaGA‑GGGGGGGCAGTGAGTTTTATTTTTAATGAGGAAACATTGTATTTTACTGTTTTTGGTGGGGATTttatt    <  2:348342/76‑1 (MQ=255)                                                                      aaaaaGGAGGGGGGGCAGTGAGTTTTATTTTTAATGAGGAAACATTGTATTTTACTGTTTTTGGTGGGGATTttat     <  1:1438915/75‑1 (MQ=255)                                                                       aaaaGA‑GGGGGGGCAGTGAGTTTTATTTTTAATGAGGAAACATTGTATTTTACTGTTTTTGGTGGGGATTttatta   <  2:504786/76‑1 (MQ=255)                                                                       aaaaGA‑GGGGGGGCAGTGAGTTTTATTTTTAATGAGGAAACATTGTATTTTACTGTTTTTGGTGGGGATTttatta   >  2:864346/1‑76 (MQ=255)                                                                        aaaGA‑GGGGGGGCAGTGAGTTTTATTTTTAATGAGGAAACATTGTATTTTACTGTTTTTGGTGGGGATTTtattat  <  1:215583/76‑1 (MQ=255)                                                                             |                                                                        ATCAAAACCCAACTTACATGGAGTAATTAAAAAAGTTGATTTGAACACAAATGAAGCGGAAATT‑AAAATAAAAGAGGGGGGGGCAGTGAGTTTTATTTTTAATGAGGAAACATTGTATTTTACTGTTTTTGGTGGGGATTTTATTAT  >  REL606/242121‑242267

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