Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 315,244 G→A intergenic (+16/‑47) eaeH → / → insE1 pseudogene, attaching and effacing protein homology,factor, Not classified/IS3 transposase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913315,2440GA100.0% 131.7 / NA 43intergenic (+16/‑47)eaeH/insE1pseudogene, attaching and effacing protein homology,factor, Not classified/IS3 transposase A
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/19);  total (24/19)

AGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGC  >  NC_000913/315108‑315366
                                                                                                                                        |                                                                                                                          
aGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATg                                                                                                                          <  2:954291/139‑1 (MQ=255)
       cgcAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACAcgc                                                                                                                   >  2:992060/1‑139 (MQ=255)
                 aTGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGa                                                                                                         >  1:280024/1‑139 (MQ=255)
                      tGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTaaa                                                                                                    <  2:298639/139‑1 (MQ=255)
                      tGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTaaa                                                                                                    >  1:298639/1‑139 (MQ=255)
                           aGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACtctc                                                                                               >  1:864271/1‑139 (MQ=255)
                             tattatGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCgg                                                                                                                 >  2:1020545/1‑119 (MQ=255)
                             tattatGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCgg                                                                                                                 <  1:1020545/119‑1 (MQ=255)
                             tattatGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCtcag                                                                                             >  1:456261/1‑139 (MQ=255)
                             tattatGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCtcag                                                                                             >  1:389751/1‑139 (MQ=255)
                              attatGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCtcagt                                                                                            >  2:285904/1‑139 (MQ=255)
                                  tGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCaga                                                                                        >  1:971615/1‑139 (MQ=255)
                                  tGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCaga                                                                                        >  1:215494/1‑139 (MQ=255)
                                       aTGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGa                                                                                   >  2:274469/1‑139 (MQ=255)
                                            gTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTcac                                                                              <  1:992060/139‑1 (MQ=255)
                                              cGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACAt                                                                            <  1:274469/139‑1 (MQ=255)
                                                      ttGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACaaaaa                                                                    <  2:555560/139‑1 (MQ=255)
                                                          tAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGt                                                                >  2:292881/1‑139 (MQ=255)
                                                               aTAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCaa                                                           >  1:602070/1‑139 (MQ=255)
                                                                tAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAAc                                                          <  1:285904/139‑1 (MQ=255)
                                                                    cgcCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGt                                                      >  1:86075/1‑139 (MQ=255)
                                                                         gAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGacac                                                                                                                     >  1:822581/1‑71 (MQ=255)
                                                                         gAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGacac                                                                                                                     <  2:822581/71‑1 (MQ=255)
                                                                         gAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTaaaaa                                                 <  2:215494/139‑1 (MQ=255)
                                                                              aCCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTcac                                                                              <  2:362870/105‑1 (MQ=255)
                                                                              aCCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTcac                                                                              >  1:362870/1‑105 (MQ=255)
                                                                              aCCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAAcccc                                            <  2:86075/139‑1 (MQ=255)
                                                                                   cATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTaaa                                       <  2:971615/139‑1 (MQ=255)
                                                                                   cATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTaaa                                       >  2:187398/1‑139 (MQ=255)
                                                                                        tATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCa                                  <  1:187398/139‑1 (MQ=255)
                                                                                                  gAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGaa                        >  1:568006/1‑139 (MQ=255)
                                                                                                    ggTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGac                                                                                                                       >  1:830621/1‑42 (MQ=255)
                                                                                                    ggTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGac                                                                                                                       <  2:830621/42‑1 (MQ=255)
                                                                                                            ataCGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTg              <  2:568006/139‑1 (MQ=255)
                                                                                                             taCGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGa             >  1:109852/1‑139 (MQ=255)
                                                                                                             taCGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGa             >  1:316267/1‑139 (MQ=255)
                                                                                                                  cAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCg                   >  2:141414/1‑128 (MQ=25)
                                                                                                                  cAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCg                   <  1:141414/128‑1 (MQ=25)
                                                                                                                   aGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTaaa                                                                                                    >  2:709908/1‑46 (MQ=255)
                                                                                                                   aGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTaaa                                                                                                    <  1:709908/46‑1 (MQ=255)
                                                                                                                     tGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGc                             >  1:161131/1‑115 (MQ=18)
                                                                                                                     tGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGc                             <  2:161131/115‑1 (MQ=18)
                                                                                                                        ctctTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGc  <  2:389751/139‑1 (MQ=11)
                                                                                                                                        |                                                                                                                          
AGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGC  >  NC_000913/315108‑315366

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: