Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 390,251 Δ14,565 bp [yaiT][dgcC] 16 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 390251 404815 14565 45 [0] [0] 45 [yaiT]–[dgcC] [yaiT],ytiB,insF2,insE2,yaiT,iprA,ampH,sbmA,yaiW,yaiY,yaiZ,ddlA,iraP,phoA,psiF,[dgcC]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 3902500 (0.000)45 (0.790) 42/488 0.3 100% intergenic (‑523/‑1) hemB/yaiT porphobilinogen synthase/putative autotransporter YaiT
?NC_000913 404816 = 0 (0.000)coding (1114/1116 nt) dgcC diguanylate cyclase DgcC

GGCGATGAGTTTGCAGTAATCATGTCCGGTACGCCAGCTGAGAGCGCCATTACCGCCATGTTACGGGTGCATGAAGGGCTAAATACATTACGTTTGCCGAATACGCCACAGGTAACTTTACGGATTAGTGTGGGGGTTGCGCCGCTGAACCCACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTACCGCCTGAGCGGCAAATTTATACGCCTGGGCGCGTGGC  >  NC_000913/404570‑405060
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
gtgatatttattatgcgtatgcttcaaaaacaaaattatctgcaacatgttgaaatactgtgcttttatgaatttgatgcgtgtttttctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctaggcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGAc                                                                                                                                                                                                                                                   <  2:266511‑M2/4‑1 (MQ=255)
            atgcgtatgcttcaaaaacaaaattatctgcaacatgttgaaatactgtgcttttatgaatttgatgcgtgtttttctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGt                                                                                                                                                                                                                                      >  2:370137‑M2/235‑251 (MQ=255)
                                   ttatctgcaacatgttgaaatagtgtgcttttatgaatttgatgcgtgtttttctccataaattatctatgtccacattcggacttaggggatagaataattgaacaattcgtctgtaccgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAgtttttt                                                                                                                                                                                                               >  2:498285‑M2/212‑250 (MQ=255)
                                   ttatctgcaacatgttgaaatactgtgcttttatgaatttgatgcgtgtttttctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTc                                                                                                                                                                                                               >  1:49909‑M2/212‑251 (MQ=255)
                                         gcaacatgttgaaatactgtgcttttatgaatttgatgcgtgtttttctccataaattatatatgtccacattcggtcttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTg                                                                                                                                                                                                             <  2:176461‑M2/42‑1 (MQ=255)
                                             catgttgaaatactgtgcttttatgaatttgatgcgtgtttttctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTcc                                                                                                                                                                                                     >  1:428610‑M2/202‑251 (MQ=255)
                                                   gaaatactgtgcttttatgaatttgatgcgtgtttttctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATaca                                                                                                                                                                                                <  2:428610‑M2/55‑1 (MQ=255)
                                                    aaatactgtgcttttatgaatttgatgcgtgtttttctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGtc                                                                                                                                                                                                                                     >  1:164277‑M2/195‑212 (MQ=255)
                                                    aaatactgtgcttttatgaatttgatgcgtgtttttctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggagcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGtc                                                                                                                                                                                                                                     <  2:164277‑M2/18‑1 (MQ=255)
                                                       tactgtgcttttatgaatttgatgcgtgtttttctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCg                                                                                                                                                                                           >  1:446708‑M2/192‑251 (MQ=255)
                                                             gcttttatgaatttgatgcgtgtttttctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCg                                                                                                                                                                                     >  1:474351‑M2/186‑251 (MQ=255)
                                                                tttatgaatttgatgcgtgtttttctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCt                                                                                                                                                                                   <  1:370137‑M2/68‑1 (MQ=255)
                                                                              gcgtgtttttctccataacttatatatgtccccattcggacttaggggaaagaataattgaaccattcgtatgtaacgccgcataatagttagcgcgaaacatcatatgttttctatgcattgataattgatggatcaacttattacgtcactgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTgcagc                                                                                                                                                                     >  2:359079‑M2/169‑250 (MQ=255)
                                                                                gtgtttttctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACg                                                                                                                                                                  <  2:195118‑M2/85‑1 (MQ=255)
                                                                                       tctccataaattatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGcc                                                                                                                                                            <  1:299327‑M2/91‑1 (MQ=255)
                                                                                               aagtatatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataagtgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAATCCTTTCTCTTc                                                                                                                                                   <  2:67701‑M2/100‑1 (MQ=255)
                                                                                                attatatatgtccccattcggacttaggggaaagactaattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgctttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTcc                                                                                                                                                  >  2:361763‑M2/151‑251 (MQ=255)
                                                                                                   atatatgtccacattcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGt                                                                                                                                               <  1:105768‑M2/104‑1 (MQ=255)
                                                                                                     atatgtccacattcggacttaggggaaagaataattcaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccatgaggagggatgccaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTAc                                                                                                                                             >  2:264139‑M2/146‑251 (MQ=255)
                                                                                                           ccacattcggacttaggggaaagaataattgaaccattagtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTAc                                                                                                                                       >  1:465085‑M2/140‑251 (MQ=255)
                                                                                                               attcggacttaggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtcactgaggagggatgacaaTGACGTCCGGCGAAAGTCTTCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCt                                                                                                                                   >  2:48197‑M2/136‑251 (MQ=255)
                                                                                                                          ggggaaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctctgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATggtggt                                                                                                                        >  2:62514‑M2/125‑251 (MQ=255)
                                                                                                                              aaagaataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCt                                                                                                                    <  2:474351‑M2/131‑1 (MQ=255)
                                                                                                                                   ataattgaaccattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGAt                                                                                                                                                                                                           <  1:509987‑M2/44‑1 (MQ=255)
                                                                                                                                   ataattgaaccattcgtctgtaacgcagaatactcgttagcgcgaaacataatatgttttctatgcattgataattgatggataaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGAt                                                                                                                                                                                                           >  2:509987‑M2/116‑159 (MQ=255)
                                                                                                                                           accattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcatcttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTGCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGa                                                                                                       <  1:498285‑M2/144‑1 (MQ=255)
                                                                                                                                             cattcgtctgtaacgcagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGAcc                                                                                                     >  1:19075‑M2/106‑251 (MQ=255)
                                                                                                                                             cattcgtctgtaacgcagcataatcgttagagcgaaccataatctgctttctatgcattgctaattgatggatccacttattacatccctgaggagggatgacaaTGACGTCCGGCGAACGTCATCAGGATTTGCTGAGTTTTTATGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGAcc                                                                                                     >  2:201415‑M2/106‑251 (MQ=255)
                                                                                                                                                            cagcataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGcc                                                                                      >  1:126398‑M2/91‑251 (MQ=255)
                                                                                                                                                                ataatcgttagcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCg                                                                                  >  2:429078‑M2/87‑251 (MQ=255)
                                                                                                                                                                   atcgttagcgcgaaacataatatgttttctatgcattgatcattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGCAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGAt                                                                               >  1:323933‑M2/84‑251 (MQ=255)
                                                                                                                                                                         agcgcgaaacataatatgttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTtctc                                                                         >  1:325540‑M2/78‑251 (MQ=255)
                                                                                                                                                                                          gttttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCAtt                                                        >  2:105945‑M2/61‑251 (MQ=255)
                                                                                                                                                                                            tttctatgcattgataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCAGGCGCAAGTCATCAGGATTTGCTGCGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTGCGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCAtttt                                                      >  1:233411‑M2/59‑251 (MQ=255)
                                                                                                                                                                                                         ataattgatggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGTATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCAtt                                         <  1:264139‑M2/206‑1 (MQ=255)
                                                                                                                                                                                                                 tggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAAc                                              >  1:11557‑M2/38‑238 (MQ=255)
                                                                                                                                                                                                                 tggatcaacttattacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAAc                                              <  2:11557‑M2/201‑1 (MQ=255)
                                                                                                                                                                                                                    atcaacttattacgtccctgaggagggatgacaaTGACGTTCGGCGAAAGTCATCAGGATTTGCTGAGTGTGTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCAGATCTTTCAGGGCCCCCGGATGTTCTCCCGGTGCCAGCACCATTTTTGCGGAACCCATTACCGCCTGAGc                              <  1:359079‑M2/217‑1 (MQ=255)
                                                                                                                                                                                                                           tattacgtccctgaggagggatgacaaTGACGTGCGGAGAAAGTCATCCGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACAGCTTCAATGGTGGTGCCTCCCGGTGTGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTttgcggaacccattaccagttcggaagagcgg                       >  2:311358‑M2/28‑236 (MQ=255)
                                                                                                                                                                                                                           tattacggccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTAcc                                      <  1:311358‑M2/209‑1 (MQ=255)
                                                                                                                                                                                                                            attacgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTACCGCCTGAGCGGCAAAt                       <  1:83646‑M2/224‑1 (MQ=255)
                                                                                                                                                                                                                               acgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGCTTTGCTGAGTTTTTCTGATTTTTCCATACCCTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTACCGCCTGAGCGGCAAATTtat                   >  1:307276‑M2/24‑251 (MQ=255)
                                                                                                                                                                                                                               acgtccctgaggagggatgacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTACCGCCTGAGCGGCAAATTtat                   >  1:358834‑M2/24‑251 (MQ=255)
                                                                                                                                                                                                                                       gaggagggatgacaaTGCCGTCCGGGGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCCTCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCACTGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTACCGCCTGAGCGGCAAATTTATACGCCTgg           >  2:484972‑M2/16‑251 (MQ=255)
                                                                                                                                                                                                                                                 gacaaTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTACCGCCTGCGCGGCAAATTTATACGCCTGGGCGCGTGGc  >  2:301352‑M2/6‑250 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
GGCGATGAGTTTGCAGTAATCATGTCCGGTACGCCAGCTGAGAGCGCCATTACCGCCATGTTACGGGTGCATGAAGGGCTAAATACATTACGTTTGCCGAATACGCCACAGGTAACTTTACGGATTAGTGTGGGGGTTGCGCCGCTGAACCCACAAATGAGTCACTATCGTGAGTGGTTGAAATCGGCAGATTTGGCGCTTTACAAAGCAAAGAAAGCCGGACGTAACCGCACCGAAGTGGCGGCCTGACGTCCGGCGAAAGTCATCAGGATTTGCTGAGTTTTTCTGATTTTTCCATACACTTCGTCATCGCTTCGATCACTGCAGCACGGAAGCCTTTCTCTTCCAGTACGCGTACCGCTTCAATGGTGGTGCCTCCCGGTGAGCAGACCATATCTTTCAGTGCCCCCGGATGTTCTCCCGTTTCCAGCACCATTTTTGCGGAACCCATTACCGCCTGAGCGGCAAATTTATACGCCTGGGCGCGTGGC  >  NC_000913/404570‑405060

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