Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 4,619,478 Δ9 bp coding (1250‑1258/1323 nt) deoA → thymidine phosphorylase

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 4619478 4619486 9 67 [0] [0] 67 deoA thymidine phosphorylase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 46194770 (0.000)64 (1.730) 47/222 0.0 100% coding (1249/1323 nt) deoA thymidine phosphorylase
?NC_000913 4619487 = 0 (0.000)coding (1259/1323 nt) deoA thymidine phosphorylase

CATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTAT  >  NC_000913/4619356‑4619622
                                                                                                                                   |                                                                                                                                       
atctgacaccatcgattacagcgtcggctttactgatatggcgcgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTaaa                                                                                                                                  <  1:94488‑M2/8‑1 (MQ=255)
 tctgacaccatcgattacagcgtcggctttactgatatggcgcgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAAc                                                                                                                                 <  2:34712‑M2/9‑1 (MQ=255)
            cgattacagcgtcggctttactgatatggcgcgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATaaa                                                                                                                      >  1:368987‑M2/120‑139 (MQ=255)
             gattacagcgtcggctttactgatatggcgcgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAg                                                                                                                     >  2:171576‑M2/119‑139 (MQ=255)
              attacagcgtcggctttactgatatggcgcgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGc                                                                                                                    >  2:216112‑M2/118‑139 (MQ=255)
                 acagcgtcggctttactgatatggcgcgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCAcc                                                                                                                 >  1:598748‑M2/115‑139 (MQ=255)
                    gcgtcggctttactgatatggcgcgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGa                                                                                                              >  2:708887‑M2/112‑139 (MQ=255)
                     cgtcggctttactgatatggcgcgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGaa                                                                                                             >  1:426772‑M2/111‑139 (MQ=255)
                     cgtcggctttactgatatggcgcgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGaa                                                                                                             >  2:413857‑M2/111‑139 (MQ=255)
                       tcggctttactgatatggcgcgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAg                                                                                                           >  1:221314‑M2/109‑139 (MQ=255)
                         ggctttactgatatggcgcgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCaa                                                                                                                                       <  1:268413‑M2/3‑1 (MQ=255)
                         ggctttactgatatggcgcgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCaa                                                                                                                                       >  2:268413‑M2/107‑109 (MQ=255)
                         ggctttactgatatggcgcgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGca                                                                                                         >  2:723066‑M2/107‑139 (MQ=255)
                          gctttactgatatggcgcgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGcac                                                                                                        <  2:175059‑M2/34‑1 (MQ=255)
                               actgatatggcgcgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCaa                                                                                                   <  1:786438‑M2/39‑1 (MQ=255)
                                ctgatatggcgcgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAAc                                                                                                  <  2:617270‑M2/40‑1 (MQ=255)
                                         cgcgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCaa                                                                                                   <  1:255780‑M2/39‑1 (MQ=255)
                                         cgcgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCaa                                                                                                   >  2:255780‑M2/91‑129 (MQ=255)
                                         cgcgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCg                                                                                         >  1:87558‑M2/91‑139 (MQ=255)
                                           cgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGcc                                                                                       >  2:783113‑M2/89‑139 (MQ=255)
                                             tctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGACGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGt                                                                                     >  1:760937‑M2/87‑139 (MQ=255)
                                                      ccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGaa                                                                            <  1:626017‑M2/62‑1 (MQ=255)
                                                      ccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGaa                                                                            >  2:238649‑M2/78‑139 (MQ=255)
                                                       caggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACCTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAAt                                                                           <  1:480912‑M2/63‑1 (MQ=255)
                                                         ggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATaa                                                                         <  1:723066‑M2/65‑1 (MQ=255)
                                                         ggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATaa                                                                         <  2:760937‑M2/65‑1 (MQ=255)
                                                           tagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACg                                                                       <  1:216112‑M2/67‑1 (MQ=255)
                                                           tagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACg                                                                       <  2:614354‑M2/67‑1 (MQ=255)
                                                              acggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGta                                                                    >  1:392715‑M2/70‑139 (MQ=255)
                                                                 gtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCaa                                                                                                   >  2:756209‑M2/67‑105 (MQ=255)
                                                                 gtcagcgtccgctggcggtgatccacgcgaaagacgaaaacaactggcagtaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCaa                                                                                                   <  1:756209‑M2/39‑1 (MQ=255)
                                                                   cagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCa                                                                                 >  2:99151‑M2/65‑121 (MQ=255)
                                                                   cagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCa                                                                                 <  1:99151‑M2/57‑1 (MQ=255)
                                                                    agcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATaa                                                                         >  1:379317‑M2/64‑128 (MQ=255)
                                                                    agcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATaa                                                                         <  2:379317‑M2/65‑1 (MQ=255)
                                                                     gcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATActgat                                                             <  2:294740‑M2/77‑1 (MQ=255)
                                                                      cgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATActgatc                                                            >  2:698885‑M2/62‑139 (MQ=255)
                                                                            ctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCa                                                      >  1:428551‑M2/56‑139 (MQ=255)
                                                                              ggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCAtt                                                    <  2:598748‑M2/86‑1 (MQ=255)
                                                                                cggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGGCGATAAAGCACCGGAAAGCACACCAACTGTCGATCGCCGTATCAGCGAATAACGGCATACTGATCTGATCATTTa                                                  >  1:520388‑M2/52‑139 (MQ=255)
                                                                                 ggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTaa                                                 <  1:413857‑M2/89‑1 (MQ=255)
                                                                                        cacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGaa                                          <  2:684760‑M2/96‑1 (MQ=255)
                                                                                          cgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGc                                        <  2:392715‑M2/98‑1 (MQ=255)
                                                                                           gcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCa                                       <  1:238649‑M2/99‑1 (MQ=255)
                                                                                            cgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACg                                                                       >  1:221045‑M2/40‑106 (MQ=255)
                                                                                            cgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACg                                                                       <  2:221045‑M2/67‑1 (MQ=255)
                                                                                            cgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCAc                                      <  2:368987‑M2/100‑1 (MQ=255)
                                                                                            cgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCAc                                      <  2:87558‑M2/100‑1 (MQ=255)
                                                                                             gaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGaa                                                                            >  2:717212‑M2/39‑100 (MQ=255)
                                                                                             gaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGaa                                                                            <  1:717212‑M2/62‑1 (MQ=255)
                                                                                              aaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTg                                    <  1:649412‑M2/102‑1 (MQ=255)
                                                                                               aagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGa                                   <  1:165429‑M2/103‑1 (MQ=255)
                                                                                               aagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGa                                   <  2:428551‑M2/103‑1 (MQ=255)
                                                                                                    gaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACg                              <  2:524961‑M2/108‑1 (MQ=255)
                                                                                                       aacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGgag                           <  2:228329‑M2/111‑1 (MQ=255)
                                                                                                                    aagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGca                                                                                                         >  2:272637‑M2/16‑48 (MQ=255)
                                                                                                                    aagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGca                                                                                                         <  1:272637‑M2/33‑1 (MQ=255)
                                                                                                                       cggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGc           <  2:426772‑M2/127‑1 (MQ=255)
                                                                                                                         gcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGaa                                                                            <  1:348944‑M2/62‑1 (MQ=255)
                                                                                                                         gcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGaa                                                                            >  2:348944‑M2/11‑72 (MQ=255)
                                                                                                                         gcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGgag                           <  2:102908‑M2/111‑1 (MQ=255)
                                                                                                                         gcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGgag                           >  1:102908‑M2/11‑121 (MQ=255)
                                                                                                                         gcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGc           >  2:426887‑M2/11‑137 (MQ=255)
                                                                                                                         gcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGc           <  1:426887‑M2/127‑1 (MQ=255)
                                                                                                                              agcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGaa                                                                            >  2:581163‑M2/6‑67 (MQ=255)
                                                                                                                              agcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGaa                                                                            <  1:581163‑M2/62‑1 (MQ=255)
                                                                                                                                cggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTTGATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACCGAGTACGGAGAACATATGAAACGTGCATTtattat  <  2:242049‑M2/136‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                       
CATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTAT  >  NC_000913/4619356‑4619622

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