Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,499,500 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_0009134,499,5000GA100.0% 127.4 / NA 41*302* (TAG→TAA) insD1IS2 transposase TnpB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/20);  total (21/20)

GCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACT  >  NC_000913/4499375‑4499634
                                                                                                                             |                                                                                                                                      
gCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAcc                                                                                                                           <  2:268854/139‑1 (MQ=11)
                 ttATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATa                                                                                                          >  2:183682/1‑139 (MQ=255)
                              cATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTtata                                                                                             <  1:1509579/139‑1 (MQ=255)
                                        gTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACtt                                                                                   >  1:1511523/1‑139 (MQ=255)
                                             cTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                                       >  1:955917/1‑114 (MQ=255)
                                             cTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                                       <  2:955917/114‑1 (MQ=255)
                                                   tATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAcc                                                                                                                           <  2:1418870/88‑1 (MQ=14)
                                                   tATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAcc                                                                                                                           >  1:1418870/1‑88 (MQ=14)
                                                   tATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAg                                                                        <  1:183682/139‑1 (MQ=255)
                                                          cgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCGACCTtgtg                                                                                                                      >  1:1042744/1‑86 (MQ=255)
                                                          cgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTCATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                                      <  2:1042744/86‑1 (MQ=255)
                                                                gggAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACtttt                                                                                 <  2:836135/117‑1 (MQ=255)
                                                                gggAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACtttt                                                                                 >  1:836135/1‑117 (MQ=255)
                                                                     tatCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTg                                                      <  1:1645922/139‑1 (MQ=255)
                                                                       tCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGcc                                                    >  2:1116879/1‑139 (MQ=255)
                                                                            ggCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCt                                                                                                                    >  1:1423736/1‑70 (MQ=255)
                                                                            ggCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCt                                                                                                                    <  2:1423736/70‑1 (MQ=255)
                                                                            ggCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGaa                                                                                                            >  1:1244262/1‑78 (MQ=255)
                                                                            ggCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGaa                                                                                                            <  2:1244262/78‑1 (MQ=255)
                                                                            ggCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGa                                                                     >  1:1236609/1‑117 (MQ=255)
                                                                            ggCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGa                                                                     <  2:1236609/117‑1 (MQ=255)
                                                                                 cGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGAt                                           >  2:610800/1‑138 (MQ=255)
                                                                                 cGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGAt                                           <  1:610800/138‑1 (MQ=255)
                                                                                     cTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTtata                                                                                             >  2:8780/1‑84 (MQ=255)
                                                                                     cTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTtata                                                                                             <  1:8780/84‑1 (MQ=255)
                                                                                           aTGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGc                                                                                                                     <  1:1090118/54‑1 (MQ=255)
                                                                                           aTGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGc                                                                                                                     >  2:1090118/1‑54 (MQ=255)
                                                                                              ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCt                                                                                                                    >  2:1348480/1‑52 (MQ=255)
                                                                                              ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCt                                                                                                                    <  1:1348480/52‑1 (MQ=255)
                                                                                                         aaCAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGa                                                                     <  2:154514/88‑1 (MQ=255)
                                                                                                         aaCAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGa                                                                     >  1:154514/1‑88 (MQ=255)
                                                                                                              aTGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGa                                                                     <  1:823976/83‑1 (MQ=255)
                                                                                                              aTGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGa                                                                     >  2:823976/1‑83 (MQ=255)
                                                                                                                gTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCtgtg                                   >  1:284709/1‑115 (MQ=255)
                                                                                                                gTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCtgtg                                   <  2:284709/115‑1 (MQ=255)
                                                                                                                 tCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGa                                                                     >  1:88317/1‑80 (MQ=255)
                                                                                                                 tCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGa                                                                     <  2:88317/80‑1 (MQ=255)
                                                                                                                   tGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATa                                                                                      <  2:1332910/61‑1 (MQ=255)
                                                                                                                   tGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATa                                                                                      >  1:1332910/1‑61 (MQ=255)
                                                                                                                      aaaTATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCaa     >  1:418858/1‑139 (MQ=255)
                                                                                                                         tataAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACt  >  1:1177399/1‑139 (MQ=255)
                                                                                                                             |                                                                                                                                      
GCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACT  >  NC_000913/4499375‑4499634

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