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 65 [0] [0] 64 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)64 (1.090) 61/492 0.0 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

CGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGATTCGCGCGACGACGGCTGGCAGCGTGAAGTGAGCGACGATCGGCGCAGACAATATGACGACCGCCGCCGCCAGTTTGAAGATCGCCGCAGGCAGC  >  NC_000913/4313613‑4314111
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                        
cGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTg                                                                                                                                                                                                                                                          <  2:438236/251‑1 (MQ=255)
  aaCCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaa                                                                                                                                                                                                                                                        <  1:60708‑M1/251‑3 (MQ=255)
   aCCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaa                                                                                                                                                                                                                                                        >  1:440632‑M1/1‑248 (MQ=255)
        aTTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccat                                                                                                                                                                                                                                                  >  1:119663‑M1/1‑243 (MQ=255)
          tCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAGTGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatg                                                                                                                                                                                                                                                 <  1:278334‑M1/250‑10 (MQ=255)
            cgcgCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgt                                                                                                                                                                                                                                              <  2:37143‑M1/251‑13 (MQ=255)
              cgcAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtg                                                                                                                                                                                                                                             >  1:31703‑M1/1‑237 (MQ=255)
              cgcAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTGTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtga                                                                                                                                                                                                                                            <  1:426994‑M1/251‑15 (MQ=255)
                 aaCTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagt                                                                                                                                                                                                                                         >  1:87653‑M1/1‑234 (MQ=255)
                     ttCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttg                                                                                                                                                                                                                                       <  2:188245‑M1/249‑20 (MQ=255)
                      tCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatc                                                                                                                                                                                                                                    >  1:13911‑M1/1‑229 (MQ=255)
                      tCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgat                                                                                                                                                                                                                                     <  1:263399‑M1/250‑22 (MQ=255)
                       cgcgGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatc                                                                                                                                                                                                                                    <  2:315995‑M1/250‑23 (MQ=255)
                           gcagcaTAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaa                                                                                                                                                                                                                               <  2:268503‑M1/251‑28 (MQ=255)
                             agcaTAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatt                                                                                                                                                                                                                             >  1:362506‑M1/1‑222 (MQ=255)
                                aTAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaattta                                                                                                                                                                                                                           >  2:51480‑M1/1‑219 (MQ=255)
                                              gTGCTGGCTTTTGGTTCAAGAGTGGTTGGCCCAGAACTGGCAAAAATGATTGGGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGTTATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagtggatcacaaatttaaacactggtagggta                                                                                                                                                                                                            <  1:391287‑M1/251‑47 (MQ=255)
                                                   ggCTTTTGGTTCACGAGTGGTTGGCCGCGAACTGGGAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGGACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacacgggtagggtaaaaa                                                                                                                                                                                                        <  2:336940‑M1/250‑51 (MQ=255)
                                                       tttGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtc                                                                                                                                                                                                    <  2:432661‑M1/250‑55 (MQ=255)
                                                        ttGGTGCACTAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtca                                                                                                                                                                                                   <  1:160207‑M1/250‑56 (MQ=255)
                                                         tGGTTCACGAGTGGTTGGCCTCAAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCTATTACGGCAATAGATCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcat                                                                                                                                                                                                  <  1:111479‑M1/250‑57 (MQ=255)
                                                                 gAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgccc                                                                                                                                                                                         <  2:216309‑M1/251‑66 (MQ=255)
                                                                      gTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattc                                                                                                                                                                                    <  2:279178‑M1/251‑71 (MQ=255)
                                                                            cTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGGACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTCCTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgagaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgt                                                                                                                                                                              >  2:397866‑M1/1‑175 (MQ=255)
                                                                                 aCTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaact                                                                                                                                                                         <  2:440632‑M1/251‑82 (MQ=255)
                                                                                    ggCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactg                                                                                                                                                                        <  1:169368‑M1/249‑83 (MQ=255)
                                                                                      cAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggttt                                                                                                                                                                    <  1:101070‑M1/251‑87 (MQ=255)
                                                                                      cAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggttt                                                                                                                                                                    <  1:196000‑M1/251‑87 (MQ=255)
                                                                                          aaTGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaatgcaggcgtcaactggtttgatt                                                                                                                                                                <  2:111464‑M1/251‑91 (MQ=255)
                                                                                              aTTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtac                                                                                                                                                            <  2:117923‑M1/251‑95 (MQ=255)
                                                                                                    gATGCGGGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattcct                                                                                                                                                      <  1:326651‑M1/251‑101 (MQ=255)
                                                                                                       gCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaa                                                                                                                                                   <  2:129892‑M1/251‑104 (MQ=255)
                                                                                                         gTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaacc                                                                                                                                                 >  1:209899‑M1/1‑146 (MQ=255)
                                                                                                           ggCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccgg                                                                                                                                               >  1:242276‑M1/1‑144 (MQ=255)
                                                                                                              tGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagg                                                                                                                                            >  1:454631‑M1/1‑141 (MQ=255)
                                                                                                                    cacaGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgta                                                                                                                                       >  1:24089‑M1/1‑135 (MQ=255)
                                                                                                                      caGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGAGTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagtggatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattttacattccttaacctgagggtgtaagc                                                                                                                                    <  2:81027‑M1/251‑119 (MQ=255)
                                                                                                                        gTACGTAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttacccggagggtgtaagca                                                                                                                                   >  2:179860‑M1/1‑131 (MQ=255)
                                                                                                                              aGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgc                                                                                                                            <  2:42044‑M1/251‑127 (MQ=255)
                                                                                                                               ggCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttatccggagggtgtaagcaaacccgct                                                                                                                           <  2:119663‑M1/251‑128 (MQ=255)
                                                                                                                                                       gAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcc                                                                                                   >  1:41500‑M1/1‑100 (MQ=255)
                                                                                                                                                            gATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgca                                                                                               >  1:67504‑M1/1‑95 (MQ=255)
                                                                                                                                                              ttACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaac                                                                                                                                                                          <  2:438124‑M1/173‑81 (MQ=255)
                                                                                                                                                              ttACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaac                                                                                                                                                                          >  1:438124‑M1/1‑93 (MQ=255)
                                                                                                                                                                 cGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattccc                                                                                         <  2:87653‑M1/251‑162 (MQ=255)
                                                                                                                                                                  ggCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccg                                                                                        >  2:449561‑M1/1‑89 (MQ=255)
                                                                                                                                                                      aTAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctacaattcccgctcc                                                                                    <  1:331925‑M1/251‑167 (MQ=255)
                                                                                                                                                                        agagCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctcccc                                                                                  >  1:357115‑M1/1‑83 (MQ=255)
                                                                                                                                                                        agagCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctcccc                                                                                  >  1:399190‑M1/1‑83 (MQ=255)
                                                                                                                                                                                          agaTTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgct                                                                <  2:95342‑M1/251‑187 (MQ=255)
                                                                                                                                                                                                       actTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcg                                                   <  2:41500‑M1/251‑200 (MQ=255)
                                                                                                                                                                                                          tGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcgcccgtcgcgctgatttttctggcgttt                                                >  1:106247‑M1/1‑49 (MQ=255)
                                                                                                                                                                                                                  tGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggattgattgtacattccttaaccggtgggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtttgcggaaat                                        >  1:309419‑M1/1‑41 (MQ=255)
                                                                                                                                                                                                                   gAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtttgcggaaat                                        <  2:209899‑M1/250‑211 (MQ=255)
                                                                                                                                                                                                                                          gTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtttgcggaaat                                        <  2:31857‑M1/227‑211 (MQ=255)
                                                                                                                                                                                                                                          gTTTCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtttgcggaaat                                        >  1:31857‑M1/1‑17 (MQ=255)
                                                                                                                                                                                                                                             tCTTTTTGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtttgcggaaatgggccaactctgcgaggaaagctatg              <  1:274921‑M1/250‑237 (MQ=255)
                                                                                                                                                                                                                                                ttttGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtttgcggaaatgggccaactctgcgaggaaagctatgc             <  2:362506‑M1/248‑238 (MQ=255)
                                                                                                                                                                                                                                                  ttGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtttgcggaaatgggccaactctgcgaggaaagctatgctgaa         <  1:98136‑M1/250‑242 (MQ=255)
                                                                                                                                                                                                                                                  ttGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctcccc                                                                                  >  2:35509‑M1/1‑9 (MQ=255)
                                                                                                                                                                                                                                                  ttGAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctcccc                                                                                  <  1:35509‑M1/177‑169 (MQ=255)
                                                                                                                                                                                                                                                    gAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttaccg                                                                                                              >  2:142070‑M1/1‑7 (MQ=255)
                                                                                                                                                                                                                                                    gAGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacag                                                                                                              <  1:142070‑M1/147‑141 (MQ=255)
                                                                                                                                                                                                                                                     aGGTTGaaatccatgtgtgaagttgatcacaaatttaaacactggtagggtaaaaaggtcattaactgcccaattcaggcgtcaactggtttgattgtacattccttaaccggagggtgtaagcaaacccgctacgcttgttacagagattgcatcctgcaattcccgctccccttttgcggccgtcgcgctgatttttctggcgtttgcggaaatgggccaactctgcgaggaaagctatgctgaaaagg     >  2:144794‑M1/1‑6 (MQ=255)
                                                                                                                                                                                                                                                         tGaaatccatgtgggaagttgatcacaaatttaaacgctggtagggtaaaaaggtcattaactgcccaattcaggcgtcagctggtttgattgtgcattccttaaccggagggtgtaaggaaacccgctacgcttgttgcagagattgcatcctgcaattcccgctcccctttggcggcgggcgcggtgatttttctggcgtttgcggaaatgggcgaactctgggaggaaagctatgctgaaaaggaaa  <  1:66850‑M1/250‑249 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                        
CGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGATTCGCGCGACGACGGCTGGCAGCGTGAAGTGAGCGACGATCGGCGCAGACAATATGACGACCGCCGCCGCCAGTTTGAAGATCGCCGCAGGCAGC  >  NC_000913/4313613‑4314111

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