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 48 [0] [0] 48 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)48 (0.840) 46/488 0.1 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

AACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGATTCGCGCGACGACGGCTGGCAGCGTGAAGTGAGCGACGATCGGCGCAGACAATATGACGACCGCCGCCGCCAGTTTGAAGATCGCCGCAGGCAGC  >  NC_000913/4313615‑4314111
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                        
aaCCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCGGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaa                                                                                                                                                                                                                                                        <  2:314330‑M1/251‑3 (MQ=255)
 aCCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGCTTAATGTGGATAGAGTTTCTTTTTGAGGTTGaa                                                                                                                                                                                                                                                        >  1:201898‑M1/1‑248 (MQ=255)
  ccTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGGGCTGGCTTTGGTGTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGTTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaat                                                                                                                                                                                                                                                      <  2:149939‑M1/251‑5 (MQ=255)
     tATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatcc                                                                                                                                                                                                                                                    <  1:312227‑M1/250‑7 (MQ=255)
      aTTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccat                                                                                                                                                                                                                                                  >  1:198782‑M1/1‑243 (MQ=255)
         ccGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtg                                                                                                                                                                                                                                               >  1:99425‑M1/1‑240 (MQ=255)
                  tttCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGTTAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttg                                                                                                                                                                                                                                       >  1:4652‑M1/1‑231 (MQ=255)
                   ttCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttga                                                                                                                                                                                                                                      <  2:249782‑M1/250‑21 (MQ=255)
                     cgcgGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatca                                                                                                                                                                                                                                   <  1:173726‑M1/251‑24 (MQ=255)
                                aaCGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaac                                                                                                                                                                                                                        >  2:417328‑M1/1‑217 (MQ=255)
                                  cGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacac                                                                                                                                                                                                                      >  1:139877‑M1/1‑215 (MQ=255)
                                   gACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacact                                                                                                                                                                                                                     >  2:17447‑M1/1‑214 (MQ=255)
                                     cacCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTGTTTGAGGTTGaaatccatgtgtgaagttgagcacaaatttaaacactgg                                                                                                                                                                                                                   <  2:386491‑M1/251‑40 (MQ=255)
                                                 ggCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaa                                                                                                                                                                                                         <  1:255072‑M1/249‑50 (MQ=255)
                                                 ggCTTTTGGATCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGCGGTTGaaatccatgtgggaagttgagcacaaatttaaacactggtagggtaaaa                                                                                                                                                                                                         >  2:88728‑M1/1‑200 (MQ=255)
                                                                 aGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgccca                                                                                                                                                                                        <  1:282645‑M1/251‑67 (MQ=255)
                                                                  gTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaa                                                                                                                                                                                       <  2:99425‑M1/251‑68 (MQ=255)
                                                                               aaCTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaac                                                                                                                                                                          <  2:198809‑M1/251‑81 (MQ=255)
                                                                                                  ggATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattcc                                                                                                                                                       <  1:401326‑M1/251‑100 (MQ=255)
                                                                                                     tGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattcctt                                                                                                                                                     >  1:246205‑M1/1‑149 (MQ=255)
                                                                                                                gcgcAAAGTACGAAGGCGGACGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGCGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTACTGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgggtgaagttgatcacaaatttaaacactgggagggtaaaaaggtcattaactgcccaatttaggcgtcaactggtttgattgtaaattcattaaccggagggtg                                                                                                                                         >  2:69934‑M1/1‑138 (MQ=255)
                                                                                                                           gAAGGCGGCCGTCATCAACAACGCGAGGCGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggccattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggggcaagcaaaccc                                                                                                                              >  2:438443‑M1/1‑127 (MQ=255)
                                                                                                                                 ggTCGTCATCAACATCGCGTGGAGGCGATTGCGGCAATAGAGCAGCGGAGAAATTGATATTCATCCACTTCTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcatacccgctacg                                                                                                                        <  1:320159‑M1/251‑131 (MQ=255)
                                                                                                                                    cgtcATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGAGTAATGTGGATAGAGTTTCTTTTGGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgttcattccttaaccggagggtgtaagcaaacccgctacgctt                                                                                                                     <  2:439849‑M1/251‑134 (MQ=255)
                                                                                                                                                 gcgTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTGGCATGGATGAGTAATGATTAATGGGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaatgcaggcgtcaactggtttgattgtacattccgtaaccggagggtgtaagcaaacccgctacgcttgttacagagattg                                                                                                        <  2:144698‑M1/251‑147 (MQ=255)
                                                                                                                                                           gATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagcgattgcatactgca                                                                                               >  2:84649‑M1/1‑95 (MQ=255)
                                                                                                                                                                    aaTAGAGCAGCGGAGCAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattc                                                                                           >  1:327810‑M1/1‑86 (MQ=255)
                                                                                                                                                                    aaTAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattc                                                                                           <  2:327810‑M1/245‑160 (MQ=255)
                                                                                                                                                                    aaTAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcagcctgcaattcccgctc                                                                                     <  1:345700‑M1/251‑166 (MQ=255)
                                                                                                                                                                                  gaAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccg                                                                       >  2:17580‑M1/1‑72 (MQ=255)
                                                                                                                                                                                         agaTTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCGTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctccgcttgttacagagattgcatcctgcatttcccgctccccttttgcggtcttcgagct                                                                <  2:415738‑M1/251‑187 (MQ=255)
                                                                                                                                                                                                    ctactTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaccccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctgg                                                     >  2:336811‑M1/1‑54 (MQ=255)
                                                                                                                                                                                                      actTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctgg                                                     <  1:418890‑M1/249‑198 (MQ=255)
                                                                                                                                                                                                            aTGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtttgcg                                             <  1:17447‑M1/251‑206 (MQ=255)
                                                                                                                                                                                                              ggATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaaatggtttgattgtacattccttaaccagaggctgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtttgcgg                                            >  1:200494‑M1/1‑44 (MQ=255)
                                                                                                                                                                                                                  gAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtttgcggaaatg                                       <  1:325831‑M1/251‑212 (MQ=255)
                                                                                                                                                                                                                       aTGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtttgcggaaatgggcc                                   >  2:451564‑M1/1‑35 (MQ=255)
                                                                                                                                                                                                                         gATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagatttcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtttgcggaaatgggccaa                                 <  1:298282‑M1/250‑218 (MQ=255)
                                                                                                                                                                                                                          aTTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaaggtgatcacaaatttaaacaaaggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttgtgcggccgtcgcggtgatttttctggagtttgcggcaatgggccaact                               >  2:368922‑M1/1‑32 (MQ=255)
                                                                                                                                                                                                                                  ggATAGAGTTTCTTTTTGAGGTTGaaatccatttgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtttgcggaaatgggccaactctgcgagg                       <  1:194908‑M1/251‑228 (MQ=255)
                                                                                                                                                                                                                                           ttCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaactctggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctcccctttttcggccgtcgcgcttagttttctggcgtttgcggaaatgtgccaactctgcgaggaaagcgat               <  2:376081‑M1/250‑236 (MQ=255)
                                                                                                                                                                                                                                             cTTTTTGAGGTTGaaatcaatgtgtgaagttgatcaaaaatttaaacacgggtagggtaataaggtcagtaactgcccaatttatacgtcaactggtttgattgtacattccttaaccggagggagtaagcaaacccgctacgcttgttacagcgattgcatcatgcaatcaacgctccccttttgcggccgtcgcgaggatttttctggcgtttgcggcaatgggccaactctgcgaggaaagctatgct            >  2:86970‑M1/1‑13 (MQ=255)
                                                                                                                                                                                                                                               ttttGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtttgcggaaatgggccaactctgcgaggaaagctatgctg           >  2:191272‑M1/1‑11 (MQ=255)
                                                                                                                                                                                                                                                 ttGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaacggcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtttgcggaaatgggccaactctgcgaggaaagctatgctgaa         <  1:162184‑M1/250‑242 (MQ=255)
                                                                                                                                                                                                                                                    agggtgaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtttgcggaaatgggccaactctgcgaggaaagctatgctgaaaagg     <  2:198782‑M1/247‑246 (MQ=255)
                                                                                                                                                                                                                                                      gTTGaaatccatgggtgaagttgatcacaaatttaaacactggtagggtaaaaagggcattaactgcccaattcaggcgtcaactggtttgattggacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtttgcggaaatgggccaactctgcgaggaaagctatgctgaaaaggaa   <  2:13024‑M1/251‑248 (MQ=255)
                                                                                                                                                                                                                                                       ttGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtttgcggaaatgggccaactctgcgaggaaagctatgctgaaaaggaa   >  1:109544‑M1/1‑3 (MQ=255)
                                                                                                                                                                                                                                                        tGaaatccatgtgtgaagtttatcacaaatntaaacactggttgggtaaaaaggtcattaactgcccaatgcaggcgtcaactggtttgattgtacatgccttaaccggaggttgtaagccaccccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtttgcggaaatgggccaactctgcgaggaaagctatgctgaaaaggaaa  <  2:224581‑M1/250‑249 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                        
AACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGATTCGCGCGACGACGGCTGGCAGCGTGAAGTGAGCGACGATCGGCGCAGACAATATGACGACCGCCGCCGCCAGTTTGAAGATCGCCGCAGGCAGC  >  NC_000913/4313615‑4314111

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