Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,530,076	(A)7→8	coding (155/483 nt)	ydcD →	putative immunity protein for RhsE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,530,069	1	.	A	96.5%	197.4 / ‑0.3	57	coding (148/483 nt)	ydcD	putative immunity protein for RhsE
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base A (33/22);  total (35/22)
Fisher's exact test for biased strand distribution p-value = 5.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑A‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGC  >  NC_000913/1529970‑1530162                                                                                                     |                                                                                                tGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTg                                                                                                  <  2:320113/100‑1 (MQ=255) tGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTg                                                                                                  <  2:209523/100‑1 (MQ=255)  gTTTATTTGTTGTTTACCTTTGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTg                                                                                                  <  2:993973/99‑1 (MQ=255)  gTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑a                                                                                                >  1:297562/1‑100 (MQ=255)   tttATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑a‑a                                                                                              >  1:583847/1‑100 (MQ=255)    ttATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑a‑aa                                                                                             <  1:397784/100‑1 (MQ=255)         tgttgtTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑a‑aaaaa‑aa                                                                                       <  2:857428/100‑2 (MQ=255)             gttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑a‑aaaaa‑aatatt                                                                                   >  1:91397/1‑95 (MQ=255)              tttACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAA‑AAAAA‑Atatta                                                                                   >  2:779724/1‑100 (MQ=255)                tACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGAACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAA‑AAAAA‑ATattatt                                                                                 <  2:975421/100‑1 (MQ=255)                   cTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAA‑AAAAA‑ATAttatttat                                                                              >  2:99561/1‑100 (MQ=255)                   cTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAA‑AAAAA‑ATAttatttat                                                                              >  1:75817/1‑100 (MQ=255)                   cTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAA‑AAAAA‑ATAttatttat                                                                              >  2:592151/1‑100 (MQ=255)                   cTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAA‑AAAAA‑ATAttatttat                                                                              >  1:285720/1‑100 (MQ=255)                    ttGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAA‑AAAAA‑ATATTATTtata                                                                             >  1:654558/1‑100 (MQ=255)                         tttGATAACAATTCAAATTCTGACATCGTATTTTATTGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAA‑AAAAA‑ATATTATTTATAAaata                                                                        >  1:1334391/1‑100 (MQ=255)                            gATAACAATTCAAATTCTGACATCTTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAA‑AAAAA‑ATATTATTTATAAAAtaatt                                                                     <  2:1016019/100‑1 (MQ=255)                                 cAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGa                                                                >  2:364453/1‑100 (MQ=255)                                  aaTTTAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGAc                                                               >  2:402432/1‑100 (MQ=255)                                   aTTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGAcc                                                              >  1:875664/1‑100 (MQ=255)                                    ttCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCa                                                             >  1:602199/1‑100 (MQ=255)                                      cAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCaaa                                                           >  1:652645/1‑100 (MQ=255)                                          ttCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAAt                                                       <  1:779728/100‑1 (MQ=255)                                            cTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATta                                                     >  1:717957/1‑100 (MQ=255)                                             tGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATtag                                                    <  1:740807/100‑1 (MQ=255)                                              gACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATtag                                                    <  2:717954/99‑1 (MQ=255)                                               aaatagTATTTTATGGACATAAATCGCCAAAGAGCCTTGAAATAAAAAAAAAAAAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATtag                                                    <  1:1361372/93‑1 (MQ=255)                                               aCATTGTATTTTATGGACATAAAACACTAAAGAGCGTTGAGATATATCTTTCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATtag                                                    <  2:280257/98‑1 (MQ=255)                                                 aTCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATTAGTaga                                                >  1:936868/1‑100 (MQ=255)                                                  tCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATTAGTagag                                               <  2:875657/100‑1 (MQ=255)                                                  tCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATTAGTagag                                               >  2:449373/1‑100 (MQ=255)                                                       ttttATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAAt                                          <  1:521667/100‑1 (MQ=255)                                                         ttATAGACATAAAACACCAAAGAGCGTTAAGATATATCTTTATA‑A‑AAAAACAGATTATTTATAAAATACTTATTGACCaa                                                            >  2:1277269/1‑80 (MQ=255)                                                             ggACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAt                                    <  1:605583/100‑1 (MQ=255)                                                                 aTAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTta                                >  2:622187/1‑100 (MQ=255)                                                                  tAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtat                               >  2:640693/1‑100 (MQ=255)                                                                  tAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtat                               >  1:327767/1‑100 (MQ=255)                                                                  tAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtat                               >  2:670974/1‑100 (MQ=255)                                                                  tAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtat                               >  2:296166/1‑100 (MQ=255)                                                                  tAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtat                               >  2:641111/1‑100 (MQ=255)                                                                  tAAAACACCAAAGAGCGTTGAGATAAATAAAAAAAAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtat                               <  2:15590/100‑1 (MQ=255)                                                                   aaaaCACCAAAGAGCGTTGAGATATATCTTTCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtata                              <  2:576480/100‑1 (MQ=255)                                                                   aaaaCACCAAAGAGCGTTGAGATATATCTTTCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtata                              >  2:201579/1‑100 (MQ=255)                                                                     aaCACCAAAGAGCGTTGAGATATATCTTTCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                            <  2:285719/100‑1 (MQ=255)                                                                     aaCACCAAAGAGCGTTAAGATATATATTTAAAAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                            <  2:583844/100‑1 (MQ=255)                                                                      acacCAAAGAGCGTTGAGATATATCTTTCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGt                           >  1:115244/1‑100 (MQ=255)                                                                        acCAAAGAGCGTTGAGATATATCTTTCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGtat                         <  1:973774/100‑1 (MQ=255)                                                                        acCAAAAAGCGTTGAGATATAACTTTCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGta                          <  2:483746/99‑1 (MQ=255)                                                                         ccAAAGAGCGTTGAGATATATCTTTCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGtata                        >  1:544823/1‑100 (MQ=255)                                                                         ccAAAGAGCGTTGAGATATATCTTTCTG‑ACAAAAA‑ATATTATTTACAAAATAAGTAATGACCAAAAAATTAATAGAGGGAATGGTCATTTTATAAGtata                        >  2:1296507/1‑100 (MQ=255)                                                                          cAAAGAGCGTTGAGATATATCTTTCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATaa                       >  1:502584/1‑100 (MQ=255)                                                                               agCGTTGAGATATATCTTTCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGtt                  >  2:19847/1‑100 (MQ=255)                                                                               agCGTTGAAATATATATTTCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGtt                  <  2:1403114/100‑1 (MQ=255)                                                                                 cGTTGAGATATATCTTTCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaa                >  1:1470486/1‑100 (MQ=255)                                                                                    tGAGATATATCTTTCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaataa             <  1:830744/100‑1 (MQ=255)                                                                                     gagaTATATCTTTCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGACTGGTCATTTTATAAGTATAATGGTTaataat            >  1:1099235/1‑100 (MQ=255)                                                                                       gaTATATCTTTCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTa          >  1:1995/1‑100 (MQ=255)                                                                                         tatatCTTTCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGATGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACa        <  1:402433/100‑1 (MQ=255)                                                                                         tatatCTTTCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACa        <  2:652642/100‑1 (MQ=255)                                                                                          atatCTTTCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAg       <  2:1177859/100‑1 (MQ=255)                                                                                             cttttCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGAc    <  1:675835/98‑1 (MQ=255)                                                                                              cTTTCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACg   >  1:573485/1‑100 (MQ=255)                                                                                              cTTTCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACg   >  1:590653/1‑100 (MQ=255)                                                                                               tttCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGAcag  >  2:675832/1‑98 (MQ=255)                                                                                               tttCTGAA‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATAAGTAGAGGTAATGATCATTTTATAAGTATAATGGTTAATAATTACAGGa     >  2:1422120/1‑97 (MQ=255)                                                                                                     |                                                                                                TGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑A‑AAAAA‑ATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGC  >  NC_000913/1529970‑1530162

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