Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,008,400 +GAAA coding (483/516 nt) b4555 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,008,3981.A100.0% 132.3 / NA 40E161K (GAA→AAA) b4555hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (18/22);  total (18/22)
*NC_0009134,008,3982.A100.0% 134.6 / NA 40E161K (GAA→AAA) b4555hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (18/22);  total (18/22)
*NC_0009134,008,3983.G100.0% 124.3 / NA 40E161? (GAA→NAA) b4555hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (18/22);  total (18/22)
*NC_0009134,008,3984.A100.0% 137.6 / NA 40E161K (GAA→AAA) b4555hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (18/22);  total (18/22)

GCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAGGCCATGCGCTGGAAGGCTGGCG‑‑‑‑AACTGGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCT  >  NC_000913/4008309‑4008457
                                                                                            ||||                                                                  
gCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGGGCTTCAAAAAGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑a                                                                   >  2:3326159/1‑93 (MQ=255)
 cGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aa                                                                  <  2:3966728/93‑1 (MQ=255)
 cGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aa                                                                  <  1:1782384/93‑1 (MQ=255)
  gATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aag                                                                 <  1:1177000/93‑2 (MQ=255)
   aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                                                >  1:46973/1‑91 (MQ=255)
   aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                                                >  1:3511183/1‑91 (MQ=255)
   aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAAGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                                                >  1:3215869/1‑91 (MQ=255)
    tCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aagaa                                                               >  2:2281928/1‑90 (MQ=255)
    tCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aagaa                                                               >  2:3445805/1‑90 (MQ=255)
    tCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aagaa                                                               <  1:3329958/93‑4 (MQ=255)
                                     tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtt                                 <  2:1756261/93‑1 (MQ=255)
                                      gtGGAAGTGATGGCGGGCGGCGTGTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGttt                                <  2:1612774/92‑2 (MQ=255)
                                      cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGttt                                <  2:2585349/93‑2 (MQ=255)
                                       tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttg                               >  1:3760002/1‑91 (MQ=255)
                                        ggAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttgg                              >  1:4328546/1‑90 (MQ=255)
                                        ggAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttgg                              >  1:1504553/1‑90 (MQ=255)
                                         gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttggt                             >  2:362797/1‑89 (MQ=255)
                                         gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttggt                             <  2:3225814/93‑5 (MQ=255)
                                               gtGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACggg                      <  1:3290333/92‑1 (MQ=255)
                                               aTGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACggg                      <  1:980410/93‑1 (MQ=255)
                                                tGGCGGGCGGCGTTTCAAAAGGCAATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGa                     <  2:2605564/93‑1 (MQ=255)
                                                 ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAt                    <  2:1250992/93‑1 (MQ=255)
                                                  gcgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATg                   >  2:4734727/1‑93 (MQ=255)
                                                   cgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGa                  <  2:1009336/93‑1 (MQ=255)
                                                   cgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGa                  <  2:2641560/93‑1 (MQ=255)
                                                    gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                 <  1:1004691/93‑1 (MQ=255)
                                                    gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                 <  2:594135/93‑1 (MQ=255)
                                                     ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac                >  2:2025077/1‑93 (MQ=255)
                                                     ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTAGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac                >  2:1242468/1‑93 (MQ=255)
                                                      gcggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacg               <  1:3802206/93‑1 (MQ=255)
                                                       cggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacga              <  1:2414018/93‑1 (MQ=255)
                                                        ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac             <  1:321200/93‑1 (MQ=255)
                                                        ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac             <  2:907022/93‑1 (MQ=255)
                                                         gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgacg            >  1:978131/1‑93 (MQ=255)
                                                         gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgacg            >  1:1848270/1‑93 (MQ=255)
                                                          cgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGc           <  2:1666186/93‑1 (MQ=255)
                                                           gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc          >  1:1095282/1‑93 (MQ=255)
                                                           gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc          <  2:4673438/93‑1 (MQ=255)
                                                              tcaaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaa       >  2:3850580/4‑93 (MQ=255)
                                                              gcaaaaGGCCCTGCGCTGTAAGCGCTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaa       <  2:3861309/90‑1 (MQ=255)
                                                               caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa      >  1:2116878/3‑93 (MQ=255)
                                                               caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa      <  2:643823/91‑1 (MQ=255)
                                                                aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTT‑‑GTTGACGGGAGGAACTACGCCGAAAt     >  2:3805337/2‑93 (MQ=255)
                                                                aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt     >  2:2436174/2‑93 (MQ=255)
                                                                aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt     >  2:208377/2‑93 (MQ=255)
                                                                aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt     <  2:2079373/92‑1 (MQ=255)
                                                                aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt     <  1:3137890/92‑1 (MQ=255)
                                                                 aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg    >  1:3644293/1‑93 (MQ=255)
                                                                 aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg    >  2:991364/1‑93 (MQ=255)
                                                                   aGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt  >  1:1869543/1‑93 (MQ=255)
                                                                                            ||||                                                                  
GCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAGGCCATGCGCTGGAAGGCTGGCG‑‑‑‑AACTGGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCT  >  NC_000913/4008309‑4008457

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