Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,499,871 G→A intergenic (+371/‑1389) insD6 → / → yjgZ IS2 insertion element protein InsB/uncharacterized protein YjgZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,499,8710GA100.0% 26.3 / NA 12intergenic (+371/‑1389)insD6/yjgZIS2 insertion element protein InsB/uncharacterized protein YjgZ
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (5/7);  total (5/7)

CCGGTATATTCCTCTCTCCCTTAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGCTAGGAGTTAAGGCTGTCACACGGATTTGGATG  >  NC_000913/4499640‑4500075
                                                                                                                                                                                                                                       |                                                                                                                                                                                                            
ccGGTATATTCCTCTCTCCCTTAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGAATTGGCACTGTCGCTGAcc                                                                                                                                                                                           <  1:99511/251‑1 (MQ=255)
     atatTCCTCTCTCCCTTAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGAATTgg                                                                                                                                                                                                         <  1:34005/232‑1 (MQ=255)
                                                                                              ggTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACAAGACCACAATGATGCACAGAATTGGCACTGTCGCTGACCGACAGCCCAAGGGCCGAAAAGTAAAGCATGGCGAAATGACGGGAGGACCGGGGAAGGATACCAACAAGCTGGCAAAGGAAGATATCCAggtttc                                                                                             >  1:34680/1‑246 (MQ=255)
                                                                                               gtaatccacgtaaaaagcccTGAGAAATTACGGGCACTGATTTTTCTACTGACGGGCTTATGAGTCGTTATATCACGGGCGGTAAGAATCAGTGTCACACTGTAGACTCCCCGCCCAGGACCATAATGATGAACAGAATTGGCACTGTCGCTGACCGACCGCCCAATGCCAGAAAAGTAATGCATGGCTAAATGACGTGAGTTACGTCGAATGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACt                                                                                            <  1:172655/233‑1 (MQ=255)
                                                                                                                      gAAATGAAGGGCACTGATTTTTCTACTCACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAATGGCACCCTGTAGTCTCCCTGCACATGGCCATCATGGTGAACAAAATTGGCACTGTCGCCGACCGACAGCCCCAGACCAGGAAAGTAAAGCATGGCGAAATGACGGGAGTACCGGCGAAGGAAACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGAttccgacggagtaacatgca                                                                      >  1:86263/1‑231 (MQ=255)
                                                                                                                                                                                         aGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGAATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAAGTTTTCGGGGATACTGCTCCAGGTAATTATTCGGCTAGGAGTTAAGGCTGTCACACGGATTTGGATg  <  1:87295/251‑1 (MQ=255)
                                                                                                                                                                                         aGTGGCACACTGTAGCCTCCCTGCACAGGACCATAATGATGAACAGAATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGCTAGGAGTTAATGCTGTCACACGGATTTGGATg  <  1:69868/251‑1 (MQ=255)
                                                                                                                                                                                              cacaCTGTAGCCTCCCTGCACATGACCATAATGATGAACAGAATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGCTAGGAGTTAAGGCTGTc                 >  1:132792/1‑231 (MQ=255)
                                                                                                                                                                                                    gTAGCCTCCCTGCACATGACCATAATGATGAACAGAATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGANAGtnnnnn                                                                                                    >  2:153796/1‑137 (MQ=255)
                                                                                                                                                                                                    gTAGCCTCCCTGCACATGACCATAATGATGAACAGAATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAggn                                                                       <  1:153796/171‑2 (MQ=255)
                                                                                                                                                                                                                  caTGACCATAATGATGAACAGAATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATCCCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAgn                                                                        <  1:149816/156‑2 (MQ=255)
                                                                                                                                                                                                                  caTGACCATAATGATGAACAGAATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAACGTTATCCAGTTTNCTGNNNNNggnn                                                                                  >  2:149816/1‑144 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                            
CCGGTATATTCCTCTCTCCCTTAAAAACCATAACAGGTTAGTTATCTTCGTCTGTCAGCGGATTGAATGGCGGTATATTTTCAGTACAAATACCGGTCATCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGCTAGGAGTTAAGGCTGTCACACGGATTTGGATG  >  NC_000913/4499640‑4500075

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