Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 3,909,311 Δ774 bp [pstA][pstC] [pstA], [pstC]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 3909311 3910084 774 50 [0] [0] 50 [pstA]–[pstC] [pstA], [pstC]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 39093100 (0.000)48 (0.950) 42/268 0.1 100% coding (129/891 nt) pstA phosphate ABC transporter permease
?NC_000913 3910085 = 0 (0.000)coding (313/960 nt) pstC phosphate ABC transporter permease

AATTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGAACAGGCCGAAGGCCATCGTCGCCATCGAGAGCGTCAGCGCAATACGGTTTTTGAGGCGGCGACGCGCCTGCATTTTGCGGCGAGATTCAGCCAGCGCCGCAGTGGTTTGCATTTCAACCATAGCCATTAGCGTGCCCCC  >  NC_000913/3909169‑3909450
                                                                                                                                             |                                                                                                                                            
aaTTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGACGGGCATCAAACTCCCGATCAGCCAGg                                                                                                                                             >  1:713767‑M1/1‑142 (MQ=255)
   ttacaGCCCGCTACCCGCCATAGCGTTCGCCAGACCAGCGCCTTCCGGATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgc                                                                                                                                          <  1:609451‑M1/141‑5 (MQ=255)
          ccGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGATATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAAgcgaaagctcg                                                                                                                                   >  1:466183‑M1/1‑132 (MQ=255)
                       gagCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacaggg                                                                                                                      <  1:753444‑M1/143‑25 (MQ=255)
                               gCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccg                                                                                                              >  1:628390‑M1/1‑111 (MQ=255)
                                       acccccTCCCGTATCGGGCGGCGGCGTCATTGCAGTGACCAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcac                                                                                                      <  2:658207‑M1/143‑41 (MQ=255)
                                          accTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcac                                                                                                      >  2:456732‑M1/1‑100 (MQ=255)
                                          accTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcac                                                                                                      <  1:456732‑M1/140‑41 (MQ=255)
                                            cTTCCGTATTGGGCGGCGGCGTCAGTTCAGGGAGCGGCTCCAGCGACATACCGGCGATACCGCGAGTGATGGTGGACATTAAAAGCCGGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccggaactcgccggga                                                                                                 <  2:621938‑M1/143‑46 (MQ=255)
                                             ttCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggac                                                                                                >  1:93202‑M1/1‑97 (MQ=255)
                                             ttCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagcgcagtcaggaccaggacgataccgcacctcaccgcgac                                                                                                >  2:649191‑M1/1‑97 (MQ=255)
                                               tcctattGGGCGGTGGCGTTTTTTTAGTGAGCTGCGCCAGCGACCTACCGTCGATACCGCGAGTGATGTTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtaaagaacaggcctatacctaaactcaccgggacgg                                                                                              <  2:332218‑M1/140‑49 (MQ=255)
                                                 gTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcg                                                                                            >  1:909641‑M1/1‑93 (MQ=255)
                                                  tagcgggcggcGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGGCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcga                                                                                           <  1:407049‑M1/139‑52 (MQ=255)
                                                  tATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcga                                                                                           >  1:101646‑M1/1‑92 (MQ=255)
                                                   atggggcggCGGCGTCATTCCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCTAGgcgcaagctcagtcaggaacagggcgataccgaagctcaccgggacggcgat                                                                                          <  1:484574‑M1/140‑53 (MQ=255)
                                                   aTTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgat                                                                                          <  1:623220‑M1/143‑53 (MQ=255)
                                                      gggCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcag                                                                                       <  1:561916‑M1/143‑56 (MQ=255)
                                                           gcggcgTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcg                                                                                  <  1:812032‑M1/143‑61 (MQ=255)
                                                             ggcgTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcg                                                                                <  1:473383‑M1/143‑63 (MQ=255)
                                                                    tggcagtGAGCAGCGCCAGCGACATACCGGCGATGCCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCAAGgcgcaagctcagtcaggaacagggcgataccgaaactaaccgggacggcgatcagcagcgcgataaacg                                                                         <  2:253131‑M1/140‑70 (MQ=255)
                                                                      tCAGTGAACAGCGCCAGCGACAGACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgagactcaccgggacggcgatcagcagcgcgataaacgaa                                                                       <  1:818682‑M1/143‑72 (MQ=255)
                                                                         gTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaaactcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgatcgacgaagtc                                                                    >  2:756123‑M1/1‑69 (MQ=255)
                                                                          ttacccgcgccagcgaCAGACCGGCGATACCGCGATCGAGTGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggtacagggcgatatcgaaactcaccgtgacggcgatcagcagcgcgataaacgaagtca                                                                   <  1:103759‑M1/137‑76 (MQ=255)
                                                                          tGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtca                                                                   <  1:248148‑M1/143‑76 (MQ=255)
                                                                            aaCAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcacccggaaggcgttcagcctcgcgataaacgaagtcacc                                                                 >  1:722260‑M1/1‑66 (MQ=255)
                                                                              cAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaa                                                               >  2:527728‑M1/1‑64 (MQ=255)
                                                                               aGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacaaggcgataattaaactcaccggcacggagcacagcagcgcaagcaaagaagtaaccaac                                                              >  2:107004‑M1/1‑63 (MQ=255)
                                                                                   ccAGCGACATACCGTCGAGACCGCGAGTGAGGGTGGACATTAAAAGCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacgacgcaccgccgcgcgataaacgaagtcaccaacgtac                                                          <  2:566187‑M1/143‑85 (MQ=255)
                                                                                      gCGACATACCGCCGATACCGCGAGTGACGGTGGACATTAAAATCCAGTTCAGCCAGccgcaagctcagtcaggaacagggcgataccgacactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgt                                                       <  1:110940‑M1/143‑88 (MQ=255)
                                                                                            tACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagcgcagtcaggaacagtgcgataccgaaactcaccgggacggcgatcagcatcgcgataaacgaagtcaccgacgtaccgtagatag                                                 >  1:444306‑M1/1‑50 (MQ=255)
                                                                                             aCCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcgg                                                <  2:101646‑M1/143‑95 (MQ=255)
                                                                                              ccGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggctataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggc                                               >  1:567397‑M1/1‑48 (MQ=255)
                                                                                                    aTACCGCGAGTGATGGCGGACATTAAAATCCAGAACAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagc                                         >  1:708962‑M1/1‑42 (MQ=255)
                                                                                                        cgcgGGGGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccc                                     <  2:470108‑M1/143‑106 (MQ=255)
                                                                                                           gAGTCATGAGGGACAGTACACGCCAGACCAGCCACgcgcaagctcagtcaggaacagcgcggtaccgaacgtcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgt                                  <  2:505363‑M1/143‑109 (MQ=255)
                                                                                                            aGTGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgta                                 <  1:102536‑M1/143‑110 (MQ=255)
                                                                                                              tGATGGTGGACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaacgcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtacggtacatcggcaccagcgcccggtgga                               >  2:675023‑M1/1‑32 (MQ=255)
                                                                                                                  tgtgtcccttacaatccaGCTCCGCCACgcgcaagctccgtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatc                           <  1:878354‑M1/135‑116 (MQ=255)
                                                                                                                    tggACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtcgatcggcaccagcgccccgtagctatcgt                         >  1:181914‑M1/1‑26 (MQ=255)
                                                                                                                    tggACATTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcaataccaaaacacaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgt                         >  1:264276‑M1/1‑26 (MQ=255)
                                                                                                                       aCATTAGAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcg                      <  2:901907‑M1/143‑121 (MQ=255)
                                                                                                                        ccttaaAATCCAGATCGCCCAGgcgcaggcgcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcgg                     <  2:131351‑M1/141‑122 (MQ=255)
                                                                                                                         aTTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaa                                                                       >  2:335040‑M1/1‑21 (MQ=255)
                                                                                                                         aTTAAAATCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaa                                                                       <  1:335040‑M1/92‑72 (MQ=255)
                                                                                                                         aTTAAAACCCAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggt                    <  1:536190‑M1/143‑123 (MQ=255)
                                                                                                                                 ccAGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatccca            >  1:365145‑M1/1‑13 (MQ=255)
                                                                                                                                   aGATCAGCCAGgcgcaagctcagtcaggaacagggcgataccgaaactcaccgggacggcgatcagcagcgcgataaacgcagtcaccaacgtcgcgtagatgggcaccagcgccccgtagaaattatttggtgcatcacgcg          >  1:465988‑M1/1‑11 (MQ=255)
                                                                                                                                          ccAGgcgcacgcccggccaggaacagagcgataccgtaactcaccgggacggcgatcagcagcgcgataaacgacgtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccactctttggt   <  2:592993‑M1/143‑140 (MQ=255)
                                                                                                                                           cAGgcgcaagctcagtcaggaacagggcgataccgaacctcacccgggcggcgatcagcagcgcgataaacgaagtcaccaacgtaccgtagatcggcaccagcgccccgtagatatcgttcggtgcatcccactctttggtc  <  1:253601‑M1/143‑141 (MQ=255)
                                                                                                                                             |                                                                                                                                            
AATTAACAGCCCGCTACCCGCCAGAGCGTTCGCCAGACCACCACCTTCCGTATTGGGCGGCGGCGTCATTTCAGTGAACAGCGCCAGCGACATACCGTCGATACCGCGAGTGATGGTGGACATTAAAATCCAGATCAGCCAGAACAGGCCGAAGGCCATCGTCGCCATCGAGAGCGTCAGCGCAATACGGTTTTTGAGGCGGCGACGCGCCTGCATTTTGCGGCGAGATTCAGCCAGCGCCGCAGTGGTTTGCATTTCAACCATAGCCATTAGCGTGCCCCC  >  NC_000913/3909169‑3909450

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: