Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,623,585 A→G 100% E63G (GAA→GGA)  yhhH → putative NTF2 fold immunity protein for polymorphic toxin RhsB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,623,5850AG100.0% 113.4 / NA 35E63G (GAA→GGA) yhhHputative NTF2 fold immunity protein for polymorphic toxin RhsB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/17);  total (18/17)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGAACTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATA  >  NC_000913/3623441‑3623730
                                                                                                                                                |                                                                                                                                                 
ggAAATGTGTCATTGTATATTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTa                                                                                                                                               <  2:619675/149‑1 (MQ=255)
 gAACTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAg                                                                                                                                              <  1:617927/149‑1 (MQ=255)
  aaCTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGAGCGCAATGTCATAATAAGATTAAAATAATGCGCAAGTACGAAAGTGAAGGCACATACACCGTAAGAAATTTGGTTAAAAATAAAGCTATAGCACTGGGATTAGc                                                                                                                                             >  2:30362/1‑149 (MQ=255)
    cTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTg                                                                                                                                           <  1:192371/149‑1 (MQ=255)
                  cTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATGTATGTTAAg                                                                                                                             >  2:545793/1‑149 (MQ=255)
                   tttgtatgTTTTTACATACAACTCATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGa                                                                                                                            <  2:140059/149‑1 (MQ=255)
                   gttttttgtttTTACATACAACACATGGGCGCAGTGTAATAACAATATTAACATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGa                                                                                                                            <  2:76033/145‑1 (MQ=255)
                         tgtttTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGtt                                                                                                                      <  2:338199/149‑1 (MQ=255)
                           ttattACATAAAACAAATTGGCGCAGGGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATACGTGAAGAATCTTTAt                                                                                                                    <  2:482387/146‑1 (MQ=255)
                                    acaacaCATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGAt                                                                                                           >  1:33888/1‑149 (MQ=255)
                                     caacaCATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGATATGGACAGGATg                                                                                                          >  1:47020/1‑149 (MQ=255)
                                            tGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTGAGAATCGTTATGGACAGGATGCCGCaga                                                                                                   <  1:470023/149‑1 (MQ=255)
                                             tggCTCAATGTAATAATAATATTAAAATAATGCGCAACCATCACACTGACCTTAAATATACAGTTAGAAATTTGCTTAAAAATAAATCTATAGCATTGGGATTAGCTGAGATATAGGTTAAGAATCGTTATGGACAGGATGCCGCagaa                                                                                                  <  2:84537/148‑1 (MQ=255)
                                                       taataataaTATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAAc                                                                                        <  1:692485/149‑1 (MQ=255)
                                                           gataataTTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGATGAGATATATGTTAATAATCGTTCTGGACAGGAAGCAGGCGAAGAAGAGAAAACATa                                                                                    >  2:621907/2‑149 (MQ=255)
                                                               atatTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGaa                                                                                <  1:541305/149‑1 (MQ=255)
                                                                   tAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCa                                                                            <  1:754290/149‑1 (MQ=255)
                                                                     aaaTAATGCGCAAGCATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCTCAGAAGAAGAAAAACCATACGAAATCACt                                                                          <  1:491505/149‑1 (MQ=255)
                                                                      aataatGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATGGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACCAAGTCaagg                                                                         >  2:119547/1‑146 (MQ=255)
                                                                           tGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACGGGATGCCGCAGAAGAAGAAAAACCATACGAAATCCCTGAGATa                                                                    >  1:99825/1‑149 (MQ=255)
                                                                           tGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTa                                                                    >  1:269654/1‑149 (MQ=255)
                                                                           tGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTACCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAGGAAAAACCATACGAAATCACTGAGTTa                                                                    >  1:545156/1‑149 (MQ=255)
                                                                             cggaaGTATGAAAGTGAAGGTATATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTaac                                                                  <  1:424322/146‑1 (MQ=255)
                                                                             cgcAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTaac                                                                  <  1:50153/149‑1 (MQ=255)
                                                                                             aaGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCGGCAGAAGAAGGAAAACCATACGAAATGAGTGTGTTAACAACAGGTTGGgttgtt                                                  >  2:43025/1‑149 (MQ=255)
                                                                                                gTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACAGAGATAACAACAAGATGGGGTGTTaaa                                               >  1:363673/1‑149 (MQ=255)
                                                                                                        aCCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTAccgt                                       >  2:103162/1‑147 (MQ=255)
                                                                                                              aGAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTc                                 <  2:33888/149‑1 (MQ=255)
                                                                                                               gAAATTTGGTTAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCa                                >  2:521300/1‑149 (MQ=255)
                                                                                                                       gTTAAAAATGAAGCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAt                        >  2:777858/1‑149 (MQ=255)
                                                                                                                         tAAAAATAAAGCTATAGCATTGGGATTAGCTGAGATATATGATAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATAAGAAATAACTGAGTTAACAACAAGTTGGGATGTTGAACGTAACATTCACTCCGACCAAAtag                      >  2:320141/1‑147 (MQ=255)
                                                                                                                                   gCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTAGGAATAACTCAGAGCAGGTGGGTGGTGGGGt            >  1:45566/1‑149 (MQ=255)
                                                                                                                                   gCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGt            >  2:308567/1‑149 (MQ=255)
                                                                                                                                   gCTATAGCATTGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGGGGGGt            >  1:740133/1‑149 (MQ=255)
                                                                                                                                             tGGGATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATtata  <  1:545793/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
GGAACTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATA  >  NC_000913/3623441‑3623730

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