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

GCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTG  >  NC_000913/1367508‑1367964
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                       
gCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttag                                                                                                                                                                                                                <  2:417787‑M1/251‑10 (MQ=255)
                 ccAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactca                                                                                                                                                                                               <  1:280704‑M1/251‑27 (MQ=255)
                        aCCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggc                                                                                                                                                                                        >  2:72396‑M1/1‑218 (MQ=255)
                                cAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcatggcatgttgct                                                                                                                                                                                <  1:433924‑M1/251‑42 (MQ=255)
                                 aaTCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctg                                                                                                                                                                               >  1:414263‑M1/1‑209 (MQ=255)
                                  aTCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttattgtaattcgctaactctgcctggcatgttgctgt                                                                                                                                                                              <  1:133596‑M1/251‑44 (MQ=255)
                                      cgcgCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCACTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgat                                                                                                                                                                          >  1:91855‑M1/1‑204 (MQ=255)
                                                  tctcCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcag                                                                                                                                                               <  1:338308‑M1/250‑59 (MQ=255)
                                                       tGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatcttt                                                                                                                                                         >  1:37077‑M1/1‑187 (MQ=255)
                                                           ccATCAGCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattctccaatcagatctttataa                                                                                                                                                     <  2:414263‑M1/251‑69 (MQ=255)
                                                               catcGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatctttataaatc                                                                                                                                                  <  1:234838‑M1/250‑72 (MQ=255)
                                                                                         gAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTACCGCCGCGCAGTTAAGGGAAATACACGGCCCTTGCCAACacttgttcgtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatctttctaaatcagaaagataaaaaattggcacgcaaat                                                                                                                       >  2:113219‑M1/1‑153 (MQ=255)
                                                                                           aaTAGCGTACCGACGGCGGCACGTTCAAATCTCCCTGGATGACTGTTTTGCGCCCCGGTAACCGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCATCacttgttagtgtaattcgctaactcaccctggcatgttgctgttgattcttcaatcagatctttataaatcaaaaagataaaaaattggcacgcaaattg                                                                                                                     <  1:402319‑M1/251‑101 (MQ=255)
                                                                                                             gCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatctttataaatcaaaaagataaaaaattggcacgcaaattgtattaacagttcagcagg                                                                                                   <  1:54140‑M1/251‑119 (MQ=255)
                                                                                                                            ccTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGTCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatctttataaatcaaaaagataaaaaattggcacgcaaattgtattaacagttcagcaggacaatcctgaacgca                                                                                    <  1:72396‑M1/251‑134 (MQ=255)
                                                                                                                                           tGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGGTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatctttataaatcaaaaagataaaaaattggcacgcaaattgtattaacagttcagcaggacaatcctgaacgcagaaatcaagaggac                                                                      <  2:255462‑M1/250‑148 (MQ=255)
                                                                                                                                                      aaaCGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatctttataaatcaaaaagataaaaaattggcacgcaaattgtattaacagttcagcaggacaatcctgaacgcagaaatcaagaggacaacattatgggt                                                          <  1:113219‑M1/251‑160 (MQ=255)
                                                                                                                                                                                                      aTTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatctttataaatcaaaaagataaaaaattggcacgcaaattgtattaacagttcagcaggacaatcctgaacgcagaaatcaagaggacaacattatgggtattttttctcgctttgccgacatcgtgaatgccaacatcaacgctctg          <  1:112618‑M1/251‑208 (MQ=255)
                                                                                                                                                                                                                gcgcAGTTAAGGGAAATAAACGGCCCTTGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatctttataaatcaaaaagataaaaaattggcacgcaaattgtattaacagttcagcaggacaatcctgaacgcagaaatcaagaggacaacattatgggtattttttctcgctttgccgacatcgtgaatgccaacatcaacgctctgttagagaa  >  2:106643‑M1/1‑34 (MQ=255)
                                                                                                                                                                                                                                          ttGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatctttataaatcaaaaagataacaaattggcacgcaaattgtattaacagttcagcaggacaatcctgaacgcagaaatcaagaggacaacattatgggtattttttctcgctttgccgacatc                                  >  2:145372‑M1/1‑8 (MQ=255)
                                                                                                                                                                                                                                          ttGCCAACacttgttagtgtaattcgctaactcatcctggcatgttgctgttgattcttcaatcagatctttataaatcaaaaagataaaaaattggcacgcaaattgtattaacagttcagcaggacaatcctgaacgcagaaatcaagaggacaacagtatgggtattttttctcgctttgccgacatc                                  <  1:145372‑M1/191‑184 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                       
GCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAACGCCCCCGCTTCGTGACCAAACAGTTCGGAATCCAGCAGATTTTCATTTAACGCCGCGCAGTTAAGGGAAATAAACGGCCCTTGCCAACGGGAGGAGAGATAATGCAGGCGGCTGGCAATCAGCTCTTTACCGGTGCCGCGTTCGCCGATGATGAGCACCGGTTTGTCCAGCGGTGCGAGATGCGAAACCTGTTCCAGCACTTCGAGAAAGCTGTTCGCCTCACCAAGTAAATTATCTTTGTATTCTGCCATGATGAAATTCGCCACTTGTTAGTGTAATTCGCTAACTCATCCTGGCATGTTG  >  NC_000913/1367508‑1367964

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 34 ≤ 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.