Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,187,415 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_0009133,187,4150GA98.2% 160.2 / ‑4.5 54*302* (TAG→TAA) insD1IS2 transposase TnpB
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (28/25);  total (29/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

AGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGG  >  NC_000913/3187292‑3187548
                                                                                                                           |                                                                                                                                     
agagGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTa                                                                                                                        <  2:799842/139‑1 (MQ=17)
       ttCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATcatcat                                                                                                                 >  1:660215/1‑139 (MQ=38)
                 ataACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGGATATCTGCGGCAGCGTGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCg                                                                                                       <  1:273711/139‑1 (MQ=255)
                 ataACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCAACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAGCAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCg                                                                                                       <  1:134351/139‑1 (MQ=255)
                  taACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGt                                                                                                      <  1:1232130/139‑1 (MQ=255)
                  taACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGt                                                                                                      >  2:1391529/1‑139 (MQ=255)
                          ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTatc                                                                                                                      <  2:1294707/115‑1 (MQ=21)
                          ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTatc                                                                                                                      >  1:1294707/1‑115 (MQ=21)
                          ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTAt                                                                                              >  1:770706/1‑139 (MQ=255)
                          ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTAt                                                                                              >  1:1677638/1‑139 (MQ=255)
                           gCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGACAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATg                                                                                             <  2:770706/139‑1 (MQ=255)
                                         cgcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTa                                                                                                                        >  1:447733/1‑98 (MQ=18)
                                         cgcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTa                                                                                                                        <  2:447733/98‑1 (MQ=18)
                                         cgcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTcaac                                                                               <  1:389372/139‑1 (MQ=255)
                                         cgcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTcaac                                                                               >  1:1557554/1‑139 (MQ=255)
                                          gcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTcaaca                                                                              >  2:419718/1‑139 (MQ=255)
                                                      ctcgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCg                                                                  >  2:1281017/1‑139 (MQ=255)
                                                          ccACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTcaaca                                                                              <  2:489804/123‑1 (MQ=255)
                                                          ccACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTcaaca                                                                              >  1:489804/1‑123 (MQ=255)
                                                                     tCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAAt                                                                    >  2:675144/1‑122 (MQ=255)
                                                                     tCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAAt                                                                    <  1:675144/122‑1 (MQ=255)
                                                                     tCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTg                                                   <  1:419718/139‑1 (MQ=255)
                                                                               cGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTatcat                                                                                                                    >  2:584373/1‑64 (MQ=255)
                                                                               cGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTatcat                                                                                                                    <  1:584373/64‑1 (MQ=255)
                                                                               cGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTa                                                                                               <  2:1585554/85‑1 (MQ=255)
                                                                               cGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTa                                                                                               >  1:1585554/1‑85 (MQ=255)
                                                                               cGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTaa                                         <  2:204186/139‑1 (MQ=255)
                                                                               cGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTaa                                         >  2:1665554/1‑139 (MQ=255)
                                                                               cGGGCTTGTAATGGGTTAAGTGATAACAGAAGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCTGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTaa                                         >  1:1708882/1‑139 (MQ=255)
                                                                                 ggCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAAcc                                       >  2:23922/1‑139 (MQ=255)
                                                                                      gTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTcc                                  <  1:56991/139‑1 (MQ=255)
                                                                                      gTAATGGGTCAAGTGATAACAGATGTCTGTAAATAGTAGGGCAAATCCAGTGATCATTATTGGTTGTGCGTCAACTTATGCCGTGGAATTCAACAAAGATTAAATCCAAGCCGAAGATCGGGAAATCGTTAACCTTTcc                                  <  1:1281017/139‑1 (MQ=37)
                                                                                            ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAAttt                            >  1:892857/1‑139 (MQ=255)
                                                                                            ggTTAAGTGATAACAGATGTCTGGAAATACAAGGGCAAATCCAGTGATCATCTTTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAAttt                            >  1:1086991/1‑139 (MQ=255)
                                                                                                aaGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGtt                                           <  2:603934/120‑1 (MQ=255)
                                                                                                aaGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGtt                                           >  1:603934/1‑120 (MQ=255)
                                                                                                 aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTgg                                                                                                              <  1:503321/52‑1 (MQ=255)
                                                                                                 aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTgg                                                                                                              >  2:503321/1‑52 (MQ=255)
                                                                                                 aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAAc                       >  1:225765/1‑139 (MQ=255)
                                                                                                 aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTAATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAAc                       >  1:1281677/1‑139 (MQ=255)
                                                                                                     aTAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGAt                   <  2:1708882/139‑1 (MQ=255)
                                                                                                      tAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATg                  <  2:1557554/139‑1 (MQ=255)
                                                                                                      tAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATg                  >  2:436565/1‑139 (MQ=255)
                                                                                                        aCAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGc                >  1:708790/1‑139 (MQ=255)
                                                                                                            aTGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAAt            >  2:1504754/1‑139 (MQ=255)
                                                                                                            aTGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAAt            >  2:821966/1‑139 (MQ=255)
                                                                                                            aTGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAAt            >  2:221607/1‑139 (MQ=255)
                                                                                                             tGTCTGGAAATATAGGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAAtt           >  1:1066424/1‑139 (MQ=255)
                                                                                                              gTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAAGAAAGATTTAATCg                                                                  >  1:1095355/1‑83 (MQ=255)
                                                                                                              gTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCg                                                                  <  2:1095355/83‑1 (MQ=255)
                                                                                                                  ggAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCg      <  2:460772/139‑1 (MQ=255)
                                                                                                                  ggAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCg      <  2:708790/139‑1 (MQ=255)
                                                                                                                      atatAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCgg  <  1:436565/139‑1 (MQ=255)
                                                                                                                      atatAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCgg  <  2:1677638/139‑1 (MQ=255)
                                                                                                                           |                                                                                                                                     
AGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGG  >  NC_000913/3187292‑3187548

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