Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 2,066,152 Δ12,778 bp insH6[yoeF] insH6, yoeA, insD3, insC3, yoeA, isrC, flu, yeeR, yeeS, yeeT, cbeA, cbtA, yeeW, [yoeF]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 2066152 2078929 12778 85 [0] [0] 85 insH6–[yoeF] insH6,yoeA,insD3,insC3,yoeA,isrC,flu,yeeR,yeeS,yeeT,cbeA,cbtA,yeeW,[yoeF]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 20661510 (0.000)85 (0.900) 75/486 0.2 100% intergenic (‑387/+154) cobU/insH6 cobinamide‑P guanylyltransferase/cobinamide kinase/IS5 family transposase and trans‑activator
?NC_000913 2078930 = 0 (0.000)coding (356/357 nt) yoeF putative uncharacterized protein YoeF

AAATAAGATGGCATTACGGAATAGCCAGACTTTCTCGTGTTAACGTCCGACAAAATGCGTCATTGACAGCGGTTTCTGTGGATTGACACATTTCATTCAGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAG  >  NC_000913/2065903‑2066397
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
aaaTAAGATGGCATTACGGAATAGCCAGACTTTCTCGTGTTAACGTCCGACAAAATGCGTCATTGACAGCGGTTTCTGTGGATTGACACATTTCATTCAGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATa                                                                                                                                                                                                                                                       >  2:35007‑M1/1‑249 (MQ=255)
      gattgcattAAGGAATAGCCAGACTTTCTCGTGTTAACGTCCGAAAAAATGCGTCATTGAAAGCGGTTTCTGTGGATTGACACATTTCATTCAGATAATGAATTAATGCTACTGCCGGAACAACCAGCAACCGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGACTGGAGGCGACTATaatttcg                                                                                                                                                                                                                                                 >  2:233953‑M1/5‑243 (MQ=255)
          gCATTACGGAATAGCCAGACTTTCTCGTGTTAACGTCCGACAAAATGCGTCATTGACAGCGGTTTCTGTGGATTGACACATTTCATTCAGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattgg                                                                                                                                                                                                                                            <  2:797508‑M1/251‑13 (MQ=255)
            aTTACGGACTAGCCAGACTTTCTCGTGTTAACGTCCGACAAAATGCGTCATTGACAGCGGTTTCTGTGGATTGACACATTTCATTCAGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggtt                                                                                                                                                                                                                                          >  1:56412‑M1/1‑237 (MQ=255)
                 ggAATAGCCAGCCTTTCTCGTGTTAACGTCCGACAAAATGCGTCATTGACAGCGGTTTCTGTGGATTGACACATTTCATTCAGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcgg                                                                                                                                                                                                                                     >  1:152024‑M1/1‑232 (MQ=255)
                 ggAATAGCCAGACTTTCTCGTGTTAACGTCCGACAAAATGCGTCATTGACAGCGGTTTCTGTGGATTGACACATTTCATTCAGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcgg                                                                                                                                                                                                                                     >  1:689545‑M1/1‑232 (MQ=255)
                   aaTAGCCAGACTTTCTCGTGTTAACGTCCGACAAAATGCGTCATTGACAGCGGTTTCTGTGGATTGACACATTTCATTCAGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggat                                                                                                                                                                                                                                   <  1:413641‑M1/251‑22 (MQ=255)
                   aaTAGCCAGACTTTCTCGTGTTAACGTCCGACAAAATGCGTCATTGACAGCGGTTTCTGTGGATTGACACATTTCATTCAGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttc                                                                                                                                                                                                                                                  >  1:252532‑M1/1‑230 (MQ=255)
                   aaTAGCCAGACTTTCTCGTGTTAACGTCCGACAAAATGCGTCATTGACAGCGGTTTCTGTGGATTGACACATTTCATTCAGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttc                                                                                                                                                                                                                                                  <  2:252532‑M1/236‑7 (MQ=255)
                         cAGACTTTCTCGTGTTACCGTCCGACAAAATGCGTCATTGACAGCGGTTTCTGTGGATTGACACATTTCATTCAGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGATTATaatttcgattggttctcgtatgaggat                                                                                                                                                                                                                             >  2:385877‑M1/1‑224 (MQ=255)
                         cAGACTTTCTCGTGTTAACGTCCGACAAAATGCGTCATTGACAGCGGTTTCTGTGGATTGACACATTTCATTCAGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAAGTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggat                                                                                                                                                                                                                             <  1:225195‑M1/251‑28 (MQ=255)
                                         aaCGTCCGACAAAATGCGTCATTGACAGCGGTTTCTGTGGATTGACACATTTCATTCAGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacac                                                                                                                                                                                                             <  1:496133‑M1/251‑44 (MQ=255)
                                         aaCGTCCGACAAAATGCGTCATTGACAGCGGTTTCTGTGGATTGACACATTTCATTCAGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacac                                                                                                                                                                                                             >  2:125392‑M1/1‑208 (MQ=255)
                                         aaCGTCCGACAAAATGCGTCATTGACAGCGGTTTCTGTGGATTGACACATTTCATTCAGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgataca                                                                                                                                                                                                              >  2:442270‑M1/1‑208 (MQ=255)
                                         aaCGTCCGACAAAATGCGTCATTGACAGCGGTTTCTGTGGATTGACACAATTCATTCAGATAATGAATTAATGCTACCGCCGGAACAACCAGCAAACGGGCATTGTGTTCTTAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAGTTATCCGGAATGGTGGCGGCTATaatttcgattggttcttggaggaggaggatgacgttgatacac                                                                                                                                                                                                             <  2:613496‑M1/251‑44 (MQ=255)
                                                gACAAAATGCGTCATTGACAGCGGTTTCTGTGGATTGACACATTTCATTCAGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataag                                                                                                                                                                                                      <  2:777304‑M1/251‑51 (MQ=255)
                                                   aaaaTGCGTCATTGACAGCGGTTTCTGTGGATTGACACATTTCATTCAGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaat                                                                                                                                                                                                   >  2:25429‑M1/1‑198 (MQ=255)
                                                    aaaTGCGTCATTGACAGCGGTTTCTGTGGATTGACACATTTCATTCAGATAATGAATTAATGCTACTGCAGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGAGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatc                                                                                                                                                                                                  >  2:67739‑M1/1‑197 (MQ=255)
                                                             aTTGACAGCGGTTTCTGTGGATTGACACATTTCATTCAGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagca                                                                                                                                                                                          >  1:419157‑M1/1‑188 (MQ=255)
                                                               tGACAGCGGTTTCTGTGGATTGACACATTTCATTCAGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaat                                                                                                                                                                                        <  1:531209‑M1/250‑65 (MQ=255)
                                                                  cAGCGGTTTCTGTGGTTTGACACATTTCATTCAGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggctgaggttgttgacggtgatacaccgataagaatcaatgagcaatgatg                                                                                                                                                                                    >  2:440922‑M1/1‑183 (MQ=255)
                                                                     cGGTTTCTGTGGATTGACACATTTCATTCAGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttcgcggatgaggatgatgacggtgatacacagataagaagcaatgagcgatgagggtt                                                                                                                                                                                 <  2:643824‑M1/251‑72 (MQ=255)
                                                                          tCTGTGGATTGACCCATTTCAGTCAGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGTGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAAGTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaattagcaatgagggttgccg                                                                                                                                                                             <  2:403129‑M1/250‑76 (MQ=255)
                                                                               ggATTGACACATTTCATTCAGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaac                                                                                                                                                                       <  1:329142‑M1/251‑82 (MQ=255)
                                                                               ggATTGACACATTTCATTCAGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGAATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaac                                                                                                                                                                       >  1:225332‑M1/1‑170 (MQ=255)
                                                                                    gACACATTTCATTCAGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctca                                                                                                                                                                  <  2:438134‑M1/251‑87 (MQ=255)
                                                                                    gACACATTTCATTCAGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctc                                                                                                                                                                   >  1:678347‑M1/1‑165 (MQ=255)
                                                                                      cacaTTTCATTCAGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgacgggcaaccctcatt                                                                                                                                                                <  1:125392‑M1/251‑89 (MQ=255)
                                                                                        caTTTCATTCAGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattg                                                                                                                                                               >  2:719268‑M1/1‑161 (MQ=255)
                                                                                             cattcaGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataa                                                                                                                                                          <  1:719268‑M1/250‑95 (MQ=255)
                                                                                              attcaGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaa                                                                                                                                                         >  1:269131‑M1/1‑155 (MQ=255)
                                                                                                    aTAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacggga                                                                                                                                                   >  1:402990‑M1/1‑149 (MQ=255)
                                                                                                         gAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTATTCTGACATCCAGGCTATTGATTCAACGACAGCGACAGATCTGCGATGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctaggatgaggatgatgacggtgatacacatatacgaatcaatgagcaatgagggttgcagggcaaccctcattgaataaaaaggtaataact                                                                                                                                             >  2:82312‑M1/1‑144 (MQ=255)
                                                                                                          aaTTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTTAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattactt                                                                                                                                            <  2:360948‑M1/251‑109 (MQ=255)
                                                                                                          aaTTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagc                                                                                                                                                                                           >  1:610988‑M1/1‑143 (MQ=255)
                                                                                                          aaTTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagc                                                                                                                                                                                           <  2:610988‑M1/204‑62 (MQ=255)
                                                                                                              aaTGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTCTCCGGAATGGTGGCGAATATaatttagattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgcataaaacgggcattacttcgct                                                                                                                                        >  2:302076‑M1/1‑139 (MQ=255)
                                                                                                              aaTGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgct                                                                                                                                        <  2:716619‑M1/251‑113 (MQ=255)
                                                                                                              aaTGCTACTGACGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgct                                                                                                                                        >  2:625870‑M1/1‑139 (MQ=255)
                                                                                                                 gCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcg                                                                                                                                     <  2:842007‑M1/251‑116 (MQ=255)
                                                                                                                      tGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtt                                                                                                                                >  1:139581‑M1/1‑131 (MQ=255)
                                                                                                                        ccGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttt                                                                                                                              <  1:124565‑M1/251‑123 (MQ=255)
                                                                                                                         cGGAACAACCAGCAAACTGTCAGTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCTGCGATAGATCTGAGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGTTGGCGACTATaaagtcgattggttctcggatgaggatgaggacggtgattcacagataagaaacaatgagcaatgagggttgccgggcatcccacattgaataaaacgggatttacttcgtttcgccgttttt                                                                                                                             <  2:588103‑M1/251‑124 (MQ=255)
                                                                                                                         cGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgttttt                                                                                                                             >  1:746037‑M1/1‑128 (MQ=255)
                                                                                                                             acaacCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatt                                                                                                                         <  1:181356‑M1/251‑128 (MQ=255)
                                                                                                                               aacCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttc                                                                                                                       >  2:488863‑M1/1‑122 (MQ=255)
                                                                                                                               aacCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttc                                                                                                                       <  1:614942‑M1/251‑130 (MQ=255)
                                                                                                                               aacCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctc                                                                                                                                                                   >  1:313120‑M1/1‑122 (MQ=255)
                                                                                                                                      aaaCGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccat                                                                                                                 >  1:508147‑M1/1‑115 (MQ=255)
                                                                                                                                                tgtgTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagct                                                                                                      >  1:618763‑M1/1‑105 (MQ=255)
                                                                                                                                                    ttCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagcttttt                                                                                                  >  1:560784‑M1/1‑101 (MQ=255)
                                                                                                                                                        gAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttgg                                                                                                >  1:189691‑M1/1‑97 (MQ=255)
                                                                                                                                                          aaTCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagat                                                                                            <  2:783530‑M1/251‑157 (MQ=255)
                                                                                                                                                             ccAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagatatc                                                                                         <  1:264696‑M1/251‑160 (MQ=255)
                                                                                                                                                             ccAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagatatc                                                                                         <  1:18583‑M1/251‑160 (MQ=255)
                                                                                                                                                               aGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagatatcgc                                                                                       <  1:488863‑M1/251‑162 (MQ=255)
                                                                                                                                                                ggCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagatatcgcg                                                                                      <  1:35007‑M1/251‑163 (MQ=255)
                                                                                                                                                                  cTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagatatcgcggc                                                                                    >  2:80738‑M1/1‑87 (MQ=255)
                                                                                                                                                                    aTTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagatatcgcggcgt                                                                                  <  2:139581‑M1/251‑167 (MQ=255)
                                                                                                                                                                    aTTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagatatcgcggcgt                                                                                  >  2:239135‑M1/1‑85 (MQ=255)
                                                                                                                                                                                      aCAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagatatcgcggcgttctttggagagctcggca                                                                <  2:419157‑M1/251‑185 (MQ=255)
                                                                                                                                                                                      aCAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagatatcgcggcgttctttggagagctcggc                                                                 >  2:374636‑M1/1‑67 (MQ=255)
                                                                                                                                                                                        aGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagatatcgcggcgttctttggagagctcggcatt                                                              >  2:778882‑M1/1‑65 (MQ=255)
                                                                                                                                                                                            cTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGACTTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagatatcgcggcgttctttggagagctcggcatttttg                                                          >  2:78880‑M1/1‑61 (MQ=255)
                                                                                                                                                                                              gcgcTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagatatcgcggcgttctttggagagctcggcattttt                                                           <  1:632648‑M1/248‑190 (MQ=255)
                                                                                                                                                                                                    tAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagatatcgcggcgttctttggagagctcggcatttttgatgatgt                                                   >  1:757132‑M1/1‑53 (MQ=255)
                                                                                                                                                                                                             ttGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagatatcgcggcgttctttggcgagctcggcatttttgatgatgtaatcgtcaac                                         >  2:127564‑M1/1‑44 (MQ=255)
                                                                                                                                                                                                              tGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcgtcgccgtttttcatttcgcccatagctttcagctttttggagatatcgcggcgttctttggagagctcggcatttttgatgatgtaatcgtcaac                                         <  2:757132‑M1/250‑208 (MQ=255)
                                                                                                                                                                                                                 aaaaCAACGTAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagatatcgcggcgtgctttggagagctcggcatttttgatgatgcaatcgtcaacgcg                                      >  1:63999‑M1/1‑40 (MQ=255)
                                                                                                                                                                                                                     aaacgGAATAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagatatcgcggcgttctttggagagctcggcattttttatgatgtaatcgtcaacgcgatctt                                 <  2:639347‑M1/250‑216 (MQ=255)
                                                                                                                                                                                                                           aaTAGAATTATCCTGGATGGTGGCGACTATaatttcgattggttctcggatgaggatgatggcggtgatacacagataagaatcaatgaacaatgagggttgccgggcaaccctcattgaataaaacggtaattactgctcttcgctgtttttcatttcgcccatagctttcagctttttggagatatcgcggctttctatggagagctcggcttctttgatgatgtaatcgtcaacgcgatcttcattgt                           <  2:107352‑M1/251‑222 (MQ=255)
                                                                                                                                                                                                                           aaTAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagatatcgcggcgttctttggagagctcggcatttttgatgatgtaatcgtcaacgcgatcttcatagt                           >  2:365927‑M1/1‑30 (MQ=255)
                                                                                                                                                                                                                            aTAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagatatcgcggcgttctttggcgagctcggcatttttgcagctgtagtcgtcaacgcgatcttcatagtc                          >  2:768032‑M1/1‑29 (MQ=255)
                                                                                                                                                                                                                            aTAGAATTATCCGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagatatcgcggcgttctttggagagctcggcatttttgatgatgtaatcgtcaacgcgatcttcatagtc                          >  1:362469‑M1/1‑29 (MQ=255)
                                                                                                                                                                                                                                    aTCCGGAATGGTGGCGACTATaagttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagatatcgcggcgttctttggagagctcggcatttttgatgatgtaatcgtcaacgcgatcttcatagtcacctttc                   <  2:189691‑M1/250‑230 (MQ=255)
                                                                                                                                                                                                                                      ccgtaatggtGGCGACTATaatttctattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagatatcgcggcgttctttggagagctctgcatttttgatgatgtaatcgtcaacgcgatcttcagagtcaccttgcat                 <  1:301511‑M1/246‑232 (MQ=255)
                                                                                                                                                                                                                                      ccGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagatatcgcggcgttctttggagagctcggcatttttgatgatgtaatcgtcaacgcgatcttcatagtcacctttcatg                <  1:374636‑M1/251‑233 (MQ=255)
                                                                                                                                                                                                                                       cGGAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagatatcgcggcgttctttggagagctcggcatttttgatgatgtaatcgtcaacgcgatcttcatagtcacctttcatgc               <  2:56412‑M1/251‑234 (MQ=255)
                                                                                                                                                                                                                                        ggAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagatatcgcggcgttctttggagagctcggcatttttgatgatgtactcgtcaacgcgatcttcatagtcacctttcatgct              >  2:691313‑M1/1‑17 (MQ=255)
                                                                                                                                                                                                                                        ggAATGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagatatcgcggcgttctttggagagctcggcatttttgatgatgtaatcgtcaacgcgatcttcatagtcacctttcatgct              <  1:80738‑M1/251‑235 (MQ=255)
                                                                                                                                                                                                                                           aTGGTGGCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagatatcgcggcgttctttggagagctcggcatttttgatgatgtaatcgtccacgcgatcttcatagtcacctttcatgctggc           >  2:474334‑M1/1‑14 (MQ=255)
                                                                                                                                                                                                                                                ggCGACTATaatttcgcttggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagatatcgcggcgttctttggagagctcggcatttttgatgatgtaatcgtcaacgcgatcttcatagtcacctttcatgctggcgatg       >  2:840920‑M1/1‑9 (MQ=255)
                                                                                                                                                                                                                                                ggCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagatatcgcggcgttctttggagagctcggcatttttgatgatgtaatcgtcaacgcgatcttcatagtcacctttcatgctggcgatga      <  1:439512‑M1/251‑243 (MQ=255)
                                                                                                                                                                                                                                                 gCGACTATaatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagggttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagatatcgcggcgttctttggagagctcggcatttttgatgatgtaatcgtcaacgcgatcttcatagtcacctttcatgctggcgatgat     <  1:442270‑M1/251‑244 (MQ=255)
                                                                                                                                                                                                                                                    agtataatttcgattggttctcggatgaggatgatgacggtgatacacagataagaatcaatgagcaatgagagttgccgggcaaccctcattgaataaaacgggaattacttcgcttcgccgtttttcatttcgcccatagctttcagctttttggagatatcgcggcgttctttggagatctcggcatttttgatgatgtaatcgtcaacgcgatctgcatagtcacctttcatgctgtcgatgatgcc  <  2:314700‑M1/249‑247 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
AAATAAGATGGCATTACGGAATAGCCAGACTTTCTCGTGTTAACGTCCGACAAAATGCGTCATTGACAGCGGTTTCTGTGGATTGACACATTTCATTCAGATAATGAATTAATGCTACTGCCGGAACAACCAGCAAACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAG  >  NC_000913/2065903‑2066397

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