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

ACCGTTCCGACGGGCGGGGAGCATGCTGATTTTTTAAGTTGGCTGGGAAAGTGGACCACGGCAGATAAACAACCCGTTTGCTGGTCAGTGACCCAACGCTGGCAAACTGTCGCGCTGGGGATGCCACGACTCTGTTCAGCGCAGCGTCTGGCGGGGGCAATGCTCGAGGAAATCTTCTCTGTAAATTTGGCGTAAATAATCAGTTACATCAATGAGTCCTAAACGAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGTTTGCGGAAATGGGCCAACTCTGCGAGGAAAGCTATGCTGAAAAGGAAA  >  NC_000913/4330044‑4330516
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       
ggcagcataacgacaccaacgtgctggcttttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggttgagaaatgattaatgtggatagagtttctttttgaggttgAAATCCATGAGTGAAGTTGATcaca                                                                                                                                                                                                                                 >  1:176322‑M2/226‑250 (MQ=255)
         acgacaccaacgtgctggcttttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAaca                                                                                                                                                                                                                       <  1:435735‑M2/35‑1 (MQ=255)
                         ggcttttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAg                                                                                                                                                                                                       <  2:226628‑M2/51‑1 (MQ=255)
                                   tcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAAc                                                                                                                                                                                              <  1:9306‑M2/60‑1 (MQ=255)
                                                 gctcgaacgggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggctggagtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcattcactactagcatggatgagtaaggattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCg                                                                                                                                                                               <  2:57581‑M2/75‑1 (MQ=255)
                                                    cgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtatgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGAAGGGTAAAAAGGTCATTAACTGCCCAATTTAGGCGTc                                                                                                                                                                             >  2:220302‑M2/174‑250 (MQ=255)
                                                       actggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACt                                                                                                                                                                         >  1:373974‑M2/171‑251 (MQ=255)
                                                          ggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTgg                                                                                                                                                                       <  1:278572‑M2/83‑1 (MQ=255)
                                                                      tgtggatgcggggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcttccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgagggtgAAATCCATGTGTGAAGCTGAGCACAAATTTAAACACTGGGAGGGTAAAAAGGGCATTAACTGCCCGATTGAGGCGTCAACTGGTTTGATTGTACa                                                                                                                                                           <  2:409715‑M2/95‑1 (MQ=255)
                                                                          gatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCt                                                                                                                                                      >  2:312385‑M2/152‑251 (MQ=255)
                                                                               gtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAAcc                                                                                                                                                 >  2:437410‑M2/147‑251 (MQ=255)
                                                                                tggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCg                                                                                                                                                >  1:28455‑M2/146‑251 (MQ=255)
                                                                                 ggctgggcgcacagtaccaagtcggtcgtactcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTTGTAGGTTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCgg                                                                                                                                               <  2:412954‑M2/107‑1 (MQ=255)
                                                                                    tgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcctccactacttgcatggatgagtaatgattaatgtggatagcgtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAgg                                                                                                                                            >  1:137645‑M2/142‑251 (MQ=255)
                                                                                    ggggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAgg                                                                                                                                            <  2:484114‑M2/110‑1 (MQ=255)
                                                                                     gggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAg                                                                                                                                             <  2:472663‑M2/109‑1 (MQ=255)
                                                                                      ggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGgt                                                                                                                                          <  2:28455‑M2/112‑1 (MQ=255)
                                                                                         gcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGgtgt                                                                                                                                        >  2:154831‑M2/137‑250 (MQ=255)
                                                                                                       cggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTAc                                                                                                                         <  1:49690‑M2/129‑1 (MQ=255)
                                                                                                       cggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTAc                                                                                                                         >  2:422886‑M2/123‑251 (MQ=255)
                                                                                                                 caacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACTCTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACTCTTGTTACa                                                                                                               <  1:308538‑M2/139‑1 (MQ=255)
                                                                                                                  aacaacgcgtggaggcgagtacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACag                                                                                                              <  2:238185‑M2/140‑1 (MQ=255)
                                                                                                                       cgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccgctacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAGCCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGAtt                                                                                                         <  2:532753‑M2/145‑1 (MQ=255)
                                                                                                                                       cggcaatagagcagcggagaaattgtgattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggtttAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCAGTAACTGCCCAATTCAGGCGTCAACTGGTTGGATTGTACATTCCTTAACCGGAGGGTGTAAGCAGACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTccc                                                                                         <  2:507010‑M2/161‑1 (MQ=255)
                                                                                                                                           aatagagcagcggcgaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACACGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCcgccc                                                                                     >  2:378701‑M2/87‑249 (MQ=255)
                                                                                                                                                agcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCtt                                                                                <  2:463162‑M2/170‑1 (MQ=255)
                                                                                                                                                    gcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGc                                                                            <  2:293638‑M2/174‑1 (MQ=255)
                                                                                                                                                         gaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCg                                                                       >  2:539843‑M2/73‑251 (MQ=255)
                                                                                                                                                           aattgagattcatccactacttgcattgatgagtaatgattaatgtggatagagtttcttttttaggttgAAAACCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAGCTGGTTTGATTGTACATGCCTTAACCGGAGGGTGGAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCTCTTTTGCCGCCGTc                                                                     <  2:31863‑M2/181‑1 (MQ=255)
                                                                                                                                                               gagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCgcgc                                                                 >  2:321867‑M2/67‑251 (MQ=255)
                                                                                                                                                                 gattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTg                                                               >  2:316725‑M2/65‑251 (MQ=255)
                                                                                                                                                                  attcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGa                                                              <  1:369146‑M2/188‑1 (MQ=255)
                                                                                                                                                                     catccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGAtt                                                            <  1:495539‑M2/190‑1 (MQ=255)
                                                                                                                                                                      atccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTACAAAGGGCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCATTAACCGGAGGTTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCACGCTCCCCTTTTGCGGCCGTCGCGCTGAtttt                                                          >  2:59477‑M2/60‑251 (MQ=255)
                                                                                                                                                                        ccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGAttttt                                                         <  1:438002‑M2/193‑1 (MQ=255)
                                                                                                                                                                              cttgcatggatgagtaatgattaatgtgaatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCAACAATTTAAACACTGGTAGGGTCAAACGGTCATTAACTGCCCAATTAAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGAAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGGTTTTTCTGGCGt                                                  >  2:99513‑M2/52‑251 (MQ=255)
                                                                                                                                                                                tgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGttt                                                >  1:203581‑M2/50‑251 (MQ=255)
                                                                                                                                                                                                 attaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGTTTGCGGAAATGGGCCAAct                               >  2:540650‑M2/33‑251 (MQ=255)
                                                                                                                                                                                                        tggatagcgtttatttttgaggctgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGTTTGCGGAAATGGGCCAACTCTGCg                          >  2:523954‑M2/26‑249 (MQ=255)
                                                                                                                                                                                                                 tttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGTTTGCGGAAATGGGCCAACTCTGCGAGGAAAGCTAt               >  1:264798‑M2/17‑251 (MQ=255)
                                                                                                                                                                                                                         tgaggttgAAATCCATGTGTGAAGTTGAACACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACCTTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTTCGGCCGTCGCGCTGATTTTTCTTGCGTATGCGGAAATGGGCCAACTCTGCGAAGAAAGCTATGCTGaaa        <  2:91395‑M2/242‑1 (MQ=255)
                                                                                                                                                                                                                             gttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGTTTGCGGAAATGGGCCAACTCTGCGAGGAAAGCTATGCTGAAAAGGaa   <  1:261527‑M2/247‑1 (MQ=255)
                                                                                                                                                                                                                               tgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACGGCCCAATTCAGGCGTCAACGGGTTTGATTGTACATTCCTTAACCGGAGGGGGTAAGCAACCCGGCTACGCTTGTTACAGAGATTGCATCCTGCAATGCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGTTTGCGGAAATGGGCCAACTCTGCGAGGAAAGCTATGCTGAAAAGGaaa  <  1:204252‑M2/248‑1 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       
ACCGTTCCGACGGGCGGGGAGCATGCTGATTTTTTAAGTTGGCTGGGAAAGTGGACCACGGCAGATAAACAACCCGTTTGCTGGTCAGTGACCCAACGCTGGCAAACTGTCGCGCTGGGGATGCCACGACTCTGTTCAGCGCAGCGTCTGGCGGGGGCAATGCTCGAGGAAATCTTCTCTGTAAATTTGGCGTAAATAATCAGTTACATCAATGAGTCCTAAACGAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGTTTGCGGAAATGGGCCAACTCTGCGAGGAAAGCTATGCTGAAAAGGAAA  >  NC_000913/4330044‑4330516

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: