Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,008,404	+G	coding (487/516 nt)	b4555 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,008,402	1	.	G	97.5%	118.3 / ‑4.5	40	L162R (CTG→CGG)	b4555	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base G (17/22);  minor base A (1/0);  total (18/22)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

TCTGGAAGTGATGGCGGGCGGCGTTTC‑AAAGGCCATGCGCTGGAAG‑GCTGGCG‑‑‑‑AACT‑GGCTACAG‑CTGAAGGATTG‑ATTGCG‑TT‑‑G‑TGACGGGATGAACGACGCCGAAATGCT  >  NC_000913/4008346‑4008457                                                                |                                                              tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑tt                                 <  2:1756261/93‑1 (MQ=255)  gtGGAAGTGATGGCGGGCGGCGTGTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑ttt                                <  2:1612774/92‑2 (MQ=255)  cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑ttt                                <  2:2585349/93‑2 (MQ=255)   tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑tttg                               >  1:3760002/1‑91 (MQ=255)    ggAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑tttgg                              >  1:1504553/1‑90 (MQ=255)    ggAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑tttgg                              >  1:4328546/1‑90 (MQ=255)     gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑tttggt                             <  2:3225814/93‑5 (MQ=255)     gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑tttggt                             >  2:362797/1‑89 (MQ=255)           gtGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACggg                      <  1:3290333/92‑1 (MQ=255)           aTGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACggg                      <  1:980410/93‑1 (MQ=255)            tGGCGGGCGGCGTTTCAAAAGGCAATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGa                     <  2:2605564/93‑1 (MQ=255)             ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGAt                    <  2:1250992/93‑1 (MQ=255)              gcgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATg                   >  2:4734727/1‑93 (MQ=255)               cgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGa                  <  2:2641560/93‑1 (MQ=255)               cgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGa                  <  2:1009336/93‑1 (MQ=255)                gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGaa                 <  1:1004691/93‑1 (MQ=255)                gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGaa                 <  2:594135/93‑1 (MQ=255)                 ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAac                >  2:2025077/1‑93 (MQ=255)                 ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTAGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAac                >  2:1242468/1‑93 (MQ=255)                  gcggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAacg               <  1:3802206/93‑1 (MQ=255)                   cggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAacga              <  1:2414018/93‑1 (MQ=255)                    ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAacgac             <  1:321200/93‑1 (MQ=255)                    ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAacgac             <  2:907022/93‑1 (MQ=255)                     gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAacgacg            >  1:978131/1‑93 (MQ=255)                     gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAacgacg            >  1:1848270/1‑93 (MQ=255)                      cgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAACGACGc           <  2:1666186/93‑1 (MQ=255)                       gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAACGACGcc          <  2:4673438/93‑1 (MQ=255)                       gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAACGACGcc          >  1:1095282/1‑93 (MQ=255)                          tcaaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAACGACGCCGaa       >  2:3850580/4‑93 (MQ=255)                          gcaaaaGGCCCTGCGCTGTAAGCGCTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAACGACGCCGaa       <  2:3861309/90‑1 (MQ=255)                           caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAACGACGCCGaaa      >  1:2116878/3‑93 (MQ=255)                           caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAACGACGCCGaaa      <  2:643823/91‑1 (MQ=255)                            aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTT‑‑GTTGACGGGAGGAACTACGCCGAAAt     >  2:3805337/2‑93 (MQ=255)                            aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAACGACGCCGAAAt     >  2:2436174/2‑93 (MQ=255)                            aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAACGACGCCGAAAt     >  2:208377/2‑93 (MQ=255)                            aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAACGACGCCGAAAt     <  2:2079373/92‑1 (MQ=255)                            aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAACGACGCCGAAAt     <  1:3137890/92‑1 (MQ=255)                             aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAACGACGCCGAAATg    >  1:3644293/1‑93 (MQ=255)                             aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAACGACGCCGAAATg    >  2:991364/1‑93 (MQ=255)                               aGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAACGACGCCGAAATGCt  >  1:1869543/1‑93 (MQ=255)                                                                |                                                              TCTGGAAGTGATGGCGGGCGGCGTTTC‑AAAGGCCATGCGCTGGAAG‑GCTGGCG‑‑‑‑AACT‑GGCTACAG‑CTGAAGGATTG‑ATTGCG‑TT‑‑G‑TGACGGGATGAACGACGCCGAAATGCT  >  NC_000913/4008346‑4008457

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