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 44 [0] [0] 44 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)44 (0.890) 42/488 0.2 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

CCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGATTCGCGCGACGACGGCTGGCAGCGTGAAGTGAGCGACGATCGGCGCAGACAATATGACGACCGCCGCCGCCAGTTTGAAGATCGCCGCAGGCA  >  NC_000913/4313617‑4314109
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
ccTTAGTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaat                                                                                                                                                                                                                                                    <  2:2706‑M1/251‑5 (MQ=255)
   tATTCCGCGTAACTTTCGCGGCAGCATAACGACACCAACGTTCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTCCTGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatc                                                                                                                                                                                                                                                   <  2:80964‑M1/249‑6 (MQ=255)
    aTTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGCGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGCGTTTCTTTTTGAGGTTGaaatccat                                                                                                                                                                                                                                                >  2:48069‑M1/1‑243 (MQ=255)
                     cgGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcac                                                                                                                                                                                                                                >  1:414850‑M1/1‑226 (MQ=255)
                            aTAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAACATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaattta                                                                                                                                                                                                                         >  1:324437‑M1/1‑219 (MQ=255)
                              aaCGACACCAACGTGCTGGCTTCTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGAGGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaac                                                                                                                                                                                                                      <  2:391123‑M1/251‑35 (MQ=255)
                               aCGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaac                                                                                                                                                                                                                      >  2:125187‑M1/1‑216 (MQ=255)
                                  acacCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactg                                                                                                                                                                                                                  <  1:48069‑M1/251‑39 (MQ=255)
                                                gCTTTTGGTTCACGAGTGGTTGGCCTCGAACTAGCAAAAATGATTGTGGACGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGCTTCATCCACTACTTGCATGGTTGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGTTTGaaatccatgtgtgaagttgatcacaaatttcaacactggtagggtaaaacgg                                                                                                                                                                                                    >  2:87059‑M1/1‑199 (MQ=255)
                                                   tttGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggt                                                                                                                                                                                                   <  1:276474‑M1/249‑54 (MQ=255)
                                                        ttCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaac                                                                                                                                                                                            >  2:280333‑M1/1‑191 (MQ=255)
                                                          cACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTATGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGTGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgacgttgatcacaaatttaaacactggtagggtaaaaaggtcattaactg                                                                                                                                                                                          >  1:158358‑M1/1‑189 (MQ=255)
                                                              aGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgccc                                                                                                                                                                                       >  1:64058‑M1/1‑185 (MQ=255)
                                                                           gAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtca                                                                                                                                                                          >  1:265191‑M1/1‑172 (MQ=255)
                                                                                ggCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggt                                                                                                                                                                    <  2:311931‑M1/251‑85 (MQ=255)
                                                                                   aaaaaTGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttg                                                                                                                                                                 <  2:133767‑M1/251‑88 (MQ=255)
                                                                                                   gCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattcctta                                                                                                                                                  >  1:267195‑M1/1‑148 (MQ=255)
                                                                                                          tGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagg                                                                                                                                          >  2:316641‑M1/1‑141 (MQ=255)
                                                                                                                cgccgtacgTAGGCGGTCGTCATCAGCAACCCGTGGAGGCGATTACGGCAATAGCGCAACTGAGAAATTGAGATTCATCCACTACTCGCACGGATGAGTAATGATTAACGTGTATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgctcacaaatttaacccctggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttacccggagggtgtaa                                                                                                                                    <  1:354772‑M1/247‑117 (MQ=255)
                                                                                                                 acaGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgta                                                                                                                                     >  2:52661‑M1/1‑134 (MQ=255)
                                                                                                                  caGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagc                                                                                                                                  <  1:311120‑M1/251‑119 (MQ=255)
                                                                                                                        gAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCAGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaaggtgatcacaaatttaaacactggtagggtaaaaagggcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaaccc                                                                                                                            >  1:151576‑M1/1‑127 (MQ=255)
                                                                                                                           ggCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgct                                                                                                                         <  1:309486‑M1/251‑128 (MQ=255)
                                                                                                                                  gtcATCAACAGCGCTTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAAGTGAGATTCATCCACTACTTGCATGGATGAGGAATGAGTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcagtcgtctactggtgtgattgtacattccttaactggagggtgtaagcaaacccgctacgcttg                                                                                                                  <  1:349394‑M1/251‑135 (MQ=255)
                                                                                                                                     atcaACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgtta                                                                                                               <  1:101475‑M1/251‑138 (MQ=255)
                                                                                                                                       caacaaCGCGTGGAGGCGATTACGGCAATAGAGCAGCTGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgtttgttac                                                                                                              <  2:151576‑M1/250‑139 (MQ=255)
                                                                                                                                           aacGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTATGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagaga                                                                                                         <  2:316023‑M1/251‑144 (MQ=255)
                                                                                                                                                                    agagCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccc                                                                                 <  1:370622‑M1/250‑168 (MQ=255)
                                                                                                                                                                          gcGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgc                                                                          <  1:430524‑M1/251‑175 (MQ=255)
                                                                                                                                                                                             tCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgaacattccttacccggagggtggaagcaaacccgctacgcttgttacagtgattgcatcctgcaattcccgctccccttttgcgtccgtcgcgctgattttt                                                       >  2:23820‑M1/1‑58 (MQ=255)
                                                                                                                                                                                               cACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttct                                                     <  1:125187‑M1/251‑196 (MQ=255)
                                                                                                                                                                                                       gCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtttg                                             <  1:301704‑M1/251‑204 (MQ=255)
                                                                                                                                                                                                             aTGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtttgcggaa                                        >  1:438966‑M1/1‑42 (MQ=255)
                                                                                                                                                                                                               gAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtttgcggaaat                                      <  2:157445‑M1/250‑211 (MQ=255)
                                                                                                                                                                                                               gAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccacttttgcggccgtcgcgctgatttttctggcgtttgcggaaatg                                     <  2:49217‑M1/251‑212 (MQ=255)
                                                                                                                                                                                                                    ctggttcatgtggatATAGTTACTTTTAGAGGTTGaaattcaagtgtgaagttgatcacaaatttaaacactggtagggcaaaaaggtcattaactgccaaattcaggcgtcaactggtttgattgtacattccttacccggagggtgtaagaaaacacgctacgattgttacagagattgcatcatgcaattcccgctccccttttccggccgtcgcgctgatttttctggcgtttgcggacatgggcc                                 >  2:369697‑M1/8‑35 (MQ=255)
                                                                                                                                                                                                                      gATTAATGTGGCTAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtttgcggaaatgggcccac                              >  2:377886‑M1/1‑33 (MQ=255)
                                                                                                                                                                                                                      gATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtttgcggaaatgggccaac                              <  1:259154‑M1/251‑219 (MQ=255)
                                                                                                                                                                                                                      gATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtttgcggaaatgggccaac                              <  2:414850‑M1/251‑219 (MQ=255)
                                                                                                                                                                                                                           aTGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtttgcggaaatgggccaactctgc                         >  1:44017‑M1/1‑28 (MQ=255)
                                                                                                                                                                                                                             gtgGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtttgcggaaatgggccaactctgcgg                       >  2:431341‑M1/1‑26 (MQ=255)
                                                                                                                                                                                                                                  tAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtttgcggaaatgggccaactctgcgaggca                   >  2:22429‑M1/1‑21 (MQ=255)
                                                                                                                                                                                                                                      gTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtttgcggaaatgggccaactctgcgaggaaagct               >  1:147724‑M1/1‑17 (MQ=255)
                                                                                                                                                                                                                                                   gTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtttgcggaaatgggccaactctgcgaggaaagctatgctgaaaagga  >  1:28678‑M1/1‑4 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
CCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGATTCGCGCGACGACGGCTGGCAGCGTGAAGTGAGCGACGATCGGCGCAGACAATATGACGACCGCCGCCGCCAGTTTGAAGATCGCCGCAGGCA  >  NC_000913/4313617‑4314109

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