Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 1,640,057 Δ7,591 bp [rrrQ][dicC] 18 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 1640057 1647647 7591 64 [0] [0] 64 [rrrQ]–[dicC] [rrrQ],[ydfR],essQ,ynfU,cspB,cspF,ynfS,ydfT,ydfU,rem,hokD,relE,relB,ydfV,flxA,ydfW,ydfX,[dicC]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 16400560 (0.000)64 (0.830) 58/494 0.1 100% coding (2/534 nt)
coding (310/312 nt)
rrrQ
ydfR
putative lysozyme
protein YdfR
?NC_000913 1647648 = 0 (0.000)coding (203/231 nt) dicC putative DNA‑binding transcriptional regulator DicC

GCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGCCTGCTCCTGTGTCATTCCTGGTCGACATATCCAGGCATCCAGAGCCTCACGAGCCTGTTCAGGAGTGATTTTCATTGTTCAACCGCCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGACGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGAT  >  NC_000913/1647399‑1647887
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                               
atcacgctcaatggcgttgacctggtcgcatttttccttcgacagtttcatattgggaaaaacggtttttccatccaccactgtggcaccccgacagatggtccatatgccagaaccatcgcggtatgccattgtgtggttaccttctttttcgtccagaaactggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcacc                                                                                                                                                                                                                                                <  1:33607‑M2/2‑1 (MQ=255)
  cacgctcaatggcgttggcctggtcgcatttttccttcgacagtttcatattgggaaaaacggtttttccatccaccactgtggcaccccgacagatggtccatatgccagaaccatcgcggtatgccattgcgtggttgccttctttttcgtccagaaacgggtcaagtatctgagtagcagattcgccatcatcaatcatcgccagaacggcagccgacaggctgtatctgatcttttgtggtcacccg                                                                                                                                                                                                                                              <  2:537059‑M2/4‑1 (MQ=255)
   acgctcaatggcgttgacctggtcgcatttttccttcgacagtttcatattgggaaaaacggtttttccatccaccactgtggcaccccgacagatggtccatatgccagaaccatcgcggtatgccattgtgtggttaccttctttttcgtccagaaactggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcacccgc                                                                                                                                                                                                                                             >  1:589531‑M2/247‑251 (MQ=255)
     gctcaatggcgttgacctggtcgcatttttccttcgacagtttcatattgggaaaaacggtttttccatccaccactgtggcaccccgacagatggtccatatgccagaaccatcgcggtatgccattgtgtggttaccttctttttcgtccagaaactggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcacccgccc                                                                                                                                                                                                                                           <  2:124552‑M2/7‑1 (MQ=255)
        caatggcgttgacctggtcgcatttttccttcgacagtttcatattgggaaaaacgttttttcccaccaccactgtggcgccccgacagatggtccatatgccagaaccatcgcggtatgccattgtgtggttaccttctttttcgtccagaaacgggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagagcggcagcctacaggccgtatctgatcttttgtgttcacccgcccg                                                                                                                                                                                                                                          <  2:37389‑M2/8‑1 (MQ=255)
        caatggcgttgacctggtcgcatttttccttcgacagtttcatattgggaaaaacggtttttccatccaccactgtggcaccccgacagatggtccatatgccagaaccatcgcggtatgccattgtgtggttaccttctttttcgtccagaaactggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCt                                                                                                                                                                                                                                        <  1:326928‑M2/10‑1 (MQ=255)
          atggcgttgacctggtcgcatttttccttcgacagtttcatattgggaaaaacggtttttccatccaccactgtggcaccccgacagatggtccatatgccagaaccatcgcggtatgccattgtgtggttaccttctttttcgtccagaaactggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTc                                                                                                                                                                                                                                      <  2:620431‑M2/12‑1 (MQ=255)
            ggcgttgacctggtcgcatttttccttcgacagtttcatattgggaaaaacgttttttccatccaccactgtggcaccccgacgtatggtccatatgccagatccatcgcggtatgccattgtgtggttaccttctttttcgtccagaatctgggcaagtatctgaggagcagatgtgccagccccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCg                                                                                                                                                                                                                                     <  1:249693‑M2/13‑1 (MQ=255)
             gcgttgacctggtcgcatttttccttcgacagtttcatattgggaaaaacggtttttccatccaccactgtggcaccccgacagatggtccatatgccagaaccatcgcggtatgccattgtgtggttaccttctttttcgtccagaaactggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTc                                                                                                                                                                                                                                   >  1:497082‑M2/237‑251 (MQ=255)
                  gacctggtcgcatttttccttcgacagtttcatattgggaaaaacggtttttccatccacccctgtggcaccccgacagatggtccatatgccagaaccatcgcggtatgccattgtgtggttaccttcttttccgtccagaaactggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgacaggccgtagctgatcttgtgtgttcaCCCGCCCGCTTCGTCTTACg                                                                                                                                                                                                                              <  2:3702‑M2/20‑1 (MQ=255)
                  gacctggtcgcatttttccttcgacagtttcatattgggaaaaacggtttttccatccaccactgtggcaccccgacagatggtccatatgccagaaccatcgcggtatgccattgtgtggttaccttctttttcgtccagaaactggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTCTTACg                                                                                                                                                                                                                              <  1:534054‑M2/20‑1 (MQ=255)
                          cgcatttttccttcgacagtttcatattgggaaaaccgggttttccatccacccctgtggcaccccgacagatggtccatatgccagaaccatcgcggtatgccattgtgtggttaccttctttttcgtccagaaactggtcaagtatctgaggagcagatgcgccagcaccaatcagcgcccgaacggcagccgacaggccgtatctgctcttttgtgttcaCCCGCACGCTTCGTCTTACGATATtcat                                                                                                                                                                                                                      >  2:574895‑M2/224‑251 (MQ=255)
                           gcatttttccttcgacagtttcatattgggaaaaacggtttttccatccaccactgtggcaccccgacagatggtccatatgccagaaccatcgcggtatgccattgtgtggttaccttctttttcgtccagaaactggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTCTTACGTTATtcatc                                                                                                                                                                                                                     >  2:395088‑M2/223‑251 (MQ=255)
                           gcatttttccttcgacagtttcatattgggaaaaacggtttttccatccaccactgtggcaccccgacagatggtccatatgccagaaccatcgcggtatgccattgtgtggttaccttctttttcgtccagaaactggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTCTTACGATATtcatc                                                                                                                                                                                                                     <  2:11560‑M2/29‑1 (MQ=255)
                             atttttccttcgacagtttcatattgggaaaaacggtttttccatccaccactgtggcaccccgacagatggtccatatgccagaaccatcgcggtatgccattgtgtggttaccttctttttcgtccagaaactggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTCTTACGATATTcatcat                                                                                                                                                                                                                   >  1:500203‑M2/221‑251 (MQ=255)
                                  tccttcgccagtttcatattgggaaaaacggtttttccatccaccactgtggcaccccgacagatggtccatatgccagaaccatcgcggtatgccattgtgtggttaccttcttttgcgtccagaaactggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTCTTACGATATTCATCATAAAc                                                                                                                                                                                                               <  2:340898‑M2/35‑1 (MQ=255)
                                  tccttcgacagtttcatattgggaaaaacggtttttccatccaccactgtggcaccccgacagatggtccatatgccagaaccatcgcggtatgccattgtgtggttaccttctttttcgtccagaaactggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTCTTACGATATTCATCATAAAc                                                                                                                                                                                                               <  1:405717‑M2/35‑1 (MQ=255)
                                   ccttcgacagtttcatattgggaaaaacggtttttccatccaccactgtggcaccccgacagatggtccatatgccagaaccatcgcggtatgccattgtgtggttaccttctttttcgtccagaaactggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTCTTACGATATTCATCATAAACt                                                                                                                                                                                                              <  1:171318‑M2/36‑1 (MQ=255)
                                    cttcgacagtttcatattgggaaaaatggtttttccatccaccactgtggcaccccgacagatggtccatatgccagaaccatcgcggtatgccattgtgtggttaccttctttttcgtccagaaactggtcaagtatctgaggagcagatgcgccagccccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTCTTACGATATTAATCATAAACttt                                                                                                                                                                                                            >  2:265878‑M2/214‑251 (MQ=255)
                                      tcgacagtttcatattgggaaaaacggtttttccatccaccactgtggcaccccgacagatggtccatatgccagaaccatcgcggtatgccattgtgtggttaccttctttttcgtccagaaactggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgctcaCCCGCCCTCTTCGTCTTACGATATTCATCATAAACTTTgg                                                                                                                                                                                                          >  1:168995‑M2/212‑251 (MQ=255)
                                         acagtttcatattgggaaaaacggtttttccatccaccactgtggcaccccgacagatggtccatatgccagaaccatcgcggtatgccattgtgtggttaccttctttttcgtccagaaactggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGAt                                                                                                                                                                                                       <  1:98183‑M2/43‑1 (MQ=255)
                                          cagtttcatattgggaaaaacggtttttccatcccccactgtggcaccccgacagatggtccatatgccagaaccatcgcggtatgccattgtgtggttaccttctttttcgtccagaaactggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgacagtccgtatctgatcttttgtgttcaCCCGCCCTCTTCGTCTTACGATATTCATCATAAACTTTGGGATc                                                                                                                                                                                                      <  2:504282‑M2/44‑1 (MQ=255)
                                            gtttcatattgggaaaaacggtttttccatccaccactgtggcaccccgacagatggtccatatgccagaaccatcgcggtatgccattgtgtggttaccttctttttcgtccagaaactggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCAt                                                                                                                                                                                                    >  1:345114‑M2/206‑251 (MQ=255)
                                               tcatattgggaaaaacggtttttccatccaccagtttggcaccccgacagatggtccatatgccagaaccatcgcggtatgccattgtgtggttaccttctttttcgtccagaaactggtcaagtatctgcggagcagatgcgccagcaccaatcagcgccagaacggcaaccgccaggccgtatctgatctcttccgttcaCCCGCCCGCTTCGTCTTACGATATTCATCATACActttgggaccagact                                                                                                                                                                                                 >  2:410109‑M2/203‑244 (MQ=255)
                                                catattgggaaaaacggtttttccatccaccactgttgcaccccgacagatggtccatatgccagaaccatcgcggtatgccattgtgtggttaccttctttttcgtccagaaactggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTg                                                                                                                                                                                                <  2:309888‑M2/50‑1 (MQ=255)
                                                      gggaaaaacggtttttccatccaccactgtggcaccccgacagatggtccatatgccagaaccatcgcggtatgccattgtgtggttaccttctttttcgtccagaaactggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTCTTACGATATTCATCATAAGCTTTGGGAt                                                                                                                                                                                                       <  1:397335‑M2/43‑1 (MQ=255)
                                                      gggaaaaacggtttttccatccaccactgtggcaccccgacagatggtccatatgccagaaccatcgcggtatgccattgtgtggttaccttctttttcgtccagaaactggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTCTTACGAAATTCATCATAACCTTTGGGAt                                                                                                                                                                                                       >  2:397335‑M2/196‑238 (MQ=255)
                                                       ggaaaaacggtttttccatccaccactgtggcaccccgacagatggtccatatgccagaaccatcgcggtatgccattgtgtggttaccttctttttcgtccagaacctggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTcc                                                                                                                                                                                         >  2:304053‑M2/195‑251 (MQ=255)
                                                        gaaaaacggtttttccatccaccactgtggcaccccgacagatggtccatatgccagaaccatcgcggtatgccattgtgtggttaccttctttttcgtccagaaactggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTccc                                                                                                                                                                                        >  1:378524‑M2/194‑251 (MQ=255)
                                                        gaaaaacggtttttccatccaccactgtggcaccccgacagatggtccatatgccagaaccatcgcggtatgccattgtgtggttaccttctttttcgtccagaaactggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTccc                                                                                                                                                                                        >  1:96593‑M2/194‑251 (MQ=255)
                                                                    ttccatcccccactgtggcaccccgacagatggtccatatgccagaaccatcgcggtatgccattgtgtggttgccttcttttgcgtccagaaactggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGcc                                                                                                                                                                             <  2:597113‑M2/69‑1 (MQ=255)
                                                                           caccactgtggcaccccgacagatggtccatatgccagaaccatcgcggtatgccattgtgtggttaccttctttttcgtccagaaactggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGt                                                                                                                                                                       >  1:107540‑M2/175‑249 (MQ=255)
                                                                             ccactgtggcaccccgacagatggtccatatgccagaaccatcgcggtatgccattgtgtggttaccttctttttcgtccagaaactggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGAc                                                                                                                                                                   >  1:619132‑M2/173‑251 (MQ=255)
                                                                                  gtggcaccccgacagatggtccatatgccagaaccatcgcggtatgccattgtgtggttaccttctttttcgtccagaaactggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGACGCATc                                                                                                                                                              <  2:96593‑M2/84‑1 (MQ=255)
                                                                                  gtggcaccccgacagatggtccatatgccagaaccatcgcggtatgccattgtgtggttaccttctttttcgtccagaaactggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGACGCATc                                                                                                                                                              >  2:20564‑M2/168‑251 (MQ=255)
                                                                                            gacagatggtccatatgccagaaccatcgcggtatgccattgtgtggttaccttctttttcgtccagaaactggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAAATGCCTCCTGTAGACGCATCGCGCGACCtt                                                                                                                                                    >  2:399730‑M2/158‑251 (MQ=255)
                                                                                                  tgctccatatgccagaaccatcgcggtaggccattgtgtggttaccttctgtttcgtccagaaacgggtcaagtatctgaggagcagaggcgccagccccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCTCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGACGCATCGCGCGACCTTCGGGa                                                                                                                                               <  2:321855‑M2/99‑1 (MQ=255)
                                                                                                          atgccagaaccatcgcggtatgccattgtgtggttcccttctttttcgtccagaaactggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGACGCATCGCGCGACCTTCGGGAACTAAATc                                                                                                                                       <  2:349525‑M2/107‑1 (MQ=255)
                                                                                                                 aaccatcgcggtatgccattgtgtggttaccttctttttcgtccagaaactggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGACGCATCGCGCGACCTTCGGGAACTCAACCCCCTTTcc                                                                                                                               >  2:598761‑M2/137‑251 (MQ=255)
                                                                                                                       cgcggtatgccattgtgtggttaccttctttttcgtccagaaactggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGACGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCta                                                                                                                          <  1:546042‑M2/120‑1 (MQ=255)
                                                                                                                           gtatgccattgtgtggttaccttctttttcgtccagaaactggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGACGCATCGCGCGACCTTCGGGAACAAAATCCCCTTTCCAGCTATaa                                                                                                                       >  2:331271‑M2/127‑249 (MQ=255)
                                                                                                                            tatgccattgtgtggttaccttctttttcgtccagaaactggtcaagtatctgaggagcagattcgccagcaccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaACAACCCGCTTCGTCTTACGCTATTCATCATAAACTTTGGGATGATACTGACGCTCCCCGCCAGATGCCCCCTGTAGACGAATCGAGCGACCTTCGGAAACAAACACCCCTTTCCAGATATAAAGc                                                                                                                    >  2:74310‑M2/126‑251 (MQ=255)
                                                                                                                               gccattgtgtggttaccttctttttcgtccagaaactggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGTGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGACGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGa                                                                                                                  >  1:631351‑M2/123‑250 (MQ=255)
                                                                                                                               gccattgtgtggttaccttctttttcgtccagaaactggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGACGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGaa                                                                                                                 <  1:331271‑M2/129‑1 (MQ=255)
                                                                                                                                    tgtgtgtgtacattcgtgttcgtccagaatctggtcaagtatctgaggagcagatccgccagcacaaagctgcgccagaccgccagccgacaggccgtatctgatcttttgtgtgcaCCCGCCCGCTTCGTCTTACTATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGACGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCa                                                                                                             <  2:278981‑M2/133‑1 (MQ=255)
                                                                                                                                     gtgtggttaccttctttttcgtccagaaactggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGGCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGACGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCaa                                                                                                            >  2:567018‑M2/117‑250 (MQ=255)
                                                                                                                                       gtggttaccttctttttcgtccagaaactggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGACGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACg                                                                                                         <  1:353751‑M2/137‑1 (MQ=255)
                                                                                                                                                          tccagaaactggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGACGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTtgt                                                                                         <  2:619132‑M2/153‑1 (MQ=255)
                                                                                                                                                            cagaaactggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGACGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAg                                                                                    <  1:34591‑M2/158‑1 (MQ=255)
                                                                                                                                                               aaacgggtcaagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGACGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAgttt                                                                                 <  2:176781‑M2/161‑1 (MQ=255)
                                                                                                                                                                        aagtatctgaggagcagatgcgccagcaccaatcagcgccagaacggcagccgccaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGACGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTg                                                                         >  1:270303‑M2/82‑250 (MQ=255)
                                                                                                                                                                         agaatcggaggagcagatgcgccagccccaatcatcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGACGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGaa                                                                       <  2:133583‑M2/171‑1 (MQ=255)
                                                                                                                                                                            atctgaggagcagatgcgccagcaccaatcagcgccagaccggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGACGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAAcc                                                                     <  2:497082‑M2/173‑1 (MQ=255)
                                                                                                                                                                             cctgagggtccgatgcccctacaccaatcagccccagaacggcagccgtcaggccgtatctgagcttttgttttccCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCATATGCCTCCTGTAGACGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGa                                                                   <  2:591679‑M2/175‑1 (MQ=255)
                                                                                                                                                                                  ggagcagatgcgccagcaccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGACGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATAc                                                              <  1:395088‑M2/180‑1 (MQ=255)
                                                                                                                                                                                   gagcagatgcgccagccccaatcagcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGACGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATTCCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACa                                                             <  1:104737‑M2/181‑1 (MQ=255)
                                                                                                                                                                                                    ccaatcagctccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGACGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTa                                            <  2:46147‑M2/198‑1 (MQ=255)
                                                                                                                                                                                                          atcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGACGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCAtt                                      <  2:252520‑M2/204‑1 (MQ=255)
                                                                                                                                                                                                           gcgccagaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGACGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCAttt                                     <  1:37745‑M2/205‑1 (MQ=255)
                                                                                                                                                                                                                agaacggcagccgacaggccgtatctgatcttttgtgttcaCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGACGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAac                                >  1:83150‑M2/42‑251 (MQ=255)
                                                                                                                                                                                                                                          gatcttttgtgttccCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGACGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCAtaactaac       <  2:14146‑M2/235‑1 (MQ=255)
                                                                                                                                                                                                                                             cttttgtgttcaCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGACGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACa      >  1:593281‑M2/13‑248 (MQ=255)
                                                                                                                                                                                                                                               tttgtgttcaCCCGCCCGCTTTGTCTTTCGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGACGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAagat  <  2:105378‑M2/240‑1 (MQ=255)
                                                                                                                                                                                                                                               tttgtgttcaCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGACGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACa      <  1:270367‑M2/236‑1 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                               
GCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGCCTGCTCCTGTGTCATTCCTGGTCGACATATCCAGGCATCCAGAGCCTCACGAGCCTGTTCAGGAGTGATTTTCATTGTTCAACCGCCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGACGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGAT  >  NC_000913/1647399‑1647887

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