Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	382,579	G→A	*302* (TAG→TAA)	insD1 →	IS2 transposase TnpB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	382,579	0	G	A	97.7%	135.2 / ‑6.8	43	*302* (TAG→TAA)	insD1	IS2 transposase TnpB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (22/20);  minor base C (0/1);  total (22/21)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTAC  >  NC_000913/382463‑382704                                                                                                                     |                                                                                                                              ttCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCa                                                                                                         >  2:557125/1‑139 (MQ=37)            taACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTgcg                                                                                              <  2:346914/139‑1 (MQ=255)                    ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCAc                                                                                      >  1:306559/1‑139 (MQ=255)                    ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCAc                                                                                      >  2:110395/1‑139 (MQ=255)                    ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCAc                                                                                      >  1:987744/1‑139 (MQ=255)                    ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCAc                                                                                      >  2:831520/1‑139 (MQ=255)                      cATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTa                                                                                    >  1:802281/1‑139 (MQ=255)                       aTCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTgcg                                                                                              <  2:309855/128‑1 (MQ=255)                       aTCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTgcg                                                                                              >  1:309855/1‑128 (MQ=255)                        tCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTAcc                                                                                  <  2:324733/139‑1 (MQ=255)                        tCCGCATAGGGCGGCGGGTTATCGCTAGCCAAGGGAATATCTACGGCAACGGGGTCGTCATGGGTTCAGTGATAACAGATGTCTGGACATATACAAGCAAATATAATTATCTATCAGGCCGTTTGCGCCAATCACTAcc                                                                                  <  2:411581/139‑1 (MQ=12)                          cGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAgg                                                                                                       <  1:375659/116‑1 (MQ=255)                          cGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAgg                                                                                                       >  2:375659/1‑116 (MQ=255)                                gTGCGTTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAAt                                                                                         <  2:368226/124‑1 (MQ=255)                                gTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAAt                                                                                         >  1:368226/1‑124 (MQ=255)                                 tGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGAAGCGGACTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATAACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGa                                                                         >  1:705840/1‑139 (MQ=255)                                          ttATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCa                                                                <  2:403304/139‑1 (MQ=255)                                            aTCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCAcc                                                              <  2:987744/139‑1 (MQ=255)                                             tcgctcgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCa                                                             >  1:444750/1‑139 (MQ=255)                                               gctcgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATg                                                           <  1:831520/139‑1 (MQ=255)                                                 tcgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGAc                                                         <  2:27871/139‑1 (MQ=255)                                                   gcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACaa                                                       >  2:535377/1‑139 (MQ=255)                                                    ccACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACaaa                                                      <  2:306559/139‑1 (MQ=255)                                                             tatCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTAt                                             >  1:885695/1‑139 (MQ=255)                                                                                      ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGc                    >  1:999788/1‑139 (MQ=255)                                                                                      ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGc                    >  2:332657/1‑139 (MQ=255)                                                                                         tAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGt                 <  2:999788/139‑1 (MQ=255)                                                                                          aaGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTg                                                                          <  2:294866/81‑1 (MQ=255)                                                                                          aaGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTg                                                                          >  1:294866/1‑81 (MQ=255)                                                                                          aaGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCAcc                                                              >  2:270592/1‑93 (MQ=255)                                                                                          aaGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCAcc                                                              <  1:270592/93‑1 (MQ=255)                                                                                          aaGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGtt                <  2:766843/139‑1 (MQ=255)                                                                                           aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGCTTGCTCCAATCACTACCGGACCCTGAATATCGGCACCATGACAAATACGGGTATTCGCGCTGATAACAACTTCACCTGCAGttt               >  2:695527/1‑139 (MQ=39)                                                                                              gATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTtcat            <  2:2664/139‑1 (MQ=255)                                                                                              gATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTtcat            <  2:444750/139‑1 (MQ=255)                                                                                                 aaCAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCTGACCCTGAATAACGGCACCATGACAAATACTGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCaa         <  1:173957/139‑1 (MQ=255)                                                                                                  aCAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCaat        >  2:1050929/1‑139 (MQ=255)                                                                                                     gATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCaataat     <  1:270420/139‑1 (MQ=255)                                                                                                       tGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTa   >  2:147171/1‑139 (MQ=255)                                                                                                       tGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTa   >  1:702923/1‑139 (MQ=255)                                                                                                        gTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTAc  <  2:137571/139‑1 (MQ=255)                                                                                                           tGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGcc                                      <  2:720979/100‑1 (MQ=255)                                                                                                           tGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGcc                                      >  1:720979/1‑100 (MQ=255)                                                                                                                     |                                                                                                                              TTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTAC  >  NC_000913/382463‑382704

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑