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,4150GA100.0% 144.6 / NA 46*302* (TAG→TAA) insD1IS2 transposase TnpB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (25/21);  total (25/21)

AGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGC  >  NC_000913/3187292‑3187549
                                                                                                                           |                                                                                                                                      
agagGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTa                                                                                                                         <  1:108436/139‑1 (MQ=17)
       ttCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGACTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATcatcat                                                                                                                  >  1:479896/1‑139 (MQ=255)
                          ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTatc                                                                                                                       <  1:1234675/115‑1 (MQ=21)
                          ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTatc                                                                                                                       >  2:1234675/1‑115 (MQ=21)
                          ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTAt                                                                                               >  2:430264/1‑139 (MQ=255)
                                       tGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGtt                                                                                                             <  2:764306/112‑1 (MQ=255)
                                       tGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGtt                                                                                                             >  1:764306/1‑112 (MQ=255)
                                         cgcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTcaac                                                                                >  2:260259/1‑139 (MQ=255)
                                         cgcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTcaac                                                                                >  1:142322/1‑139 (MQ=255)
                                          gcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTCTCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTcaaca                                                                               <  1:878618/139‑1 (MQ=255)
                                           cTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTgg                                                                                                               >  2:943650/1‑106 (MQ=255)
                                           cTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTgg                                                                                                               <  1:943650/106‑1 (MQ=255)
                                                   tcgctcgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTaa                                                                      >  1:582101/1‑139 (MQ=255)
                                                        cgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGaa                                                                 >  2:283167/1‑139 (MQ=255)
                                                            aCGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCg                                                             <  2:582101/139‑1 (MQ=255)
                                                                gAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGa                                                         <  1:283167/139‑1 (MQ=255)
                                                                     tCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTg                                                    >  2:864132/1‑139 (MQ=255)
                                                                        gCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCg                                                                   <  1:765230/121‑1 (MQ=255)
                                                                        gCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCg                                                                   >  2:765230/1‑121 (MQ=255)
                                                                         cGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGtt                                                                                                             <  2:952342/78‑1 (MQ=255)
                                                                         cGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGtt                                                                                                             >  1:952342/1‑78 (MQ=255)
                                                                               cGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCAtt                                                                                                                 <  2:1041794/68‑1 (MQ=255)
                                                                               cGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCAtt                                                                                                                 >  1:1041794/1‑68 (MQ=255)
                                                                               cGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAAt                                                                     >  2:1173979/1‑112 (MQ=255)
                                                                               cGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAAt                                                                     <  1:1173979/112‑1 (MQ=255)
                                                                               cGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCg                                                             >  2:1315397/1‑120 (MQ=255)
                                                                               cGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCg                                                             <  1:1315397/120‑1 (MQ=255)
                                                                                      gTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTcc                                   <  2:479896/139‑1 (MQ=255)
                                                                                        aaTGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCa                                 <  1:260259/139‑1 (MQ=255)
                                                                                            ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCg                                                                                                        <  1:359028/64‑1 (MQ=255)
                                                                                            ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCg                                                                                                        >  2:359028/1‑64 (MQ=255)
                                                                                            ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAAttt                             >  1:97170/1‑139 (MQ=255)
                                                                                                aaGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTatc                                                                                                                       >  1:4596/1‑45 (MQ=255)
                                                                                                aaGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTatc                                                                                                                       <  2:4596/45‑1 (MQ=255)
                                                                                                 aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAAc                        >  1:1073092/1‑139 (MQ=255)
                                                                                                       aaCAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATgg                  >  1:930422/1‑139 (MQ=255)
                                                                                                           gATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGc                                                                                                   <  2:659494/54‑1 (MQ=255)
                                                                                                           gATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGc                                                                                                   >  1:659494/1‑54 (MQ=255)
                                                                                                            aTGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCg                                                                   >  2:1296325/1‑85 (MQ=255)
                                                                                                            aTGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCg                                                                   <  1:1296325/85‑1 (MQ=255)
                                                                                                            aTGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAAt             >  2:676453/1‑139 (MQ=255)
                                                                                                             tGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAAtt            >  1:44400/1‑139 (MQ=255)
                                                                                                                 tGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGa                                                   >  1:921869/1‑96 (MQ=255)
                                                                                                                 tGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGa                                                   <  2:921869/96‑1 (MQ=255)
                                                                                                                  ggAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCg       <  2:44400/139‑1 (MQ=255)
                                                                                                                       tataAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGc  <  2:97170/139‑1 (MQ=255)
                                                                                                                           |                                                                                                                                      
AGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGC  >  NC_000913/3187292‑3187549

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