Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,067,300	C→T	51.6%	G8R (GGA→AGA)	insH1 ←	IS5 transposase and trans‑activator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,067,300	0	C	T	51.6%	‑4.5 / 86.3	64	G8R (GGA→AGA)	insH1	IS5 transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base C (17/14);  new base T (20/13);  total (37/27)
Fisher's exact test for biased strand distribution p-value = 8.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

ATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGTCACAGC  >  NC_000913/2067203‑2067398                                                                                                  |                                                                                                   aTGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCAt                                                                                                  <  2:1382064/100‑1 (MQ=32)    cGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGc                                                                                               <  2:1665687/100‑1 (MQ=32)        aCAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCa                                                                                           <  1:679992/100‑1 (MQ=33)         cAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAg                                                                                          >  2:2637596/1‑100 (MQ=33)          aaGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGTCGAAGGTAAGTTGATGACTCATTATGAACCCAGTTCCATCGCTCCAGa                                                                                         >  2:1067988/1‑100 (MQ=11)           aGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGAt                                                                                        >  2:1166437/1‑100 (MQ=33)                  ctctTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAAc                                                                                 <  2:31408/100‑1 (MQ=33)                   tctTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTCTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTACAGATGACAAACa                                                                                >  1:2893312/1‑100 (MQ=16)                        tCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGAtc                                                                           <  2:2815857/100‑1 (MQ=33)                                     cgcTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGAc                                                              <  2:1035759/100‑1 (MQ=33)                                     cgcTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGAc                                                              <  1:2672452/100‑1 (MQ=33)                                       cTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACtt                                                            >  2:467859/1‑100 (MQ=32)                                         tACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGt                                                          >  2:1762042/1‑100 (MQ=33)                                           ctgctgAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAACCATGATCTCATATCAGGGACTTGTTc                                                        >  1:1000318/1‑100 (MQ=18)                                           ctgctgAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTc                                                        <  2:1393848/100‑1 (MQ=33)                                             gctgAATTCACTGTCGGCGAAGGTAAGTTGATGCCTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGc                                                      >  1:1992589/1‑100 (MQ=18)                                             gctgAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGc                                                      >  1:1992600/1‑100 (MQ=32)                                              ctgAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCa                                                     <  2:1365162/100‑1 (MQ=33)                                              ctgAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCa                                                     <  2:2801248/100‑1 (MQ=33)                                                gAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCAcc                                                   >  2:952431/1‑100 (MQ=33)                                                gAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCAcc                                                   >  2:1918084/1‑100 (MQ=32)                                                         gTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGa                                          >  1:2871365/1‑100 (MQ=11)                                                           cGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAac                                        <  2:677978/100‑1 (MQ=14)                                                           cGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAac                                        <  1:1774473/100‑1 (MQ=255)                                                            ggCGAAGGTAAGTTGATGGCTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAac                                        <  2:2259484/99‑1 (MQ=255)                                                            ggCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAaca                                       >  2:2462381/1‑100 (MQ=18)                                                            ggCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAaca                                       >  2:588412/1‑100 (MQ=18)                                                            ggCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAaca                                       <  2:92063/100‑1 (MQ=40)                                                              cGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAAcacc                                     >  2:2541504/1‑100 (MQ=255)                                                               gAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAAcacca                                    <  2:2130918/100‑1 (MQ=38)                                                                  ggTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACaa                                 >  1:2705111/1‑100 (MQ=255)                                                                   gTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAAt                                >  1:590199/1‑100 (MQ=255)                                                                   gTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAAt                                >  2:2379534/1‑100 (MQ=255)                                                                   gTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTACCTAGAACACCACAAt                                >  1:1702284/1‑100 (MQ=255)                                                                       gTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCg                            <  2:1341819/100‑1 (MQ=255)                                                                        ttGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATTAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGc                           <  2:92003/100‑1 (MQ=255)                                                                        ttGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGc                           <  1:294964/100‑1 (MQ=255)                                                                        ttGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGc                           <  1:295081/100‑1 (MQ=255)                                                                        ttGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGc                           <  2:862602/100‑1 (MQ=255)                                                                        ttGATGACGCATGATGAACTCTGTTCTATGGCTCCAGAGGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGATCACCACCATTTCGc                           <  1:684968/100‑1 (MQ=14)                                                                         tgatgaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGct                          <  1:1251008/100‑1 (MQ=255)                                                                          gatgaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGctc                         <  2:2738499/100‑1 (MQ=255)                                                                          gatgaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGc                           <  1:433267/98‑1 (MQ=255)                                                                             gaCTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGctctct                      >  2:341410/1‑100 (MQ=255)                                                                               cTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCg                    >  1:502312/1‑100 (MQ=255)                                                                               cTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCg                    >  1:551293/1‑100 (MQ=255)                                                                               cTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCg                    >  1:901737/1‑100 (MQ=255)                                                                               cTCATCATGAACCCTGTTCTATGGCTCCAGATGACAAACATGCTCTCATATGAGCGACTCGTTCGCACCTTCCCTACAACACCGCTTTCTCGCtctctc                     >  1:1509076/1‑99 (MQ=14)                                                                                 cATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGt                  >  2:594665/1‑100 (MQ=255)                                                                                   tgatgaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTaa                >  2:1597875/1‑100 (MQ=255)                                                                                   tgatgaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTaa                >  1:2445391/1‑100 (MQ=255)                                                                                   tgatgaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTa                 >  1:2523069/1‑99 (MQ=255)                                                                                     attaaCCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAgg              <  2:912424/97‑1 (MQ=255)                                                                                     atgaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAg               >  2:654356/1‑99 (MQ=255)                                                                                       gaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCa            >  1:2574279/1‑100 (MQ=255)                                                                                        aaCCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCAc           <  2:805906/100‑1 (MQ=255)                                                                                         aCCCTGTTCTATGGCTCCAGATGACCAACATGATATCATATAAGGGACTTGTTCGCACCTTCCCTAGAACACAAAAATTTCGCTCTCTCGGTAAGGCACt          >  1:1423184/1‑100 (MQ=18)                                                                                         aCCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAATGCACt          >  1:279960/1‑100 (MQ=255)                                                                                          cccTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTg         >  1:262676/1‑100 (MQ=255)                                                                                          cccTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTg         >  1:2623300/1‑100 (MQ=255)                                                                                           ccTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGt        >  2:2879900/1‑100 (MQ=255)                                                                                            cTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCACTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGTc       >  2:1386611/1‑100 (MQ=255)                                                                                              gTTCCATTGCTCCATATGACTAACATCATCTCAAATCAGGGACTTGTCAACACTTACCCTAGAACACCACCATTTCGCTCTCTCGGTAAGGCACTGTcac     <  1:353612/100‑1 (MQ=255)                                                                                                tCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGTCACAg   >  2:248325/1‑100 (MQ=255)                                                                                                 ccATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGTCACAGc  >  2:1898753/1‑100 (MQ=255)                                                                                                 ccATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGTCACAGc  >  2:1899282/1‑100 (MQ=255)                                                                                                  |                                                                                                   ATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGTCACAGC  >  NC_000913/2067203‑2067398

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