Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 1,095,503 Δ10,045 bp IS3‑mediated ycdU[ymdA] ycdU, serX, ghrA, ycdX, ycdY, ycdZ, csgG, csgF, csgE, csgD, csgB, csgA, csgC, [ymdA]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 1094436–1095503 1105547 10045–11112 5 [4] [0] 26 [insF4]–[ymdA] [insF4],insE4,ycdU,serX,ghrA,ycdX,ycdY,ycdZ,csgG,csgF,csgE,csgD,csgB,csgA,csgC,[ymdA]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 315229 =NA (NA)26 (0.650)
+T
24/480 0.7 100% noncoding (1/1255 nt) IS3 repeat region
?NC_000913 1105548 = 0 (0.000)coding (135/312 nt) ymdA uncharacterized protein YmdA

AATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGCCAAATTGTGGAACCGGCATGTGACGTCAGCACCCAGTCATCACCCGTAGAAATGAACTGCCCACAGAATGGTTCTATTCCGGGCAAAACCTACTCCAGCAAAGCGTTAATGAGCGGCAATGTCAAAAACGCGCAAATAGCCTCAGTAAAAGTGCAGTATCTTGATAAACAGAAAAAGCTGGCTGTAATGAACATCGAATATAACTAAGTTTCTGAACAACTCACGTTGCTGAGCAGAAAAATGCGATTTAACCAAAAAGCCTGCT  >  NC_000913/1105302‑1105783
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                           
ttcagaagacgtctgctgcttttgctgtttactgcgccagttgtagagttgtgattcatacaggctgagttcacgggctgcggcagtaacaccgatgcgttcagcaagcttcagggcttcactgcgaaattcaggcgaatgctgtttacggggttttttactggttgatactgtttttgtcatgtgagtcacctctgactgagagtttactcacttagccgcgtgtccactattgctgggtaagatCAccc                                                                                                                                                                                                                                         >  1:149977‑M2/247‑251 (MQ=255)
      agacgtctgctgattttgctgtttactgcgccagttgtagagttgtgattcatacaggctgagttcacgggctgcggcagtaacaccgatgcgttcagcaagcttcagggcttcactgcgaaattcaggcgaatgctgtttacggggttttttactggttgatactgtttttgtcatgtgagtcacctctgactgagagtttactcacttagccgcgtgtccactattgctgggtaagatCACCCAGtcat                                                                                                                                                                                                                                   >  1:315297‑M2/241‑251 (MQ=255)
          gtctgctgattttgctgtttactgcgccagttgtagagttgtgattcatacaggctgagttcacgggctgcggcagtaacaccgatgcgttcagcaagcttcagggcttcactgcgaaattcaggcgaatgctgtttacggggttttttactggttgatactgtttttgtcatgtgagtcacctctgactgagagtttactcacttagccgcgtgtccactattgctgggtaagatCACCCAGTCATCAcc                                                                                                                                                                                                                               >  1:360326‑M2/237‑251 (MQ=255)
                          gtttactgcgccagttgtagagttgtgattcatacaggctgagttcacgggctgcggcagtaacaccgatgcgttcagcaagcttcagggcttcactgcgaaattcaggcgaatgctgtttacggggttttttactggttgatactgtttttgtcatgtgagtcacctctgactgagagtttactcacttagccgcgtgtccactattgctgggtaagatCACCCAGTCATCACCCGTAGAAATGAACTg                                                                                                                                                                                                                <  2:72689‑M2/30‑1 (MQ=255)
                                   gccagttgtagagttgtgattcatacaggctgagttcacgggctgcggcagtaacaccgatgcgttcagcaagcttcagggcttcactgcgaaattcaggcgaatgctgtttacggggttttttactggttgatactgtttgtgtcatgtgagtcacctctgactgagagtttactcacttagccgcgtgtccactattgctgggtaagatCACCCAGTCATCACCCGTAGAAATGACCTGCCCACAGAAt                                                                                                                                                                                                      <  2:329008‑M2/40‑1 (MQ=255)
                                     cagttgtagagttgtgattcatacaggctgagttcacgggctgcggcagtaacaccgatgcgttcagcaagcttcagggcttcactgcgaaattcaggcgaatgctgtttacggggttttttactggttgatactgtttttgtcatgtgagtcacctctgactgagagtttactcacttagccgcgtgtccactattgctgggtaagatCACCCAGTCATCACCCGTAGAAATGAACTGCCCACAGAATgg                                                                                                                                                                                                    <  1:158011‑M2/42‑1 (MQ=255)
                                              agttgtgattcatacaggctgagttcacgggctgcggcagtaacaccgatgcgttcagcacgcgtcagggcttcactgcgaaattcaggcgaatgctgttttcggggttttttactggttgatactgtttttgtcatgtgagtcacctctgactgagagtttactcacttagccgcgtgtccacgattgctgtgtaagatCACCCAGTCAGCACCCGTAGAAATTAACTGCACACAGAGTAGTTCTATTcc                                                                                                                                                                                           <  2:247137‑M2/51‑1 (MQ=255)
                                              agttgtgattcatacaggctgagttcacgggctgcggcagtaacaccgatgcgttcagcaagcttcagggcttcactgcgaaattcaggcgaatgctgtttacggggttttttactggttgatactgtttttgtcatgtgagtcacctctgactgagagtttactcacttagccgcgtgtccactattgctgggtaagatCACCCAGTCATCACCCGTAGAAATGAACTGCCCACAGAATGGTTCTATTcc                                                                                                                                                                                           <  1:379511‑M2/51‑1 (MQ=255)
                                                                          gggctgcggcagtaacaccgatgcgttcagcaagcttcagggcttcactgcgaaattcaggcgaatgctgtttacggggttttttactggttgatactgtttttgtcatgtgagtcacctctgactgagagtttactcacttagccgcgtgtccactattgctgggtaagatCACCCAGTCATCACCCGTAGAAATGAACTGCCCACAGAATGGTTCTATTCCGGGCAAAACCTACTCCAGCAAAGCGTTa                                                                                                                                                               <  1:119952‑M2/79‑1 (MQ=255)
                                                                                 ggcagtaacaccgatgcgttcagcaagcttcagggcttcactgcgaaattcaggcgaatgctgtttacggggttttttactggttgatactgtttttgtcatgtgagtcacctctgactgagagtttactcacttagccgcgtgtccactattgctgggtaagatCACCCAGTCATCACCCGTAGAAATGAACTGCCCACAGAATGGTTCTATTCCGGGCAAAACCTACTCCAGCAAAGCGTTAATGAGCg                                                                                                                                                        <  2:5681‑M2/86‑1 (MQ=255)
                                                                                 ggcagtaacaccgatgcgttcagcaagcttcagggcttcactgcgaaattcaggcgaatgctgtttacggggttttttactggttgatactgtttttgtcatgtgagtcacctctgactgagagtttactcacttagccgcgtgtccactattgctgggtaagatCACCCAGTCATCACCCGTAGAAATGAACTGCCCACAGAATGGTTCTATTCCGGGCAAAACCTACTCCAGCAAAGCGTTAATGAGCg                                                                                                                                                        <  1:24151‑M2/86‑1 (MQ=255)
                                                                                             gatgcgttcagcaagcttcagggcttcactgcgaaattcaggcgaatgctgtttacggggttttttactggttgatactgtttttgtcatgtgagtcacctctgactgagagtttactcacttagccgcgtgtccactattgctgggtaagatCACCCAGTCATCACCCGTAGAAATGAACTGCCCACAGAATGGTTCTATTCCGGGCAAAACCTACTCCAGCAAAGCGTTAATGAGCGGCAATGTCaaa                                                                                                                                             >  1:82981‑M2/154‑250 (MQ=255)
                                                                                                   ttcagcaagcttcagggcttcactgcgaaattcaggcgaatgctgtttacggggttttttactggttgatactgtttttgtcatgtgagtcacctctgactgagagtttactcacttagccgcgtgtccactattgctgggtaagatCACCCAGTCATCACCCGTAGAAATGAACTGCCCACAGAATGGTTCTATTCCGGGCAAAACCTACTCCAGCAAAGCGTTAATGAGCGGCAATGTCAAAAACgcgc                                                                                                                                      >  1:171913‑M2/148‑251 (MQ=255)
                                                                                                                   gcttcactgcgaaattcaggcgaatgctgtttacggggttttttactggttgatactgtttttgtcatgtgagtcacctctgactgagagtttactcacttagccgcgtgtccactattgctgggtaagatCACCCAGTCATCACCCGTAGAAATGAACTGCCCACAGAATGGTTCTATTCCGGGCAAAACCTACTCCAGCAAAGCGTTAATGAGCGGCAATGTCAAAAACGCGCAAATAGCCTCAGTaaa                                                                                                                      <  2:411606‑M2/120‑1 (MQ=255)
                                                                                                                       cactgcgaaattcaggcgaatgctgtttacggggttttttactggatgatactgtttttgtcatgtgagtcacctctgactgagagtttactcacttagccgcgtgtccactattgctgggtaagatCACCCAGTCATCACCCGTAGAAATGAACTGCCCACAGAATGGTTCTATTCCGGGCAAAACCTACTCCAGCAAAGCGTTAATGAGCGGCAATGTCAAAAACGCGCAAATAGCCTCAGTAAAAGTg                                                                                                                  >  1:81843‑M2/128‑251 (MQ=255)
                                                                                                                           gcgaaattcaggcgaatgctgtttacggggttttttactggttgatactgtttttgtcatgtgagtcacctctgactgagagtttactcacttagccgcgtgtccactattgctgggtaagatCACCCAGTCATCACCCGTAGAAATGAACTGCCCACAGAATGGTTCTATTCCGGGCAAAACCTACTCCAGCAAAGCGTTAATGAGCGGCAATGTCAAAAACGCGCAAATAGCCTCAGTAAAAGTGCa                                                                                                                <  2:82981‑M2/126‑1 (MQ=255)
                                                                                                                                   caggcgaatgctgtttacggggttttttactggttgatactgtttttgtcatgtgagtcacctctgactgagagtttactcacttagccgcgtgtccactattgctgggtaagatCACCCAGTCATCACCCGTAGAAATGAACTGCCCACAGAATGGTTCTATTCCGGGCAAAACCTACTCCAGCAAAGCGTTAATGAGCGGCAATGTCAAAAACGCGCAAATAGCCTCAGTAAAAGTGCAGTATCTTGAt                                                                                                      <  1:222060‑M2/136‑1 (MQ=255)
                                                                                                                                         aatgctgtttacggggttttttactggttgatactgtttttgtcatgtgagtcacctctgactgagagtttactcacttagccgcgtgtccactattgctgggtaagatCACCCAGTCATCACCCGTAGAAATGAACTGCCCACAGAATGGTTCTATTCCGGGCAAAACCTACTTCAGCAAAGCGTTAATGAGCGGCAATGTCAAAAACGCGCAAATAGCCTCAGTAAAAGTGCAGTATCTTGATAAACAg                                                                                                <  2:201447‑M2/142‑1 (MQ=255)
                                                                                                                                             ctgtttacggggttttttactggttgatactgtttttgtcatgtgagtcacctctgactgagagtttactcacttagccgcgtgtccactattgctgggtaagatCACCCAGTCATCACCCGTAGAAATGAACTGCCCACAGAATGGTTCTATTCCGGGCAAAACCTACTCCAGCAAAGCGTTAATGAGCGGCAATGTCAAAAACGCGCAAATAGCCTCAGTAAAAGTGCAGTATCTTGATAAACAGaaa                                                                                             >  1:296381‑M2/106‑250 (MQ=255)
                                                                                                                                                      gggttttttactggttgatactgtttttgtcatgtgagtcacctctgactgagagtttactcacttagccgcgtgtccactattgctgggtaagatCACCCAGTCATCACCCGTAGAAATGAACTGCCCACAGAATGGTTCTATTCCGGGCAAAACCTACTCCAGCAAAGCGTTAATGAGCGGCAATGTCAAAAACGCGCAAATAGCCTCAGTAAAAGTGCAGTATCTTGATAAACAGAAAAAgctggctg                                                                                   <  1:233561‑M2/155‑1 (MQ=255)
                                                                                                                                                                                   tcatgtgagtcacctctgactgagagtttactcacttagccgcgtgtccactattgctgggtaagatCACCCAGTCATCACCCGTAGAAATGAACTGCCCACAGAATGGTTCTATTCCGGGCAAAACCTACTCCAGCAAAGCGTTAATGAGCGGCAATGTCAAAAACGCGCAAATAGCCTCAGTAAAAGTGCAGTATCTTGATAAACAGAAAAAGCTGGCTGTAATGAACATCGAATATAACTAAGTTTCt                                                      >  1:236112‑M2/68‑251 (MQ=255)
                                                                                                                                                                                                        gagagtttactcacttagccgcgtgtccactattgctgggtaagatCACCCAGTCATCACCCGTAGAAATGAACTGCCCACAGAATGGTTCTATTCCGGGCAAAACCTACTCCAGCAAAGCGTTAATGAGCGGCAATGTCAAAAACGCGCAAATAGCCTCAGTAAAAGTGCAGTATCTTGATAAACAGAAAAAGCTGGCTGTAATGAACATCGAATATAACTAAGTTTCTGAACAACTCACGTTGCTGAg                                  <  2:372384‑M2/204‑1 (MQ=255)
                                                                                                                                                                                                                actcacttagcctcgtgtccactattgctgggtaagatCACCCAGTCATCCCCCGTAGAAATGAACTGCCCACAGAATGGTTCTATTCCGGGCAAAACCTACTCCAGCAAAGCGTTAATGAGCGGCAATGTCAAAAACGCGCAAATAGCCTCAGTAAAAGTGCAGTATCTTGATAAACAGAAAAAGCTGGCTGTAATGAATATCGAATAGAACTAAGTTTCTGAACAACTCACGTTGCTGAGCAGaaaaa                          <  1:179977‑M2/212‑1 (MQ=255)
                                                                                                                                                                                                                                gtccactattgctgggtaagatCACCCAGTCATCACCCGTAGAAATGAACTGCCCACAGAATGGTTCTATTCCGGGCAAAACCTACTCCAGCAAAGCGTTAATGAGCGGCAATGTCAAAAACGCGCAAATAGCCTCAGTAAAAGTGCAGTATCTTGATAAACAGAAAAAgctgg                                                                                      >  2:15911‑M2/23‑174 (MQ=255)
                                                                                                                                                                                                                                gtccactattgctgggtaagatCACCCAGTCATCACCCGTAGAAATGAACTGCCCACAGAATGGTTCTATTCCGGGCAAAACCTACTCCAGCAAAGCGTTAATGAGCGGCAATGTCAAAAACGCGCAAATAGCCTCAGTAAAAGTGCAGTATCTTGATAAACAGAAAAAgctgg                                                                                      <  1:15911‑M2/152‑1 (MQ=255)
                                                                                                                                                                                                                                       attgctgggtaagatCACCCAGTCATCACCCGTAGAAATGAACTGCCCACAGAATGGTTCTATTCCGGGCAAAACCTACTCCAGCAAAGCGTTAATGAGCGGCAATGTCAAAAACGCGCAAATAGCCTCAGTAAAAGTGCAGTATCTTGATAAACAGAAAAAGCTGGCTGTAATGAACATCGAATATAACTAAGTTTCTGAACAACTCACGTTGCGGAGCAGAAAAATGCGAGTTAACCAAAAAGCctgct  <  1:159721‑M2/236‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                           
AATATCGAGCATAAGATGAATCTGGAGAGAATGGTCTGCTGCGAATCAGCCAACCTGAAAGTATGGATAACACAACCCTCAAGGATGACTAATCATTGAGGAAATAGAATAAATGTTCAGACCTTTTTTAAACTCTCTTATGCTCGGCAGTTTGTTTTTTCCTTTTATTGCCATTGCTGGAAGCACCGTGCAAGGGGGCGTGATCCATTTTTATGGCCAAATTGTGGAACCGGCATGTGACGTCAGCACCCAGTCATCACCCGTAGAAATGAACTGCCCACAGAATGGTTCTATTCCGGGCAAAACCTACTCCAGCAAAGCGTTAATGAGCGGCAATGTCAAAAACGCGCAAATAGCCTCAGTAAAAGTGCAGTATCTTGATAAACAGAAAAAGCTGGCTGTAATGAACATCGAATATAACTAAGTTTCTGAACAACTCACGTTGCTGAGCAGAAAAATGCGATTTAACCAAAAAGCCTGCT  >  NC_000913/1105302‑1105783

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