Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	798,255	T→C	I341V (ATC→GTC)	ECB_00739 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	798,255	0	T	C	94.7%	139.0 / ‑1.3	38	I341V (ATC→GTC)	ECB_00739	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (20/16);  total (21/17)
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.54e-01

TTAAACGTTTACCTATATTTATCATTGGCTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGGAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCTGCCTTATAA  >  REL606/798183‑798322                                                                         |                                                                    ttAAACGTTTACCTATATTTATCATTGGCTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGacaat                                                                  <  2:789722/76‑1 (MQ=255) ttAAACGTTTACCTATATTTATCATTGGCTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGacaat                                                                  <  1:250222/76‑1 (MQ=255) ttAAACGTTTACCTATATTTATCATTGGCTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGacaat                                                                  <  1:298086/76‑1 (MQ=255)  tAAACGTTTACCTATATTTATCATTGGCTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGacaata                                                                 >  1:941596/1‑76 (MQ=255)   aaaCGTTTACCTATATTTATCATTGGCTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGAcaataa                                                                >  2:782027/1‑76 (MQ=255)   aaaCGTTTACCTATATTTATCATTGGCTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGAcaataa                                                                <  1:262174/76‑1 (MQ=255)    aaCGTTTACCTATATTTANCATTGGCTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGACAATAAc                                                               >  1:1113145/1‑76 (MQ=255)     aCGTTTACCTATATTTATCATTGGCTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGACAATAAcc                                                              <  1:592063/76‑1 (MQ=255)        tttACCTATATTTATCATTGGCTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGACAATAACCAAg                                                           >  2:888277/1‑76 (MQ=255)          tACCTATATTTATCATTGGCTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGACAATAACCAAGGa                                                         <  1:203014/76‑1 (MQ=255)             cTATATTTATCATTGGCTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGACAATAACCAAGGaaaa                                                      <  2:75379/76‑1 (MQ=255)              tataTTTATCATTGGCTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGACAATAACCAAGGaaaaa                                                     <  1:362484/76‑1 (MQ=255)               atatTTATCATTGGCTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGCaaaaaa                                                    >  2:1133284/1‑76 (MQ=255)                    tATCATTGGCTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGACAATAACCAAGGAAAAAAAcaca                                               >  2:410605/1‑76 (MQ=255)                       cATTGGCTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGACAATAACCAAGGAAAAAAACACAGaa                                            <  1:876386/76‑1 (MQ=255)                       cATTGGCTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGACAATAACCAAGGAAAAAAACACAGaa                                            >  1:644948/1‑76 (MQ=255)                        aTTGGCTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGACAATAACCAAGGAAAAAAACACAGaaa                                           >  2:265670/1‑76 (MQ=255)                         ttGGCTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGACAATAACCAAGGAAAAAAACACAGaaaa                                          <  1:507866/76‑1 (MQ=255)                         ttGGCTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGACAATAACCAAGGAAAAAAACACAGaaaa                                          >  1:783549/1‑76 (MQ=255)                             cTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGACAATAACCAAGGAAAAAAACACAGAAAATcca                                      >  2:839935/1‑76 (MQ=255)                                ttCGATGAGAATATATATTGGATATGATGCAGTTATTGAGACAATAACCAAGGAAAAAAACACAGAAAATccacca                                   >  1:401126/1‑76 (MQ=255)                                ttCAATGAGAATATATATTGGATATGATGCAGTTATTGAGACAATAACCAAGGAAAAAAACACAGAAAATccacca                                   <  2:194151/76‑1 (MQ=255)                                  cAATGAGAATATATATTGGATATGATGCAGTTATTGAGACAATAACCAAGGAAAAAAACACAGAAAATCCACCAGt                                 >  2:475859/1‑76 (MQ=255)                                  cAATGAGAATATATATTGGATATGATGCAGTTATTGAGACAATAACCAAGGAAAAAAACACAGAAAATCCACCAGt                                 >  2:475855/1‑76 (MQ=255)                                   aaTGAGAATATATATTGGATATGATGCANTTATTGAGACGATAACCAGGGAAAAAAACACAGAAAATCCACCAGtt                                <  1:106584/76‑1 (MQ=255)                                    aTGAGAATATATATTGGATATGATGCAGTTATTGAGACAATAACCAAGNAAAAAAACACAGAAAATCCACCAGTTa                               <  1:1160377/76‑1 (MQ=255)                                    aTGAGAATATATATTGGATATGATGCAGTTATTGAGACAATAACCAAGGAAAAAAACACAGAAAATCCACCAGTTa                               >  2:59530/1‑76 (MQ=255)                                        gaATATATATTGGATATGATGCAGTTATTGAGACAATAACCAAGGAAAAAAACACAGAAAATCCACCAGTTATGcc                           >  1:672214/1‑76 (MQ=255)                                              atatTGGATATGATGCAGTTATTGAGACAATAACCACGGAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTc                     >  2:611937/1‑76 (MQ=255)                                                atTGGATATGATGCAGTTATTGAGACAATAACCAAGGAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAg                   <  1:365914/76‑1 (MQ=255)                                                atTGGATATGATGCAGTTATTGAGACAATAACCAAGGAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAg                   <  1:365896/76‑1 (MQ=255)                                                  tGGATATGATGCAGTTATTGAGACAATAACCAAGGAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCa                 <  1:252692/76‑1 (MQ=255)                                                   ggATATGATGCAGTTATTGAGACAATAACCAAGGAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCaa                >  2:21539/1‑76 (MQ=255)                                                    gATATGATGCAGTTATTGAGACAATAACCAAGGAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAg               >  1:115563/1‑76 (MQ=255)                                                      tatGATGCAGTTATTGAGACAATAACCAAGGAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGt             <  2:353441/76‑1 (MQ=255)                                                      tatGATGCAGTTATTGAGACAATAACCAAGGAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGt             >  2:424936/1‑76 (MQ=255)                                                      aatgatgCCATTATTGAGACATTACNCAAGGAAAAACACANAGAAACTGCANCNGTTATGCGAATTTCAGCAAGGt             >  2:20098/2‑76 (MQ=255)                                                        tgatgCAGTTATTGAGACAATAACCAAGGAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCt           <  1:1009125/76‑1 (MQ=255)                                                            gCAGTTATTGAGATAATAACCAAGGAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCTGCCt       <  1:703911/76‑1 (MQ=255)                                                            gCAGTTATTGAGACAATAACCAAGGAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCTGCCt       >  1:984684/1‑76 (MQ=255)                                                             cAGTTATTGAGACAATAACCAAGGAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCTGCCtt      <  2:768340/76‑1 (MQ=255)                                                               gTTATTGAGACAATAACCAAGGAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCTGCCTtat    >  2:540402/1‑76 (MQ=255)                                                               gTTATTGAGACAATAACCAAGGAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCTGCCTtat    >  2:540413/1‑76 (MQ=255)                                                               gTTATTCAGACAATAACCAAGGAAAAAAACACAGAAAATCCACCATTTATGCCAATTTCAGCAAGGTCTGCCTtat    <  1:1066559/76‑1 (MQ=255)                                                                 tATTGAGACAATAACCAAGGAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCTGCCTTataa  >  1:530935/1‑76 (MQ=255)                                                                         |                                                                    TTAAACGTTTACCTATATTTATCATTGGCTTTTCAATGAGAATATATATTGGATATGATGCAGTTATTGAGATAATAACCAAGGAAAAAAACACAGAAAATCCACCAGTTATGCCAATTTCAGCAAGGTCTGCCTTATAA  >  REL606/798183‑798322

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