Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 1,367,750 Δ176 bp [pspF] [pspF]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 1367750 1367925 176 21 [0] [0] 22 [pspF] [pspF]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 13677490 (0.000)21 (0.500) 21/480 0.7 100% coding (164/978 nt) pspF DNA‑binding transcriptional dual regulator PspF
?NC_000913 1367926 = 0 (0.000)intergenic (‑14/‑153) pspF/pspA DNA‑binding transcriptional dual regulator PspF/phage shock protein A

GTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATAAAAAATTGGCACGCAAATTGTATTAACAGTTCAGCAGGACAATCCTGAACGCAGAAATCAAGAGGACAACATTATGGGTATTTTTTCTCGCTTTGCCGACATCGTGAATGCCAACATCAACGCTCTGTTAGAGAAAGCGGAAGATCCACAGAAACTGGTTCGTCTGATGAT  >  NC_000913/1367684‑1368176
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                          
gctgccgccaacgcgctccagttcaccgtactcaatcacgcgcaataatttctcctgcaccatcatcggtgccgtagcgagttcatcaagaaatagcgtaccgccgtcggcacgttcaaatctccctggatgacgtttttgcgcaccggtaaacgcccccgcttcgtgaccaaacagttcggaatccagcagattttcatttaacgccgcgcagttaagggaaataaacggcccttgccaacACTTGTTAg                                                                                                                                                                                                                                                    <  2:417787‑M2/9‑1 (MQ=255)
                 ccagttcaccgtactcaatcacgcgcaataatttctcctgcaccatcatcggtgccgtagcgagttcatcaagaaatagcgtaccgccgtcggcacgttcaaatctccctggatgacgtttttgcgcaccggtaaacgcccccgcttcgtgaccaaacagttcggaatccagcagattttcatttaacgccgcgcagttaagggaaataaacggcccttgccaacACTTGTTAGTGTAATTCGCTAACTCa                                                                                                                                                                                                                                   <  1:280704‑M2/26‑1 (MQ=255)
                        accgtactcaatcacgcgcaataatttctcctgcaccatcatcggtgccgtagcgagttcatcaagaaatagcgtaccgccgtcggcacgttcaaatctccctggatgacgtttttgcgcaccggtaaacgcccccgcttcgtgaccaaacagttcggaatccagcagattttcatttaacgccgcgcagttaagggaaataaacggcccttgccaacACTTGTTAGTGTAATTCGCTAACTCATCCTGGc                                                                                                                                                                                                                            >  2:72396‑M2/219‑251 (MQ=255)
                                caatcacgcgcaataatttctcctgcaccatcatcggtgccgtagcgagttcatcaagaaatagcgtaccgccgtcggcacgttcaaatctccctggatgacgtttttgcgcaccggtaaacgcccccgcttcgtgaccaaacagttcggaatccagcagattttcatttaacgccgcgcagttaagggaaataaacggcccttgccaacACTTGTTAGTGTAATTCGCTAACTCATCATGGCATGTTGCt                                                                                                                                                                                                                    <  1:433924‑M2/41‑1 (MQ=255)
                                 aatcacgcgcaataatttctcctgcaccatcatcggtgccgtagcgagttcatcaagaaatagcgtaccgccgtcggcacgttcaaatctccctggatgacgtttttgcgcaccggtaaacgcccccgcttcgtgaccaaacagttcggaatccagcagattttcatttaacgccgcgcagttaagggaaataaacggcccttgccaacACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTg                                                                                                                                                                                                                   >  1:414263‑M2/210‑251 (MQ=255)
                                  atcacgcgcaataatttctcctgcaccatcatcggtgccgtagcgagttcatcaagaaatagcgtaccgccgtcggcacgttcaaatctccctggatgacgtttttgcgcaccggtaaacgcccccgcttcgtgaccaaacagttcggaatccagcagattttcatttaacgccgcgcagttaagggaaataaacggcccttgccaacACTTGTTATTGTAATTCGCTAACTCTGCCTGGCATGTTGCTGt                                                                                                                                                                                                                  <  1:133596‑M2/43‑1 (MQ=255)
                                      cgcgcaataatttctcctgcaccatcatcggtgccgtagcgagttcatcaagaaatagcgtaccgccgtcggcacgttcaaatctccctggatgacgtttttgcgcaccggtaaacgcccccgcttcgtgaccaaacagttcggaatccagcagattttcatttaacgccgcgcagttaagggaaataaacggcacttgccaacACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGAt                                                                                                                                                                                                              >  1:91855‑M2/205‑251 (MQ=255)
                                                  tctcctgcaccatcatcggtgccgtagcgagttcatcaagaaatagcgtaccgccgtcggcacgttcaaatctccctggatgacgtttttgcgcaccggtaaacgcccccgcttcgtgaccaaacagttcggaatccagcagattttcatttaacgccgcgcagttaagggaaataaacggcccttgccaacACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAg                                                                                                                                                                                                   <  1:338308‑M2/58‑1 (MQ=255)
                                                       tgcaccatcatcggtgccgtagcgagttcatcaagaaatagcgtaccgccgtcggcacgttcaaatctccctggatgacgtttttgcgcaccggtaaacgcccccgcttcgtgaccaaacagttcggaatccagcagattttcatttaacgccgcgcagttaagggaaataaacggcccttgccaacACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCttt                                                                                                                                                                                             >  1:37077‑M2/188‑251 (MQ=255)
                                                           ccatcagcggtgccgtagcgagttcatcaagaaatagcgtaccgccgtcggcacgttcaaatctccctggatgacgtttttgcgcaccggtaaacgcccccgcttcgtgaccaaacagttcggaatccagcagattttcatttaacgccgcgcagttaagggaaataaacggcccttgccaacACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTCCAATCAGATCTTTATaa                                                                                                                                                                                         <  2:414263‑M2/68‑1 (MQ=255)
                                                               catcggtgccgtagcgagttcatcaagaaatagcgtaccgccgtcggcacgttcaaatctccctggatgacgtttttgcgcaccggtaaacgcccccgcttcgtgaccaaacagttcggaatccagcagattttcatttaacgccgcgcagttaagggaaataaacggcccttgccaacACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATc                                                                                                                                                                                      <  1:234838‑M2/71‑1 (MQ=255)
                                                                                         gaaatagcgtaccgccgtcggcacgttcaaatctccctggatgacgtttttgcgcaccggtaaacgcccccgcttcgtgaccaaacagttcggaatccagcagattttcatttaccgccgcgcagttaagggaaatacacggcccttgccaacACTTGTTCGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTCTAAATCAGAAAGATAAAAAATTGGCACGCAAAt                                                                                                                                                           >  2:113219‑M2/154‑251 (MQ=255)
                                                                                           aatagcgtaccgacggcggcacgttcaaatctccctggatgactgttttgcgccccggtaaccgcccccgcttcgtgaccaaacagttcggaatccagcagattttcatttaacgccgcgcagttaagggaaataaacggcccttgccatcACTTGTTAGTGTAATTCGCTAACTCACCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATAAAAAATTGGCACGCAAATTg                                                                                                                                                         <  1:402319‑M2/100‑1 (MQ=255)
                                                                                                             gcacgttcaaatctccctggatgacgtttttgcgcaccggtaaacgcccccgcttcgtgaccaaacagttcggaatccagcagattttcatttaacgccgcgcagttaagggaaataaacggcccttgccaacACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATAAAAAATTGGCACGCAAATTGTATTAACAGTTCAGCAgg                                                                                                                                       <  1:54140‑M2/118‑1 (MQ=255)
                                                                                                                            cctggatgacgtttttgcgcaccggtaaacgcccccgcttcgtgaccaaacagttcggaatccagcagattttcatttaacgccgcgcagttaagggaaataaacgtcccttgccaacACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATAAAAAATTGGCACGCAAATTGTATTAACAGTTCAGCAGGACAATCCTGAACGCa                                                                                                                        <  1:72396‑M2/133‑1 (MQ=255)
                                                                                                                                           tgcgcaccggtaaacgcccccgcttcgtgaccaaacaggtcggaatccagcagattttcatttaacgccgcgcagttaagggaaataaacggcccttgccaacACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATAAAAAATTGGCACGCAAATTGTATTAACAGTTCAGCAGGACAATCCTGAACGCAGAAATCAAGAGGac                                                                                                          <  2:255462‑M2/147‑1 (MQ=255)
                                                                                                                                                      aaacgcccccgcttcgtgaccaaacagttcggaatccagcagattttcatttaacgccgcgcagttaagggaaataaacggcccttgccaacACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATAAAAAATTGGCACGCAAATTGTATTAACAGTTCAGCAGGACAATCCTGAACGCAGAAATCAAGAGGACAACATTATGGGt                                                                                              <  1:113219‑M2/159‑1 (MQ=255)
                                                                                                                                                                                                      atttaacgccgcgcagttaagggaaataaacggcccttgccaacACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATAAAAAATTGGCACGCAAATTGTATTAACAGTTCAGCAGGACAATCCTGAACGCAGAAATCAAGAGGACAACATTATGGGTATTTTTTCTCGCTTTGCCGACATCGTGAATGCCAACATCAACGCTCTg                                              <  1:112618‑M2/207‑1 (MQ=255)
                                                                                                                                                                                                                gcgcagttaagggaaataaacggcccttgccaacACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATAAAAAATTGGCACGCAAATTGTATTAACAGTTCAGCAGGACAATCCTGAACGCAGAAATCAAGAGGACAACATTATGGGTATTTTTTCTCGCTTTGCCGACATCGTGAATGCCAACATCAACGCTCTGTTAGAGaa                                      >  2:106643‑M2/35‑249 (MQ=255)
                                                                                                                                                                                                                                          ttgccaacACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATAACAAATTGGCACGCAAATTGTATTAACAGTTCAGCAGGACAATCCTGAACGCAGAAATCAAGAGGACAACATTATGGGTATTTTTTCTCGCTTTGCCGACATc                                                                      >  2:145372‑M2/9‑191 (MQ=255)
                                                                                                                                                                                                                                          ttgccaacACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATAAAAAATTGGCACGCAAATTGTATTAACAGTTCAGCAGGACAATCCTGAACGCAGAAATCAAGAGGACAACAGTATGGGTATTTTTTCTCGCTTTGCCGACATc                                                                      <  1:145372‑M2/183‑1 (MQ=255)
                                                                                                                                                                                                                                                  aCTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCGGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATAAAAAATTGGCACGCAAATTGTATTAACAGTTCAGCAGGACAATCCTGAACGCAGAAATCAAGAGGACAACATTATGGGTATTTTTTCTCGCTTTGCCGACATCGTGAATGCCAACATCAACGCTCTGTTAGAGAAAGCGGAAGATCCACAGAAACTGGTTCGTCTgatgat  <  1:267121/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                          
GTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTGCTGTTGATTCTTCAATCAGATCTTTATAAATCAAAAAGATAAAAAATTGGCACGCAAATTGTATTAACAGTTCAGCAGGACAATCCTGAACGCAGAAATCAAGAGGACAACATTATGGGTATTTTTTCTCGCTTTGCCGACATCGTGAATGCCAACATCAACGCTCTGTTAGAGAAAGCGGAAGATCCACAGAAACTGGTTCGTCTGATGAT  >  NC_000913/1367684‑1368176

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.