Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,015,178	A→T	D200V (GAT→GTT)	ECB_02816 →	KpsD protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,015,178	0	A	T	98.7%	304.7 / ‑3.3	79	D200V (GAT→GTT)	ECB_02816	KpsD protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base T (40/38);  total (40/39)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.75e-01

TAATCCTGGTCTGTATGGCGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGATCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTAACCGCGTGCGCTCCAACGTCAACCTGTA  >  REL606/3015107‑3015250                                                                        |                                                                         tAATCCTGGTCTGTATGGCGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGTTCCa                                                                      <  1:1628460/76‑1 (MQ=255)  aaTCCTGGTCTGTATGGCGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGTTCCag                                                                     <  1:1109554/76‑1 (MQ=255)   aTCCTGGTCTGTATGGCGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGTTCCaga                                                                    >  1:1335835/1‑76 (MQ=255)     ccTGGTCTGTATGGCGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGTTCCAGAgc                                                                  <  1:1058572/76‑1 (MQ=255)      cTGGTCTGTATGGCGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGTTCCAGAgcg                                                                 <  2:954865/76‑1 (MQ=255)       tGGTCTGTATGGCGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGTTCCAGAgcgc                                                                >  2:2000849/1‑76 (MQ=255)        ggTCTGTATGGCGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGTTCCAGAGcgcg                                                               <  2:1588741/76‑1 (MQ=255)        ggTCTGTATGGCGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGTTCCAGAGcgcg                                                               <  2:801526/76‑1 (MQ=255)        ggTCTGTATGGCGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGTTCCAGAGcgcg                                                               >  2:709907/1‑76 (MQ=255)          tCTGTATGGCGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGa                                                             <  1:1260438/76‑1 (MQ=255)           cTGTATGGCGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGaa                                                            >  1:1127057/1‑76 (MQ=255)               aTGGCGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTa                                                        <  1:132452/76‑1 (MQ=255)                  gCGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGt                                                     <  1:1644194/76‑1 (MQ=255)                     gtgtGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGa                                                  >  1:1626542/1‑76 (MQ=255)                      tgtgACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGat                                                 >  2:653662/1‑76 (MQ=255)                          aCGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATtg                                             >  1:472742/1‑76 (MQ=255)                            gTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATtgtg                                           <  1:1166938/76‑1 (MQ=255)                            gTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATtgtg                                           >  1:2139845/1‑76 (MQ=255)                             tCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATtgtag                                          >  2:1318079/1‑74 (MQ=255)                              cTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGt                                         >  2:1613659/1‑76 (MQ=255)                                gATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCa                                       <  1:986265/76‑1 (MQ=255)                                 aTTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCaa                                      <  1:1824133/76‑1 (MQ=255)                                 aTTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCaa                                      >  1:637831/1‑76 (MQ=255)                                 aTTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCaa                                      >  1:245554/1‑76 (MQ=255)                                  ttCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAg                                     >  1:1273303/1‑76 (MQ=255)                                   tCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAgc                                    <  1:1893808/76‑1 (MQ=255)                                    cGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAgcg                                   >  1:1495718/1‑76 (MQ=255)                                    cGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAgcg                                   <  1:3276/76‑1 (MQ=255)                                    cGTCACTCAATTATCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAgcg                                   >  2:798680/1‑76 (MQ=255)                                      ttACTCAATTATCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGcgcg                                 >  1:1955751/1‑76 (MQ=255)                                       tACTCAATTATCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCgg                                >  1:298456/1‑76 (MQ=255)                                        aCTCAATTATCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGt                               >  2:1570246/1‑76 (MQ=255)                                        aCTCAATTATCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGt                               >  1:1545549/1‑76 (MQ=255)                                        aCTCAATTATCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGt                               >  2:1325983/1‑76 (MQ=255)                                        aCTCAATTATCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGt                               >  2:2169370/1‑76 (MQ=255)                                         cTCAATTATCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTa                              <  1:841756/76‑1 (MQ=255)                                         cTCAATTATCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTa                              <  2:2191146/76‑1 (MQ=255)                                           cAATTATCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTAAc                            <  1:122773/76‑1 (MQ=255)                                           cAATTATCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTAAc                            <  1:736427/76‑1 (MQ=255)                                           cAATTATCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTAAc                            <  1:2089954/76‑1 (MQ=255)                                             aTTATCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTAACcg                          >  1:448701/1‑76 (MQ=255)                                             aTTATCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGGGGTCAAGCGCGGTACCcg                          >  2:1090550/1‑76 (MQ=255)                                              ttATCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTAACcgc                         <  2:1785686/76‑1 (MQ=255)                                               tATCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTAACcgcg                        <  1:482014/76‑1 (MQ=255)                                               tATCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTAACcgcg                        <  1:445485/76‑1 (MQ=255)                                               tATCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTAACcgcg                        >  2:1545121/1‑76 (MQ=255)                                               tATCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTAACcgcg                        <  2:2205528/76‑1 (MQ=255)                                                aTCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTAACcgcgg                       >  2:30309/1‑75 (MQ=255)                                                aTCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTAACCACGt                       >  1:828687/1‑76 (MQ=255)                                                aTCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGGAACcgggc                       >  1:510150/1‑73 (MQ=255)                                                 tCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTAACCGCGTg                      >  1:931442/1‑76 (MQ=255)                                                 tCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTAACCGCGTg                      <  2:1761334/76‑1 (MQ=255)                                                 tCTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGGAACCGCGTg                      >  2:1576369/1‑76 (MQ=255)                                                 tCTGATCAAGGCTGGCGGCGTTGATCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTAACCGCGTg                      <  1:31313/76‑1 (MQ=255)                                                  cTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTAACCGCGTgc                     >  2:88606/1‑76 (MQ=255)                                                  cTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTAACCGCGTgc                     <  2:1520478/76‑1 (MQ=255)                                                  cTGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTAACCGCGTgc                     <  2:1452782/76‑1 (MQ=255)                                                   tGATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTAACCGCGTgcg                    >  1:990474/1‑76 (MQ=255)                                                    gATCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGCTACGTTGATATTGTGGGCAAGCGCGGTAACCGCGTgcgc                   >  1:1471337/1‑76 (MQ=255)                                                      tCAAGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTAACCGCGTGCGCTc                 <  1:2076498/76‑1 (MQ=255)                                                        aaGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTAACCGCGTGCGCTCCa               >  2:2066734/1‑76 (MQ=255)                                                        aaGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTAACCGCGTGCGCTCCa               <  1:1082103/76‑1 (MQ=255)                                                        aaGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTAACCGCGTGCGCTCCa               <  2:1483558/76‑1 (MQ=255)                                                        aaGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTAACCGCGTGCGCTCCa               >  1:156795/1‑76 (MQ=255)                                                        aaGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTAACCGCGTGCGCTCCa               <  1:1686367/76‑1 (MQ=255)                                                        aaGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTAACCGCGTGCGCTCCa               <  1:1686373/76‑1 (MQ=255)                                                        aaGGCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTAACCGCGTGCGCTCCa               >  1:37097/1‑76 (MQ=255)                                                          ggCTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTAACCGCGTGCGCTCCAAc             <  1:2105467/76‑1 (MQ=255)                                                            cTGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTAACCGCGTGCGCTCCAACGt           <  2:1532593/76‑1 (MQ=255)                                                             tGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTAACCGCGTGCGCTCCAACGTc          <  2:1417054/76‑1 (MQ=255)                                                             tGGCGGCGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTAACCGCGTGCGCTCCAACGTc          >  1:1358829/1‑76 (MQ=255)                                                              ggcggcGTTGTTCCAGAGCGCGGAAGTTATGTTGATATTGTGGTCAAGCGCGGTAACCGCGTGCGCTCCAACGTCa         <  2:371602/76‑1 (MQ=255)                                                              ggcggcGTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTAACCGCGTGCGCTCCAACGTCa         <  2:2174832/76‑1 (MQ=255)                                                                cggcgTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTAACCGCGTGCGCTCCAACGTCAAc       <  1:514696/76‑1 (MQ=255)                                                                cggcgTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTAACCGCGTGCGCTCCAACGTCAAc       >  2:688016/1‑76 (MQ=255)                                                                cggcgTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTAACCGCGTGCGCTCCAACGTCAAc       >  1:899702/1‑76 (MQ=255)                                                                 ggcgTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTAACCGCGTGCGCTCCAACGTCAAcc      <  2:505709/76‑1 (MQ=255)                                                                  gcgTTGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTAACCGCGTGCGCTCCAACGTCAACCt     >  1:780566/1‑76 (MQ=255)                                                                    gttgttCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTAACCGCGTGCGCTCCAACGTCAACCTGt   >  1:1108433/1‑76 (MQ=255)                                                                     ttGTTCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTAACCGCGTGCGCTCCAACGTCAACCTGTa  >  1:684907/1‑76 (MQ=255)                                                                        |                                                                         TAATCCTGGTCTGTATGGCGGTGTGACGTCTGATTCGTTACTCAATTATCTGATCAAGGCTGGCGGCGTTGATCCAGAGCGCGGAAGTTACGTTGATATTGTGGTCAAGCGCGGTAACCGCGTGCGCTCCAACGTCAACCTGTA  >  REL606/3015107‑3015250

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