Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 4,313,864 Δ16,405 bp yjdPrdcB 18 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 4313864 4330268 16405 43 [0] [0] 43 yjdP–rdcB yjdP,phnP,phnO,phnN,phnM,phnL,phnK,phnJ,phnI,phnH,phnG,phnF,phnD,phnC,yjdN,yjdM,rdcA,rdcB

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 43138630 (0.000)43 (0.630) 43/492 0.3 100% intergenic (+64/‑5) rpiB/yjdP allose‑6‑phosphate isomerase/ribose‑5‑phosphate isomerase B/protein YjdP
?NC_000913 4330269 = 0 (0.000)intergenic (+31/‑233) rdcB/proP regulator of diguanylate cyclase RdcB/osmolyte:H(+) symporter ProP

GGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGATTCGCGCGACGACGGCTGGCAGCGTGAAGTGAGCGACGATCGGCGCAGACAATATGACGACCGCCGCCGCCAGTTTGAAGATCGCCGCAGGCAGC  >  NC_000913/4313639‑4314111
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                        
ggCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGTTGAGAAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgagtgaagttgatcaca                                                                                                                                                                                                                                 >  1:176322‑M1/1‑225 (MQ=255)
         aCGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaaca                                                                                                                                                                                                                       <  1:435735‑M1/251‑36 (MQ=255)
                         ggCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaag                                                                                                                                                                                                       <  2:226628‑M1/251‑52 (MQ=255)
                                   tCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaac                                                                                                                                                                                              <  1:9306‑M1/250‑61 (MQ=255)
                                                 gctCGAACGGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCTGGAGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATTCACTACTAGCATGGATGAGTAAGGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcg                                                                                                                                                                               <  2:57581‑M1/250‑76 (MQ=255)
                                                    cGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTATGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggaagggtaaaaaggtcattaactgcccaatttaggcgtc                                                                                                                                                                             >  2:220302‑M1/1‑173 (MQ=255)
                                                       aCTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaact                                                                                                                                                                         >  1:373974‑M1/1‑170 (MQ=255)
                                                          ggCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactgg                                                                                                                                                                       <  1:278572‑M1/250‑84 (MQ=255)
                                                                      tgtgGATGCGGGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCTTCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGGTGaaatccatgtgtgaagctgagcacaaatttaaacactgggagggtaaaaagggcattaactgcccgattgaggcgtcaactggtttgattgtaca                                                                                                                                                           <  2:409715‑M1/250‑96 (MQ=255)
                                                                          gATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattcct                                                                                                                                                      >  2:312385‑M1/1‑151 (MQ=255)
                                                                               gTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaacc                                                                                                                                                 >  2:437410‑M1/1‑146 (MQ=255)
                                                                                tGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccg                                                                                                                                                >  1:28455‑M1/1‑145 (MQ=255)
                                                                                 ggCTGGGCGCACAGTACCAAGTCGGTCGTACTCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacacttgtaggttaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccgg                                                                                                                                               <  2:412954‑M1/251‑108 (MQ=255)
                                                                                    tGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCCTCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGCGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagg                                                                                                                                            >  1:137645‑M1/1‑141 (MQ=255)
                                                                                    ggggcGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagg                                                                                                                                            <  2:484114‑M1/250‑111 (MQ=255)
                                                                                     gggCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggag                                                                                                                                             <  2:472663‑M1/249‑110 (MQ=255)
                                                                                      ggCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggt                                                                                                                                          <  2:28455‑M1/251‑113 (MQ=255)
                                                                                         gcACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgt                                                                                                                                        >  2:154831‑M1/1‑136 (MQ=255)
                                                                                                       cGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctac                                                                                                                         <  1:49690‑M1/251‑130 (MQ=255)
                                                                                                       cGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctac                                                                                                                         >  2:422886‑M1/1‑122 (MQ=255)
                                                                                                                 caacaaCGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaactctggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctactcttgttaca                                                                                                               <  1:308538‑M1/251‑140 (MQ=255)
                                                                                                                  aacaacGCGTGGAGGCGAGTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacag                                                                                                              <  2:238185‑M1/251‑141 (MQ=255)
                                                                                                                       cgcgTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCGCTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttagccggagggtgtaagcaaacccgctacgcttgttacagagatt                                                                                                         <  2:532753‑M1/251‑146 (MQ=255)
                                                                                                                                       cGGCAATAGAGCAGCGGAGAAATTGTGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTTaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcagtaactgcccaattcaggcgtcaactggttggattgtacattccttaaccggagggtgtaagcagacccgctacgcttgttacagagattgcatcctgcaattccc                                                                                         <  2:507010‑M1/251‑162 (MQ=255)
                                                                                                                                           aaTAGAGCAGCGGCGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacacgctacgcttgttacagagattgcatcctgcaattcccgccc                                                                                     >  2:378701‑M1/1‑86 (MQ=255)
                                                                                                                                                agcagcGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctcccctt                                                                                <  2:463162‑M1/251‑171 (MQ=255)
                                                                                                                                                    gcGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgc                                                                            <  2:293638‑M1/251‑175 (MQ=255)
                                                                                                                                                         gaAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccg                                                                       >  2:539843‑M1/1‑72 (MQ=255)
                                                                                                                                                           aaTTGAGATTCATCCACTACTTGCATTGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTTAGGTTGaaaaccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcagctggtttgattgtacatgccttaaccggagggtggaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctcctcttttgccgccgtc                                                                     <  2:31863‑M1/251‑182 (MQ=255)
                                                                                                                                                               gagaTTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgc                                                                 >  2:321867‑M1/1‑66 (MQ=255)
                                                                                                                                                                 gaTTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctg                                                               >  2:316725‑M1/1‑64 (MQ=255)
                                                                                                                                                                  aTTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctga                                                              <  1:369146‑M1/251‑189 (MQ=255)
                                                                                                                                                                     cATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatt                                                            <  1:495539‑M1/250‑191 (MQ=255)
                                                                                                                                                                      aTCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtacaaagggcattaactgcccaattcaggcgtcaactggtttgattgtacattcattaaccggaggttgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcacgctccccttttgcggccgtcgcgctgatttt                                                          >  2:59477‑M1/1‑59 (MQ=255)
                                                                                                                                                                        ccACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgattttt                                                         <  1:438002‑M1/250‑194 (MQ=255)
                                                                                                                                                                              ctTGCATGGATGAGTAATGATTAATGTGAATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcaacaatttaaacactggtagggtcaaacggtcattaactgcccaattaaggcgtcaactggtttgattgtacattccttaaccggagggtgaaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctggtttttctggcgt                                                  >  2:99513‑M1/1‑51 (MQ=255)
                                                                                                                                                                                tGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgttt                                                >  1:203581‑M1/1‑49 (MQ=255)
                                                                                                                                                                                                 aTTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtttgcggaaatgggccaact                               >  2:540650‑M1/1‑32 (MQ=255)
                                                                                                                                                                                                        tgGATAGCGTTTATTTTTGAGGCTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtttgcggaaatgggccaactctgcg                          >  2:523954‑M1/1‑25 (MQ=255)
                                                                                                                                                                                                                 tttCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtttgcggaaatgggccaactctgcgaggaaagctat               >  1:264798‑M1/1‑16 (MQ=255)
                                                                                                                                                                                                                         tGAGGTTGaaatccatgtgtgaagttgaacacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtaccttccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctcccctttttcggccgtcgcgctgatttttcttgcgtatgcggaaatgggccaactctgcgaagaaagctatgctgaaa        <  2:91395‑M1/250‑243 (MQ=255)
                                                                                                                                                                                                                             gTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtttgcggaaatgggccaactctgcgaggaaagctatgctgaaaaggaa   <  1:261527‑M1/251‑248 (MQ=255)
                                                                                                                                                                                                                               tGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaacggcccaattcaggcgtcaacgggtttgattgtacattccttaaccggagggggtaagcaacccggctacgcttgttacagagattgcatcctgcaatgcccgctccccttttgcggccgtcgcgctgatttttctggcgtttgcggaaatgggccaactctgcgaggaaagctatgctgaaaaggaaa  <  1:204252‑M1/250‑249 (MQ=255)
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                        
GGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGATTCGCGCGACGACGGCTGGCAGCGTGAAGTGAGCGACGATCGGCGCAGACAATATGACGACCGCCGCCGCCAGTTTGAAGATCGCCGCAGGCAGC  >  NC_000913/4313639‑4314111

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: