Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 2,614,614 Δ3 bp coding (681‑683/849 nt) focB → putative formate transporter

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 2614614 2614616 3 42 [0] [0] 35 focB putative formate transporter

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 26146130 (0.000)35 (0.800) 26/212 0.5 100% coding (680/849 nt) focB putative formate transporter
?NC_000913 2614617 = 0 (0.000)coding (684/849 nt) focB putative formate transporter

GCGGGCGTTCGTTATGCGACAAAATCGTCGCCATGATTTTGCCCATCACCCTGTTTGTCGCCAGTGGCTTTGAGCACTGTATCGCCAATTTGTTTGTGATTCCGTTCGCCATTGCCATTCGCCATTTCGCCCCTCCCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATGCTGGGTAATATTATCGGCGGTGCGGTGCTGGTGAGTA  >  NC_000913/2614481‑2614744
                                                                                                                                        |                                                                                                                               
attgcgggcgttcgttatgcgacaaaatcgtcgccatgattttgcccatcaccctgtttgtcgccagtggctttgagcactgtatcgccaatttgtttgtgattccgttcgccattgccattcgccatttcgccccccc                                                                                                                               >  2:10323‑M2/137‑139 (MQ=255)
attgcgggcgttcgttatgcgacaaaatcgtcgccatgattttgcccatcaccctgtttgtcgccagtggctttgagcactgtatcgccaatttgtttgtgattccgttcgccattgccattcgccatttcgccccccc                                                                                                                               >  1:35045‑M2/137‑139 (MQ=255)
           tcgttatgcgacaaaatcgtcgccatgattttgcccatcaccctgtttgtcgccagtggctttgagcactgtatcgccaatttgtttgtgattccgttcgccattgccattcgccatttcgccccCCCCTTCTGGCAGc                                                                                                                    <  2:468864‑M2/14‑1 (MQ=255)
             gttatgcgacaaaatcgtcgccatgattttgcccatcaccctgtttgtcgccagtggctttgagcactgtatcgccaatttgtttgtgattccgttcgccattgccattcgccatttcgccccCCCCTTCTGGCAGCTg                                                                                                                  <  1:404755‑M2/16‑1 (MQ=255)
                atgcgacaaaatcgtcgccatgattttgcccatcaccctgtttgtcgccagtggctttgagcactgtatcgccaatttgtttgtgattccgttcgccattgccattcgccatttcgccccCCCCTTCTGGCAGCTGgcg                                                                                                               <  2:768818‑M2/19‑1 (MQ=255)
                     acaaaatcgtcgccatgattttgcccatcaccctgtttgtcgccagtggctttgagcactgtatcgccaatttgtttgtgattccgttcgccattgccattcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAg                                                                                                          >  1:537632‑M2/116‑139 (MQ=255)
                               cgccatgattttgcccatcaccctgtttgtcgccagtggctttgagcactgtatcgccaatttgtttgtgattccgttcgccattgccattcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGAc                                                                                                <  2:35045‑M2/34‑1 (MQ=255)
                                     gattttgcccatcaccctgtttgtcgccagtggctttgagcactgtatcgccaatttgtttgtgattccgttcgccattgccattcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAAtttt                                                                                          >  1:233051‑M2/100‑139 (MQ=255)
                                                       gtttgtcgccagtggctttgagcactgtatcgccaatttgtttgtgattccgttcgccattgccattcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGc                                                                        <  2:351783‑M2/58‑1 (MQ=255)
                                                           gtcgccagtggctttgagcactgtatcgccaatttgtttgtgattccgttcgccattgccattcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCAtt                                                                    >  2:693481‑M2/78‑139 (MQ=255)
                                                              gccagtggctttgagcactgtatcgcctatttgtttgtgattccgttcgccattgccattcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGc                                                                        >  1:601779‑M2/75‑132 (MQ=255)
                                                              gccagtggctttgagcactgtatcgccaatttgtttgtgattccgttcgccattgccattcgccatttcgccccccc                                                                                                                               <  1:706091‑M2/3‑1 (MQ=255)
                                                              gccagtggctttgagcactgtatcgccaatttgtttgtgattccgttcgccattgccattcgccatttcgccccccc                                                                                                                               >  2:706091‑M2/75‑77 (MQ=255)
                                                              gccagtggctttgagcactgtatcgccaatttgtttgtgattccgttcgccattgccattcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGc                                                                        <  2:601779‑M2/58‑1 (MQ=255)
                                                                                tgtatcgccaatttgtttgtgattccgttcgccattgccattcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCAtt                                                                    >  1:281091‑M2/57‑118 (MQ=255)
                                                                                tgtatcgccaatttgtttgtgattccgttcgccattgccattcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCAtt                                                                    <  2:281091‑M2/62‑1 (MQ=255)
                                                                                  tatcgccaatttgtttgtgattccgttcgccattgccattcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAAt                                                      <  1:781692‑M2/76‑1 (MQ=255)
                                                                                  tatcgccaatttgtttgtgattccgttcgccattgccattcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAAt                                                      >  2:781692‑M2/55‑130 (MQ=255)
                                                                                     cgccaatttgtttgtgattccgttcgccattgccattcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTg                                          >  1:418017‑M2/52‑139 (MQ=255)
                                                                                     cgccaatttgtttgtgattccgttcgccattgccattcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGa                                         >  1:141360‑M2/52‑139 (MQ=255)
                                                                                       ccaatttgtttgtgattccgttcgccattgccattcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATg                                       <  2:141360‑M2/90‑1 (MQ=255)
                                                                                        caatttgtttgtgattccgttcgccattgccattcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATg                                       <  1:26492‑M2/91‑1 (MQ=255)
                                                                                        caatttgtttgtgattccgttcgccattgccattcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATg                                       <  2:708277‑M2/91‑1 (MQ=255)
                                                                                        caatttgtttgtgattccgttcgccattgccattcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATg                                       <  2:418017‑M2/91‑1 (MQ=255)
                                                                                                              gccattgccattcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATGCTGGGTAATATTATCGgcggtg                 <  2:233051‑M2/113‑1 (MQ=255)
                                                                                                                      cattcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATGCTGGGTAATATTATCGGCGGTGCGGTGCTg         <  2:42061‑M2/121‑1 (MQ=255)
                                                                                                                      cattcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATGCTGGGTAATATTATCGGCGGTGCGGTGCTg         <  1:693481‑M2/121‑1 (MQ=255)
                                                                                                                       attcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATGCTGGGTAATATTATCGGCGGTGCGGTGCTgg        <  2:36480‑M2/122‑1 (MQ=255)
                                                                                                                        ttcgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATGCTGGGTAATATTATCGGCGGTGCGGTGCTGGt       <  2:993672‑M2/123‑1 (MQ=255)
                                                                                                                          cgccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATGCTGGGTAATATTATCGGCGGTGCGGTGCTGGTGa     >  1:480665‑M2/15‑139 (MQ=255)
                                                                                                                            ccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATGCTGGGTAATATTATCGGCGGTGCGGTGCTg         >  1:958969‑M2/13‑133 (MQ=255)
                                                                                                                            ccatttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATGCTGGGTAATATTATCGGCGGTGCGGTGCTg         <  2:958969‑M2/121‑1 (MQ=255)
                                                                                                                             catttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATGCTGGGTAATATTATCGGCGGTGCGGTGCTGGTGAGTa  <  2:260395‑M2/128‑1 (MQ=255)
                                                                                                                                ttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAgcg                                                                                                     >  2:500829‑M2/9‑37 (MQ=255)
                                                                                                                                ttcgccccCCCCTTCTGGCAGCTGGCGCACAGTAgcg                                                                                                     <  1:500829‑M2/29‑1 (MQ=255)
                                                                                                                                        |                                                                                                                               
GCGGGCGTTCGTTATGCGACAAAATCGTCGCCATGATTTTGCCCATCACCCTGTTTGTCGCCAGTGGCTTTGAGCACTGTATCGCCAATTTGTTTGTGATTCCGTTCGCCATTGCCATTCGCCATTTCGCCCCTCCCCCCTTCTGGCAGCTGGCGCACAGTAGCGCAGACAATTTTCCGGCACTGACGGTCAGCCATTTTATTACCGCCAATCTGCTCCCGGTGATGCTGGGTAATATTATCGGCGGTGCGGTGCTGGTGAGTA  >  NC_000913/2614481‑2614744

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ 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.