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

GTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTTGAAATACTGTGCTTTTATGAATTTGATGCGTGTTTTTCTCCATAAATTATATATGTCCACATTCGGACTTAGGGGAAAGAATAATTGAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATAC  >  NC_000913/390005‑390495
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
gtgATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTTGAAATACTGTGCTTTTATGAATTTGATGCGTGTTTTTCTCCATAAATTATATATGTCCACATTCGGACTTAGGGGAAAGAATAATTGAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTAGGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgac                                                                                                                                                                                                                                                   <  2:266511‑M1/250‑5 (MQ=255)
            atGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTTGAAATACTGTGCTTTTATGAATTTGATGCGTGTTTTTCTCCATAAATTATATATGTCCACATTCGGACTTAGGGGAAAGAATAATTGAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagt                                                                                                                                                                                                                                      >  2:370137‑M1/1‑234 (MQ=255)
                                   ttATCTGCAACATGTTGAAATAGTGTGCTTTTATGAATTTGATGCGTGTTTTTCTCCATAAATTATCTATGTCCACATTCGGACTTAGGGGATAGAATAATTGAACAATTCGTCTGTACCGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttt                                                                                                                                                                                                               >  2:498285‑M1/1‑211 (MQ=255)
                                   ttATCTGCAACATGTTGAAATACTGTGCTTTTATGAATTTGATGCGTGTTTTTCTCCATAAATTATATATGTCCACATTCGGACTTAGGGGAAAGAATAATTGAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttc                                                                                                                                                                                                               >  1:49909‑M1/1‑211 (MQ=255)
                                         gCAACATGTTGAAATACTGTGCTTTTATGAATTTGATGCGTGTTTTTCTCCATAAATTATATATGTCCACATTCGGTCTTAGGGGAAAGAATAATTGAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctg                                                                                                                                                                                                             <  2:176461‑M1/247‑43 (MQ=255)
                                             cATGTTGAAATACTGTGCTTTTATGAATTTGATGCGTGTTTTTCTCCATAAATTATATATGTCCACATTCGGACTTAGGGGAAAGAATAATTGAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttcc                                                                                                                                                                                                     >  1:428610‑M1/1‑201 (MQ=255)
                                                   gAAATACTGTGCTTTTATGAATTTGATGCGTGTTTTTCTCCATAAATTATATATGTCCACATTCGGACTTAGGGGAAAGAATAATTGAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccataca                                                                                                                                                                                                <  2:428610‑M1/250‑56 (MQ=255)
                                                    aaaTACTGTGCTTTTATGAATTTGATGCGTGTTTTTCTCCATAAATTATATATGTCCACATTCGGACTTAGGGGAAAGAATAATTGAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtc                                                                                                                                                                                                                                     >  1:164277‑M1/1‑194 (MQ=255)
                                                    aaaTACTGTGCTTTTATGAATTTGATGCGTGTTTTTCTCCATAAATTATATATGTCCACATTCGGACTTAGGGGAAAGAATAATTGAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGAGCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtc                                                                                                                                                                                                                                     <  2:164277‑M1/212‑19 (MQ=255)
                                                       tACTGTGCTTTTATGAATTTGATGCGTGTTTTTCTCCATAAATTATATATGTCCACATTCGGACTTAGGGGAAAGAATAATTGAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcg                                                                                                                                                                                           >  1:446708‑M1/1‑191 (MQ=255)
                                                             gCTTTTATGAATTTGATGCGTGTTTTTCTCCATAAATTATATATGTCCACATTCGGACTTAGGGGAAAGAATAATTGAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcg                                                                                                                                                                                     >  1:474351‑M1/1‑185 (MQ=255)
                                                                tttATGAATTTGATGCGTGTTTTTCTCCATAAATTATATATGTCCACATTCGGACTTAGGGGAAAGAATAATTGAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgct                                                                                                                                                                                   <  1:370137‑M1/250‑69 (MQ=255)
                                                                              gCGTGTTTTTCTCCATAACTTATATATGTCCCCATTCGGACTTAGGGGAAAGAATAATTGAACCATTCGTATGTAACGCCGCATAATAGTTAGCGCGAAACATCATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCACTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagc                                                                                                                                                                     >  2:359079‑M1/1‑168 (MQ=255)
                                                                                gtgtTTTTCTCCATAAATTATATATGTCCACATTCGGACTTAGGGGAAAGAATAATTGAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacg                                                                                                                                                                  <  2:195118‑M1/251‑86 (MQ=255)
                                                                                       tctcCATAAATTATATATGTCCACATTCGGACTTAGGGGAAAGAATAATTGAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcc                                                                                                                                                            <  1:299327‑M1/250‑92 (MQ=255)
                                                                                               aagtatatatgtCCACATTCGGACTTAGGGGAAAGAATAATTGAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAAGTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaatcctttctcttc                                                                                                                                                   <  2:67701‑M1/248‑101 (MQ=255)
                                                                                                aTTATATATGTCCCCATTCGGACTTAGGGGAAAGACTAATTGAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGCTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttcc                                                                                                                                                  >  2:361763‑M1/1‑150 (MQ=255)
                                                                                                   atatatGTCCACATTCGGACTTAGGGGAAAGAATAATTGAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagt                                                                                                                                               <  1:105768‑M1/251‑105 (MQ=255)
                                                                                                     atatGTCCACATTCGGACTTAGGGGAAAGAATAATTCAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCATGAGGAGGGATGCCAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtac                                                                                                                                             >  2:264139‑M1/1‑145 (MQ=255)
                                                                                                           ccACATTCGGACTTAGGGGAAAGAATAATTGAACCATTAGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtac                                                                                                                                       >  1:465085‑M1/1‑139 (MQ=255)
                                                                                                               aTTCGGACTTAGGGGAAAGAATAATTGAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCACTGAGGAGGGATGACAAtgacgtccggcgaaagtcttcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgct                                                                                                                                   >  2:48197‑M1/1‑135 (MQ=255)
                                                                                                                          ggggAAAGAATAATTGAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTCTGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggtggt                                                                                                                        >  2:62514‑M1/1‑124 (MQ=255)
                                                                                                                              aaaGAATAATTGAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggtggtgcct                                                                                                                    <  2:474351‑M1/251‑132 (MQ=255)
                                                                                                                                   ataatTGAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgat                                                                                                                                                                                                           <  1:509987‑M1/159‑45 (MQ=255)
                                                                                                                                   ataatTGAACCATTCGTCTGTAACGCAGAATACTCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATAAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgat                                                                                                                                                                                                           >  2:509987‑M1/1‑115 (MQ=255)
                                                                                                                                           aCCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCATCTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctctgccagtacgcgtaccgcttcaatggtggtgcctcccggtgagcaga                                                                                                       <  1:498285‑M1/251‑145 (MQ=255)
                                                                                                                                             cATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggtggtgcctcccggtgagcagacc                                                                                                     >  1:19075‑M1/1‑105 (MQ=255)
                                                                                                                                             cATTCGTCTGTAACGCAGCATAATCGTTAGAGCGAACCATAATCTGCTTTCTATGCATTGCTAATTGATGGATCCACTTATTACATCCCTGAGGAGGGATGACAAtgacgtccggcgaacgtcatcaggatttgctgagtttttatgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggtggtgcctcccggtgagcagacc                                                                                                     >  2:201415‑M1/1‑105 (MQ=255)
                                                                                                                                                            cagcaTAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggtggtgcctcccggtgagcagaccatatctttcagtgcc                                                                                      >  1:126398‑M1/1‑90 (MQ=255)
                                                                                                                                                                aTAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggtggtgcctcccggtgagcagaccatatctttcagtgcccccg                                                                                  >  2:429078‑M1/1‑86 (MQ=255)
                                                                                                                                                                   aTCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATCATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgcaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggtggtgcctcccggtgagcagaccatatctttcagtgcccccggat                                                                               >  1:323933‑M1/1‑83 (MQ=255)
                                                                                                                                                                         aGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggtggtgcctcccggtgagcagaccatatctttcagtgcccccggatgttctc                                                                         >  1:325540‑M1/1‑77 (MQ=255)
                                                                                                                                                                                          gTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggtggtgcctcccggtgagcagaccatatctttcagtgcccccggatgttctcccgtttccagcaccatt                                                        >  2:105945‑M1/1‑60 (MQ=255)
                                                                                                                                                                                            tttCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtcaggcgcaagtcatcaggatttgctgcgtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtgcgcgtaccgcttcaatggtggtgcctcccggtgagcagaccatatctttcagtgcccccggatgttctcccgtttccagcaccatttt                                                      >  1:233411‑M1/1‑58 (MQ=255)
                                                                                                                                                                                                         aTAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggtggtgcctcccggtgagcagaccatatctttcagtgcccccgtatgttctcccgtttccagcaccatttttgcggaacccatt                                         <  1:264139‑M1/251‑207 (MQ=255)
                                                                                                                                                                                                                 tGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggtggtgcctcccggtgagcagaccatatctttcagtgcccccggatgttctcccgtttccagcaccatttttgcggaac                                              >  1:11557‑M1/1‑37 (MQ=255)
                                                                                                                                                                                                                 tGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggtggtgcctcccggtgagcagaccatatctttcagtgcccccggatgttctcccgtttccagcaccatttttgcggaac                                              <  2:11557‑M1/238‑202 (MQ=255)
                                                                                                                                                                                                                    aTCAACTTATTACGTCCCTGAGGAGGGATGACAAtgacgttcggcgaaagtcatcaggatttgctgagtgtgtctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggtggtgcctcccggtgagcagaccagatctttcagggcccccggatgttctcccggtgccagcaccatttttgcggaacccattaccgcctgagc                              <  1:359079‑M1/251‑218 (MQ=255)
                                                                                                                                                                                                                           tattaCGTCCCTGAGGAGGGATGACAAtgacgtgcggagaaagtcatccggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtacagcttcaatggtggtgcctcccggtgtgcagaccatatctttcagtgcccccggatgttctcccgtttccagcaccatttttgcggaacccattaccagttcggaagagcgg                       >  2:311358‑M1/1‑27 (MQ=255)
                                                                                                                                                                                                                           tattaCGGCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggtggtgcctcccggtgagcagaccatatctttcagtgcccccggatgttctcccgtttccagcaccatttttgcggaacccattacc                                      <  1:311358‑M1/236‑210 (MQ=255)
                                                                                                                                                                                                                            attaCGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggtggtgcctcccggtgagcagaccatatctttcagtgcccccggatgttctcccgtttccagcaccatttttgcggaacccattaccgcctgagcggcaaat                       <  1:83646‑M1/250‑225 (MQ=255)
                                                                                                                                                                                                                               aCGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggctttgctgagtttttctgatttttccatacccttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggtggtgcctcccggtgagcagaccatatctttcagtgcccccggatgttctcccgtttccagcaccatttttgcggaacccattaccgcctgagcggcaaatttat                   >  1:307276‑M1/1‑23 (MQ=255)
                                                                                                                                                                                                                               aCGTCCCTGAGGAGGGATGACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggtggtgcctcccggtgagcagaccatatctttcagtgcccccggatgttctcccgtttccagcaccatttttgcggaacccattaccgcctgagcggcaaatttat                   >  1:358834‑M1/1‑23 (MQ=255)
                                                                                                                                                                                                                                       gaggagGGATGACAAtgccgtccggggaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcctcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcactggtggtgcctcccggtgagcagaccatatctttcagtgcccccggatgttctcccgtttccagcaccatttttgcggaacccattaccgcctgagcggcaaatttatacgcctgg           >  2:484972‑M1/1‑15 (MQ=255)
                                                                                                                                                                                                                                                 gACAAtgacgtccggcgaaagtcatcaggatttgctgagtttttctgatttttccatacacttcgtcatcgcttcgatcactgcagcacggaagcctttctcttccagtacgcgtaccgcttcaatggtggtgcctcccggtgagcagaccatatctttcagtgcccccggatgttctcccgtttccagcaccatttttgcggaacccattaccgcctgcgcggcaaatttatacgcctgggcgcgtggc  >  2:301352‑M1/1‑5 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
GTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTTGAAATACTGTGCTTTTATGAATTTGATGCGTGTTTTTCTCCATAAATTATATATGTCCACATTCGGACTTAGGGGAAAGAATAATTGAACCATTCGTCTGTAACGCAGCATAATCGTTAGCGCGAAACATAATATGTTTTCTATGCATTGATAATTGATGGATCAACTTATTACGTCCCTGAGGAGGGATGACAAATGCACTCCTGGAAAAAGAAACTTGTAGTATCACAATTAGCATTGGCTTGCACTCTGGCTATCACCTCTCAGGCTAATGCAGCAAACTATGATACCTGGACTTATATCGATAATCCCGTTACAGCACTTGATTGGGATCATATGGATAAGGCAGGCACTGTAGATGGCAACTATGTAAACTATAGTGGTTTTGTCTATTACAACAACACCAATGGTGATTTCGATCAGTCCTTTAACGGCGATAC  >  NC_000913/390005‑390495

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: