Predicted mutation
evidence seq id position mutation annotation gene description
MC JC 1 1,420,817 (TGCCT)6→5 intergenic (+189/‑63) unknown → / → unknown Threonyl‑tRNA synthetase (EC 6.1.1.3)/Type III restriction‑modification system methylation subunit (EC 2.1.1.72)

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ 1 1420790 1420794 5 40 [0] [0] 40 unknown/unknown Threonyl‑tRNA synthetase (EC 6.1.1.3)/Type III restriction‑modification system methylation subunit (EC 2.1.1.72)

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? 1 1420795 =0 (0.000)36 (0.470) 30/376 1.3 100% intergenic (+167/‑89) unknown/unknown Threonyl‑tRNA synthetase (EC 6.1.1.3)/Type III restriction‑modification system methylation subunit (EC 2.1.1.72)
?1 = 1420789 0 (0.000)intergenic (+161/‑95) unknown/unknown Threonyl‑tRNA synthetase (EC 6.1.1.3)/Type III restriction‑modification system methylation subunit (EC 2.1.1.72)

CGTGGAAGAATTTTTAGACATTCTGAAAAAACAAGTGAGAAATCGCGAATTGAAATTATTGGGTGAGGCGTAGTTCTGACAATTAAATGAGCTGAAATTACGCTAGCAGGAATCATCATTCAATAGATTTATTGATTTTTCCGACCGGCGTAGTTTCAGCTTTTTAACACTCGATATTGTATAAAACAAAGTACACTTTGTACTTGACAACCAGTAATAGGAGTAAAAAATATGCTTGCCTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCGGAGGATAAATGATCTCCGAAGAACTTTTCTCTCGTGACGGCGTTTCCGCTAACTCTGCACAACTTACCCAACTTACCCAACTTAAAGCTCTTTTCCCGAATTGTTTCGATGTAAGCGGTCAATTTTTATTAGAAAAATTCCAATCTGAAATTGCAGAACAAGCC  >  1/1420556‑1421039
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    
ttcaccgtggaagaatttttagacattctgaaaaaacaagtgagaaatcgcgaattgaaattattgggtgaggcgtagttctgacaattaaatgagctgaaattacgctagcaggaatcatcattcaatagatttattgatttttccgaccggcgtagtttcagctttttaacactcgatattgtataaaacaaagtacactttgtacttgacaaccattaataggagtaaaaaatatgcttgccttgcct                                                                                                                                                                                                                                           >  1:483973‑M1/240‑251 (MQ=255)
                      acattctgaaaaaacaagtgagaaatcgcgaattgaaattattgggtgaggcgtagttctgacaattaaatgagctgaaattacgctagcaggaatcatcattcaatagatttattgatttttccgaccggcgtagtttcagctttttaacactcgatattgtataaaacaaagtacactttgtacttgacaaccagtaataggagtaaaaaatatgCTTGCCTTGCCTTGCCTTGCCTTGCctaaacta                                                                                                                                                                                                                      >  1:109985‑M1/218‑250 (MQ=255)
                                                            ttattgggtgaggcgtagttctgacaattaaatgagctgaaattacgctagcaggaatcatcattcaatagatttattgatttttccgaccggcgtagtttcagctttttaacactcgatattgtataaaacaaagtacactttgtacttgacaaccagtaataggagtaaaaaatatgCTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTaaaa                                                                                                                                                                               <  1:762537‑M1/72‑1 (MQ=255)
                                                                        gcgtagttctgacaattaaatgagctgaaattacgctagcaggaatcatcattcaatagatttattgatttttccgaccggcgtagtttcagctttttaacactcgatattgtataaaacaaagtacactttgtacttgacaaccagtaataggagtaaaaaatatgCTTGCCTTGCCTTGCCTTGCCTTGCctaaact                                                                                                                                                                                                                       >  1:563741‑M1/168‑199 (MQ=255)
                                                                        gcgtagttctgacaattaaatgagctgaaattacgctagcaggaatcatcattcaatagatttattgatttttccgaccggcgtagtttcagctttttaacactcgatattgtataaaacaaagtacactttgtacttgacaaccagtaataggagtaaaaaatatgCTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCggag                                                                                                                                                                   <  1:651693‑M1/84‑1 (MQ=255)
                                                                        gcgtagttctgacaattaaatgagctgaaattacgctagcaggaatcatcattcaatagatttattgatttttccgaccggcgtagtttcagctttttaacactcgatattgtataaaacaaagtacactttgtacttgacaaccagtaataggagtaaaaaatatgCTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCggag                                                                                                                                                                   <  1:515156‑M1/84‑1 (MQ=255)
                                                                        gcgtagttctgacaattaaatgagctgaaattacgcaagcaggaatcatcattcaatagatttattgatttttccgaccggcgtagtttcagctttttaacactcgatattgtataaaacaaagtacactttgtacttgacaaccagtaataggagtaaaaaatatgCTTGCCTTGCCTTGCCTTGCCTTGCctaaac                                                                                                                                                                                                                        <  1:446700‑M1/31‑1 (MQ=255)
                                                                            agttctgacaattaaatgagctgaaattacgctagcaggaatcatcattcaatagatttattgatttttccgaccggcgtagtttcagctttttaacactcgatattgtataaaacaaagtacactttgtacttgacaaccagtaataggagtaaaaaatatgCTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCAc                                                                                                                                                                                          <  1:230110‑M1/61‑1 (MQ=255)
                                                                             gttctgacaattaaatgagctgaaattacgctagcaggaatcatcattcaatagatttattgatttttccgaccggcgtagtttcagctttttaacactcgatattgtataaaacaaagtacactttgtacttgacaaccagtaataggagtaaaaaatatgCTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCt                                                                                                                                                                                               <  1:869215‑M1/56‑1 (MQ=255)
                                                                                     aattaaatgagctgaaattacgctagcaggaatcatcattcaatagatttattgatttttccgaccggcgtagtttcagctttttaacactcgatattgtataaaacaaagtacactttgtacttgacaaccagtaataggagtaaaaaatatgcttgccttgcc                                                                                                                                                                                                                                            <  1:416748‑M1/11‑1 (MQ=255)
                                                                                                  gaaattacgctagcaggaatcatcattcaatagatttattgatttttccgaccggcgtagtttcagctttttaacactcgatattgtataaaacaaagtacactttgtacttgacaaccagtaataggagtaaaaaatatgCTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCGGAGGATAAATGATCTCCGAAGAACtttt                                                                                                                                          <  1:99769‑M1/109‑1 (MQ=255)
                                                                                                            tagcaggaatcatcattcaatagatttattgatttttccgaccggcgtagtttcagctttttaacactcgatattgtataaaacaaagtacactttgtacttgacaaccagtaataggagtaaaaaatatgCTTGCCTTGCCTTGCCTTGCCTTGCctaaac                                                                                                                                                                                                                        <  1:82494‑M1/31‑1 (MQ=255)
                                                                                                                  gaatcatcattcaatagatttattgatttttccgaccggcgtagtttcagctttttaacactcgatattgtataaaacaaagtacactttgtacttgacaaccagtaataggagtaaaaaatatgCTtgccttgccttgccttgccttgcct                                                                                                                                                                                                                            <  1:125110‑M1/27‑1 (MQ=255)
                                                                                                                  gaatcatcattcaatagatttattgatttttccgaccggcgtagtttcagctttttaacactcgatattgtataaaacaaagtacactttgtacttgacaaccagtaataggagtaaaaaatatgCTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCAc                                                                                                                                                                                          >  1:509339‑M1/126‑186 (MQ=255)
                                                                                                                      catcattcaatagatttattgatttttccgaccggcgtagtttcagctttttaacactcgatattgtataaaacaaagtacactttgtacttgacaaccagtaataggagtaaaaaatatgCTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCGGAGGATAAATGATCTCCGAAGAACTTTTCTCTCGTGACGGCGTTTCCGc                                                                                                                     >  1:320664‑M1/122‑251 (MQ=255)
                                                                                                                             caatagatttattgatttttccgaccggcgtagtttcagctttttaacactcgatattgtataaaacaaagtacactttgtacttgacaaccagtaataggagtaaaaaatatgCTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCGGAGGATAAATGATCTCCGAAGAACTTTTCTCTCGTGACGGCGTTTCCGCTAAct                                                                                                                <  1:330672‑M1/135‑1 (MQ=255)
                                                                                                                                    tttattgatttttccgaccggcgtagtttcagctttttaacactcgatattgtataaaacaaagtacactttgtacttgacaaccagtaataggagtaaaaaatatgCTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCGGAGGATAAATGATCTCCGAAGAACTTTTCTCTCGTGACGGCGTTTCCGCTAACTCTGCAcaac                                                                                                       >  1:697985‑M1/108‑251 (MQ=255)
                                                                                                                                                               ttcagctttttaacactcgatattgtataaaacaaagtacactttgtacttgacaaccagtaataggagtaaaaaatatgCTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCGGAGGATAAATGATCTCCGAAGAACTTTTCTCTCGTGACGGCGTTTCCGCTAACTCTGCACAACTTACCCAACTTACCCAACTTAAAGctc                                                                            >  1:137326‑M1/81‑251 (MQ=255)
                                                                                                                                                                             actcgatattgtataaaacaaagtacactttgtacttgacaaccagtaataggagtaaaaaatatgCTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCGGAGGATAAATGATCTCCGAAGAACTTTTCTCTCGTGACGGCGTTTCCGCTAACTCTGCACAACTTACCCAACTTACCCAACTTAAAGCTCTTTTCCCGAATTGt                                                              <  1:137441‑M1/185‑1 (MQ=255)
                                                                                                                                                                              ctcgatattgtataaaacaaagtacactttgtacttgacaaccagtaataggagtaaaaaatatgCTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCGGAGGATAAATGATCTCCGAAGAACTTTTCTCTCGTGACGGCGTTTCCGCTAACTCTGCACAACTTACCCAACTTACCCAACTTAAAGc                                                                              >  1:275169‑M1/66‑234 (MQ=255)
                                                                                                                                                                                  atattgtataaaacaaagtacactttgtacttgataaccagtaataggagtaaaaaatatgCTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTCATTCGGAGGATAAATGATCTCCGAAGAACTTTTCTCTCGTGACGGCGTTTCCGCTAACTCTGCACAACTTACCCAACTTACCCAACTTAAAGc                                                                              <  1:34500‑M1/169‑1 (MQ=255)
                                                                                                                                                                                  atattgtataaaacaaagtacactttgtacttgacaaccagtaataggagtaaaaaatatgCTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCGGAGGATaaa                                                                                                                                                             <  1:42710‑M1/90‑1 (MQ=255)
                                                                                                                                                                                  atattgtataaaacaaagtacactttgtacttgacaaccagtaataggagtaaaaaatatgCTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCGGAGGATAAATGATCTCCGAAGAACTTTTCTCTCGTGACGGCGTTTCCGCTAACTCTGCACAACTTACCCAACTTACCCAACTTAAAGCTCTTTTCCCGAATTGTTTCGa                                                         <  1:22655‑M1/190‑1 (MQ=255)
                                                                                                                                                                                       gtataaaacaaagtacactttgtacttgacaaccagtaataggagtaaaaaatatgcttgccttgccttgccttgcct                                                                                                                                                                                                                                 >  1:79286‑M1/57‑78 (MQ=255)
                                                                                                                                                                                       gtataaaacaaagtacactttgtacttgacaaccagtaataggagtaaaaaatatgCTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCGGAGGATAAATGATCTCCGAAGAACTTTTCTCTCGTGACGGCGTTTCCGCTAACTCTGCACAACTTACCCAACTTACCCAACTTAAAGCTCTTTTcc                                                                      >  1:734897‑M1/57‑233 (MQ=255)
                                                                                                                                                                                       gtataaaacaaagtacactttgtacttgacaaccagtaataggagtaaaaaatatgCTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCGGAGGATAAATGATCTCCGAAGAACTTTTCTCTCGTGACGGCGTTTCCGCTAACTCTGCACAACTTACCCAACTTACCCAACTTAAAGCTCTTTTCCCGAATTGTTTCGATGTaa                                                    <  1:170108‑M1/195‑1 (MQ=255)
                                                                                                                                                                                               caaagtacactttgtacttgacaaccagtaataggagtaaaaaatatgCTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCGGAGGATAAATGATCTCCGAAGAACTTTTCTCTCGTGACGGCGTTTCCGCTAACTCTGCACAACTTACCCAACTTACCCAACTTAAAGCTCTTTTCCCGAATTGTTTCGATGTAAGCGGTCaa                                            <  1:511262‑M1/203‑1 (MQ=255)
                                                                                                                                                                                                   gtacactttgtacttgacaaccagtaataggagtaaaaaatatgCTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCGGAGGATAAATGATCTCCGAAGAACTTTTCTCTCGTGACGGCGTTTCCGCTAACTCTGCAcaacttacc                                                                                                  <  1:71705‑M1/149‑1 (MQ=255)
                                                                                                                                                                                                       actttgtacttgacaaccagtaataggagtaaaaaatatgCTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCGGAGGATAAATGATCTCCGAAGAACTTTTCTCTCGTGACGGCGTTTCCGCTAACTCTGCACAACTTACCCAACTTACCCAACTTAAAGCTCTTTTCCCGAATTGTTTCGATGTAAGCGGTCAATTTttatt                                    >  1:843693‑M1/41‑251 (MQ=255)
                                                                                                                                                                                                            gtacttgacaaccagtaataggagtaaaaaatatgCTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAA‑TTTATTCGGAGGATAAATGATCTCCGAAGAACTTTTCTCTCGTGACGGCGTTTCCGCTAACTCTGCACAACTTACCCAACTTACCCAACTTAAAGCTCTTTTCCCGAATTGTTTCGATGTAAGCGGTCAATTTtta                                      >  1:636407‑M1/36‑243 (MQ=255)
                                                                                                                                                                                                                    caaccagtaataggagtaaaaaatatgCTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCGGAGGATAAATGATCTCCGAAGAACTTTTCTCTCGTGACGGCGTTTCCGCTAACTCTGCACAACTTACCCAACTTACCCAACTTAAAGCTCTTTTCCCGAATTGTTTCGATGTAAGCGGTCAATTTTTATTAGAAAAATTCCaa                       <  1:542583‑M1/224‑1 (MQ=255)
                                                                                                                                                                                                                      accagtaataggagtaaaaaatatgCTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCGGAGGATAAATGATCTCCGAAGAACTTTTCTCTCGTGACGGCGTTTCCGCTAACTCTGCACAACTTACCCAACTTACCCAACTTAAAGCTCTTTTCCCGAATTGTTTCGATGTAAGCGGTCAATTTTTATTAGAAAAATTCCAAt                      <  1:542898‑M1/225‑1 (MQ=255)
                                                                                                                                                                                                                      accagtaataggagtaaaaaatatgCTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCGGAGGATAAATGATCTCCGAAGAACTTTTCTCTCGTGACGGCGTTTCCGCTAACTCTGCACAACTTACCCAACTTACCCAACTTAAAGCTCTTTTCCCGAATTGTTTCGATGTAAGCGGTCAATTTTTATTAGAAAAATTCCAATc                     >  1:441213‑M1/26‑251 (MQ=255)
                                                                                                                                                                                                                       ccagtaataggagtaaaaaatatgCTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCGGAGGATAAATGATCTCCGAAGAACTTTTCTCTCGTGACGGCGTTTCCGCTAACTCTGCACAACTTACCCAACTTACCCAACTTAAAGCTCTTTTCCCGAATTGTTTCGATGTAAGCGGTCAATTTttatt                                    <  1:858580‑M1/211‑1 (MQ=255)
                                                                                                                                                                                                                       ccagtaataggagtaaaaaatatgCTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCGGAGGATAAATGATCTCCGAAGAACTTTTCTCTCGTGACGGCGTTTCCGCTAACTCTGCACAACTTACCCAACTTACCCAACTTAAAGCTCTTTTCCCGAATTGTTTCGATGTAAGCGGTCAATTTTTATTAGAAAAATTCCAAt                      <  1:469339‑M1/225‑1 (MQ=255)
                                                                                                                                                                                                                        cagtaataggagtaaaaaatatgCTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCGGAGGATAAATGATCTCCGAAGAACTTTTCTCTCGTGACGGCGTTTCCGCTAACTCTGCACAACTTACCCAACTTACCCAACTTAAAGCTCTTTTCCCGAATTGTTTCGATg                                                       <  1:66360‑M1/192‑1 (MQ=255)
                                                                                                                                                                                                                          gtaataggagtaaaaaatatgCTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCGGAGGATAAATGATCTCCGAAGAACTTTTCTCTCGTGACGGCGTTTCCGCTAACTCTGCACAACttacccaacttacccaac                                                                                     <  1:145351‑M1/162‑1 (MQ=255)
                                                                                                                                                                                                                                ggagtaaaaaatatgCTTGCCTTGCCTTGCCTTGCCTTGCctaaact                                                                                                                                                                                                                       >  1:885899‑M1/16‑47 (MQ=255)
                                                                                                                                                                                                                                       aaaatatgCTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCGGAGGATAAATGATCTCCGAAGAACTTTTCTCTCGTGACGGCGTTTCCGCTAACTCTGCACAACTTACCCAACTTACCCAACTTAAAGCTCTTTTCCCGAATTGTTTCGATGTAAGCGGTCAATTTTTATTAGAAAAATTCCAATCTGAAATTGCAGAAc       >  1:870213‑M1/9‑248 (MQ=255)
                                                                                                                                                                                                                                         aatatgCTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCGGAGGATAAATGATCTCCGAAGAACTTTTCTCTCGTGACGGCGTTTCCGCTAACTCTGCACAACTTACCCAACTTACCCAACTTAAAGCTCTTTTCCCGAATTGTTTCGATGTAAGCGGTCAATTTTTATTAGAAAAATTCCAATCTGAAATTGCAGAACAAGcc  >  1:305433‑M1/7‑251 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    
CGTGGAAGAATTTTTAGACATTCTGAAAAAACAAGTGAGAAATCGCGAATTGAAATTATTGGGTGAGGCGTAGTTCTGACAATTAAATGAGCTGAAATTACGCTAGCAGGAATCATCATTCAATAGATTTATTGATTTTTCCGACCGGCGTAGTTTCAGCTTTTTAACACTCGATATTGTATAAAACAAAGTACACTTTGTACTTGACAACCAGTAATAGGAGTAAAAAATATGCTTGCCTTGCCTTGCCTTGCCTTGCCTTGCCTAAACTAAACGTTAATTGCAAGCAAAATCTATCACTTTTCTTAAAAATTTATTCGGAGGATAAATGATCTCCGAAGAACTTTTCTCTCGTGACGGCGTTTCCGCTAACTCTGCACAACTTACCCAACTTACCCAACTTAAAGCTCTTTTCCCGAATTGTTTCGATGTAAGCGGTCAATTTTTATTAGAAAAATTCCAATCTGAAATTGCAGAACAAGCC  >  1/1420556‑1421039

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑