Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NZ_CP009273 2,399,465 Δ7 bp coding (650‑656/939 nt) lrhA ← transcriptional regulator LrhA

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NZ_CP009273 2399465 2399471 7 18 [0] [0] 18 lrhA transcriptional regulator LrhA

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 23994640 (0.000)18 (0.890) 13/166 0.4 100% coding (657/939 nt) lrhA transcriptional regulator LrhA
?NZ_CP009273 2399472 = 0 (0.000)coding (649/939 nt) lrhA transcriptional regulator LrhA

GTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATC  >  NZ_CP009273/2399383‑2399544
                                                                                 |                                                                                
gTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagt                                                                          >  2:384081‑M1/1‑82 (MQ=255)
    gggCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgcc                                                                      >  1:294511‑M1/1‑78 (MQ=255)
    gggCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgcc                                                                      <  2:155997‑M1/90‑13 (MQ=255)
      gCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccat                                                                    <  2:508538‑M1/90‑15 (MQ=255)
              tctcAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatc                                                            <  1:269410‑M1/90‑23 (MQ=255)
                   aCCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaaaatctgctt                                                       >  1:196637‑M1/1‑63 (MQ=255)
                      ggCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttat                                                    >  1:447277‑M1/1‑60 (MQ=255)
                            gCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcg                                              <  1:230840‑M1/90‑37 (MQ=255)
                            gCCGTCACGCCAAGACCGGCTTTCACTACCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcg                                              <  2:210561‑M1/90‑37 (MQ=255)
                                 cACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgcc                                         <  2:90784‑M1/90‑42 (MQ=255)
                                             ggCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcg                             >  2:37310‑M1/1‑37 (MQ=255)
                                             ggCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcg                             >  1:477719‑M1/1‑37 (MQ=255)
                                                   cACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaac                       <  1:562683‑M1/90‑60 (MQ=255)
                                                                     cGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcatccagc     <  2:447277‑M1/90‑78 (MQ=255)
                                                                     cGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcatccagc     <  2:500593‑M1/90‑78 (MQ=255)
                                                                     cGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcatccagc     <  1:427461‑M1/90‑78 (MQ=255)
                                                                        aaGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcatccagcaat  <  2:294511‑M1/90‑81 (MQ=255)
                                                                        aaGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcatccagcaat  <  2:332768‑M1/90‑81 (MQ=255)
                                                                                 |                                                                                
GTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATC  >  NZ_CP009273/2399383‑2399544

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