Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,733,703 G→A *1539* (TAG→TAA lhr → putative ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,733,7030GA100.0% 141.4 / NA 44*1539* (TAG→TAAlhrputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/20);  total (24/20)

CGTCGCGAACCACGGCTGCGCTTTACGCTAACAGAAGTGAATGATCTACCGGTCCGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAGAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCGGCTCTTCAGACTTGTATTTAGCGGCAGTTTCTTTGATCAGCTGCTGCAGTTCGCCACG  >  NC_000913/1733575‑1733825
                                                                                                                                  |                                                                                                                          
gcgtccgAAACACGGCCTCGGTTTTAGCTAAAGAGAGGGAAAGGTCTACCGGTCCGGCAAAACCCGAATGTTTCCCCGGCTCGGGAGGCGGGGATTTTAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACt                                                                                                                  <  2:238822/134‑1 (MQ=255)
                  cgcTTTACGCTAACAGAAGTGAATGATCTACCGGTCCGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCg                                                                                                  <  2:1606550/137‑1 (MQ=255)
                  cgcTTTACGCTAACAGAAGTGAATGATCTACCGGTCCGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCa                                                                                                >  1:202444/1‑139 (MQ=255)
                  cgcTTTACGCTAACAGAAGTGAATGATCTACCGGTCCGGCAAACGCCGATGCTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCg                                                                                                  >  1:1606550/1‑137 (MQ=255)
                       tACGCTAACAGAAGTGAATGATCTACCGGTCCGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGt                                                                                           <  1:172079/139‑1 (MQ=255)
                               cAGAAGTGAATGATCTACCGGTCCGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCg                                                                                                  <  2:1419123/124‑1 (MQ=255)
                               cAGAAGTGAATGATCTACCGGTCCGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCg                                                                                                  >  1:1419123/1‑124 (MQ=255)
                                        aaTGATCTACCGGTCCGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAg                                                                                        <  1:1094625/126‑1 (MQ=255)
                                        aaTGATCTACCGGTCCGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAg                                                                                        >  2:1094625/1‑126 (MQ=255)
                                         aTGATCTACCGGTCCGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGa                                                                          >  2:424566/1‑139 (MQ=255)
                                             tCTACCGGTCCGGCAAACGCCGATGTTTACGCGGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGattat                                                                      >  2:985853/1‑139 (MQ=255)
                                                aCCGGTCCGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTcc                                                                   >  2:1546831/1‑139 (MQ=255)
                                                   ggTCCGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCcgcg                                                                <  1:424566/139‑1 (MQ=255)
                                                      ccGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCa                                                                                                <  1:296835/104‑1 (MQ=255)
                                                      ccGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCa                                                                                                >  2:296835/1‑104 (MQ=255)
                                                       cGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAaaga                                                                                                                         >  2:803336/1‑78 (MQ=255)
                                                       cGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAaaga                                                                                                                         <  1:803336/78‑1 (MQ=255)
                                                       cGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGt                                                                                                     <  1:96558/98‑1 (MQ=255)
                                                       cGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGt                                                                                                     >  2:96558/1‑98 (MQ=255)
                                                       cGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCa                                                                                 >  2:508769/1‑118 (MQ=255)
                                                       cGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCa                                                                                 <  1:508769/118‑1 (MQ=255)
                                                       cGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCg                                                            >  1:659555/1‑139 (MQ=255)
                                                       cGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCg                                                            >  1:1238551/1‑139 (MQ=255)
                                                         gCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCGGc                                                          >  1:1288681/1‑139 (MQ=255)
                                                                       tttACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCGGCTCTTCAGACTTGt                                             >  2:1181864/1‑139 (MQ=255)
                                                                               gctgcGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCGGCTCTTCAGACTTGTATTTAGCg                                     >  2:186135/1‑139 (MQ=255)
                                                                                ctgcGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCGGCTCTTCAGACTTGTATTTAGCgg                                    >  2:311636/1‑139 (MQ=255)
                                                                                ctgaGTGAGGCGGGATTTTCAAGTTCGACACAAGAGCTGACTTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATGCGTCAGTATTGCATGACGGATTATTCCTCGTCCGGCTCTTCAGTCTTGTATGTAGTgg                                    >  2:1489271/1‑139 (MQ=255)
                                                                                   cGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGGAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCGGCTCTTCAGACTTGTATTTAGCGGCAg                                 <  2:334521/139‑1 (MQ=255)
                                                                                        ggCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCGGCTCTTCAGACTTGTATTTAGCGGCAGTTTCt                            >  2:827472/1‑139 (MQ=255)
                                                                                             gATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCGGCTCTTCAGGCTTGTATTTAGCGGCAGTTTCTTTGAt                       >  1:1153781/1‑139 (MQ=255)
                                                                                                tttCAAGTTCGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCGGCTCTTCAGACTTGTATTTAGCGGCAGTTTCTTTGATCAg                    <  1:1178236/139‑1 (MQ=255)
                                                                                                       ttCGCCACTAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGCATTATTCCGCGTCCGGCTCTTCAGACTTGTATTTAGCGGCAGTTTCTTTGATCAGCtgctgc             <  1:985853/139‑1 (MQ=255)
                                                                                                         cGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCGGCTCTTCAGACTTGTATTTAGCGGCAGttt                              >  1:509011/1‑120 (MQ=255)
                                                                                                         cGCCACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCGGCTCTTCAGACTTGTATTTAGCGGCAGttt                              <  2:509011/120‑1 (MQ=255)
                                                                                                           ccACAAGGGCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCGGCTCTTCAGACTTGTATTTAGCGGCAGTTTCTTTGATCAGCTGCTGCAGtt         <  2:202444/139‑1 (MQ=255)
                                                                                                                  ggCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCGGCTCTTCAGGCTTGTATTTAGCGGCAGTTTCTTTGATCAGCTGCTGCAGTTCGCCACg  <  2:1153781/139‑1 (MQ=255)
                                                                                                                  ggCTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCGGCTCTTCAGACTTGTATTTAGCGGCAGTTTCTTTGATCAGCTGCTGCAGTTCGCCACg  <  2:1238551/139‑1 (MQ=255)
                                                                                                                    cTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCgg                                                           >  2:133643/1‑80 (MQ=255)
                                                                                                                    cTGGATTGGGGATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCgg                                                           <  1:133643/80‑1 (MQ=255)
                                                                                                                            gggATAAAGAAAGGACTGACGGATGCCCGTTCGCa                                                                                                <  1:1453983/35‑1 (MQ=255)
                                                                                                                            gggATAAAGAAAGGACTGACGGATGCCCGTTCGCa                                                                                                >  2:1453983/1‑35 (MQ=255)
                                                                                                                            gggATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAg                                                                                        <  1:1544245/43‑1 (MQ=255)
                                                                                                                            gggATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAg                                                                                        >  2:1544245/1‑43 (MQ=255)
                                                                                                                            gggATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCa                                                                                 >  1:577187/1‑50 (MQ=255)
                                                                                                                            gggATAAAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCa                                                                                 <  2:577187/50‑1 (MQ=255)
                                                                                                                                  |                                                                                                                          
CGTCGCGAACCACGGCTGCGCTTTACGCTAACAGAAGTGAATGATCTACCGGTCCGGCAAACGCCGATGTTTACGCTGCTGCGTGAGGCGGGATTTTCAAGTTCGCCACAAGGGCTGGATTGGGGATAGAGAAAGGACTGACGGATGCCCGTTCGCATCCGTCAGTATTGCAGGACGGATTATTCCGCGTCCGGCTCTTCAGACTTGTATTTAGCGGCAGTTTCTTTGATCAGCTGCTGCAGTTCGCCACG  >  NC_000913/1733575‑1733825

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