Predicted mutation
evidence seq id position mutation freq annotation gene description
MC JC NC_000913 4,001,645 Δ5 bp 100% coding (220‑224/951 nt) corA → magnesium/nickel/cobalt transporter

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 4001645 4001649 5 78 [0] [0] 78 corA magnesium/nickel/cobalt transporter

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 40016440 (0.000)70 (1.010) 53/254 NT 100% coding (219/951 nt) corA magnesium/nickel/cobalt transporter
?NC_000913 4001650 = 0 (0.000)coding (225/951 nt) corA magnesium/nickel/cobalt transporter

GTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGATG  >  NC_000913/4001504‑4001797
                                                                                                                                                  |                                                                                                                                                   
atgcagtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatct                                                                                                                                                    <  5:21486‑M2/2‑1 (MQ=255)
atgcagtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaacctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatct                                                                                                                                                    >  3:185782‑M2/147‑148 (MQ=255)
   cagtatggattgatcttgtcgaacccgacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccccccggaactggaagacatcgaagcatcggcatgtttctttgaagacgacgacggcctgcatCTCct                                                                                                                                                 <  6:41647‑M2/5‑1 (MQ=255)
       atggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttctt                                                                                                                                             >  5:275668‑M2/140‑148 (MQ=255)
       atggattgatcttctcgaagcggacgacgaccagccactgcgcgtacaatctgaacttggccagagccatgcaccccgcccggaactggaacacaacgaagcatcggcacgtttctttgaagacgacgacggcctgcgtCtccttcaa                                                                                                                                             >  6:219591‑M2/140‑146 (MQ=255)
         ggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttct                                                                                                                                           >  5:279802‑M2/138‑148 (MQ=255)
          gaatgatcttgtctaaccggacgacgtctagcgacttcgcgtacaacctgaactttgcaagagccttgcaacccgcacggaactttaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatGTCcttcttct                                                                                                                                           <  4:185782‑M2/11‑1 (MQ=255)
            ttgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttcttct                                                                                                                                        >  3:49303‑M2/135‑148 (MQ=255)
            ttgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttcttct                                                                                                                                        >  7:127572‑M2/135‑148 (MQ=255)
             tgatcttgtcgatccggacgacgacgagcgactgcgcggacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctccctCTCcttcttcttctt                                                                                                                                       >  4:169762‑M2/134‑148 (MQ=255)
                  ttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAg                                                                                                                                  <  4:100477‑M2/20‑1 (MQ=255)
                    gtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGa                                                                                                                                    <  6:85998‑M2/18‑1 (MQ=255)
                    gtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGa                                                                                                                                    <  5:56772‑M2/18‑1 (MQ=255)
                    gtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAg                                                                                                                                  <  4:128507‑M2/20‑1 (MQ=255)
                    gtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAg                                                                                                                                  <  2:106535‑M2/20‑1 (MQ=255)
                    gtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcagcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGa                                                                                                                                 <  8:318687‑M2/21‑1 (MQ=255)
                         accggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGa                                                                                                                           <  3:158895‑M2/27‑1 (MQ=255)
                            ggcagacgccgagcgactgcgcgtacaagctgaacttggcccgagcctggcaacccgcccggaacttgaagacagcgaagcctcggcacgttgctttgaagacgacgccggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGa                                                                                                                        <  5:219787‑M2/30‑1 (MQ=255)
                                 acgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACg                                                                                                                   <  2:13218‑M2/35‑1 (MQ=255)
                                 acgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACg                                                                                                                   <  4:49303‑M2/35‑1 (MQ=255)
                                       agcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTa                                                                                                             <  6:279802‑M2/41‑1 (MQ=255)
                                       agcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTa                                                                                                             <  1:278508‑M2/41‑1 (MQ=255)
                                        gcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCg                                                                                                                <  3:118882‑M2/38‑1 (MQ=255)
                                        gcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCg                                                                                                                >  4:118882‑M2/107‑144 (MQ=255)
                                          gactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACt                                                                                                          <  3:295621‑M2/44‑1 (MQ=255)
                                           actgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTc                                                                                                         >  1:55616‑M2/104‑148 (MQ=255)
                                            ctgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTcc                                                                                                        <  4:43364‑M2/46‑1 (MQ=255)
                                             tgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcagCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCa                                                                                                       >  2:233843‑M2/102‑148 (MQ=255)
                                                 cgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACtgt                                                                                                   <  6:275668‑M2/51‑1 (MQ=255)
                                                       atctgaacttggccagagcctggcaccccgcccggaactggaggacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCAtt                                                                                             >  7:213837‑M2/92‑148 (MQ=255)
                                                        tctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCAtt                                                                                             >  7:50749‑M2/91‑147 (MQ=255)
                                                         ctgaacttggccagagcctggcaacccgcccggaactggaagccatggcagcctcggcccgggtctttgaagacgacgacggcctgcctCTCCTTCTTCTTCTTTGAAGAGGCGGAAGAGCACGCCGGTAACTCCCCGGTGGCATTTa                                                                                           >  6:108906‑M2/90‑148 (MQ=255)
                                                            aacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCa                                                                                        >  1:279029‑M2/87‑148 (MQ=255)
                                                                tggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggtacgtttctttgaagacgacgacggcctgcatct                                                                                                                                                    >  7:60883‑M2/83‑84 (MQ=255)
                                                                tggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggtacgtttctttgaagacgacgacggcctgcatct                                                                                                                                                    <  8:60883‑M2/2‑1 (MQ=255)
                                                                      gagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACt                                                                                                          >  8:287770‑M2/77‑120 (MQ=255)
                                                                      gagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACt                                                                                                          <  7:287770‑M2/44‑1 (MQ=255)
                                                                           tggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGgtcgtc                                                                         >  1:301393‑M2/72‑148 (MQ=255)
                                                                            ggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTg                                                                                  <  8:102901‑M2/68‑1 (MQ=255)
                                                                            ggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTg                                                                                  >  7:102901‑M2/71‑138 (MQ=255)
                                                                            ggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCt                                                                        >  8:132573‑M2/71‑148 (MQ=255)
                                                                             gcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGgtcgtc                                                                         <  5:81177‑M2/77‑1 (MQ=255)
                                                                             gcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTg                                                                       <  7:137227‑M2/79‑1 (MQ=255)
                                                                             gcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTg                                                                       >  5:29901‑M2/70‑148 (MQ=255)
                                                                              caacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGt                                                                      <  1:89627‑M2/80‑1 (MQ=255)
                                                                              caacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGt                                                                      <  1:56088‑M2/80‑1 (MQ=255)
                                                                                acccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGttt                                                                    <  3:128508‑M2/82‑1 (MQ=255)
                                                                                 cccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgccgacggcctgcatATCCTTCTTCATCTTTGAAGATGCGGAAGATCACGCCGGTATATCCCCTGTGGCTTTTACCATCCGTGATGGTCGTCTGttt                                                                    >  8:209781‑M2/66‑147 (MQ=255)
                                                                                  ccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTAc                                                                  <  1:207972‑M2/84‑1 (MQ=255)
                                                                                    gcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTActc                                                                >  6:50976‑M2/63‑148 (MQ=255)
                                                                                       cggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgc                                                             <  3:70066‑M2/89‑1 (MQ=255)
                                                                                            ctggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgag                                                        <  7:91401‑M2/94‑1 (MQ=255)
                                                                                                 agacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgagcgtga                                                   <  2:55616‑M2/99‑1 (MQ=255)
                                                                                                 agacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgggagcgtga                                                   <  5:219591‑M2/99‑1 (MQ=255)
                                                                                                   acatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAAc                                                 <  6:29901‑M2/101‑1 (MQ=255)
                                                                                                      tcgaagcatcggcacgtttctttgaagacgacgccggcctgcatCTCCTTCTTGTTCTTTGAAGATGAGGAAGATCACGCCGGTAACTCAACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCTTGAGAGTGACCTGc                                              >  4:173229‑M2/45‑148 (MQ=255)
                                                                                                      tcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGc                                              >  7:24422‑M2/45‑148 (MQ=255)
                                                                                                        gaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGcc                                             <  2:279029‑M2/105‑1 (MQ=255)
                                                                                                          agcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGc                                          >  4:107906‑M2/41‑148 (MQ=255)
                                                                                                            catcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCAGCtt                                        >  5:164830‑M2/39‑148 (MQ=255)
                                                                                                              tcggcacgtttctttgaagacgacgacggcctgcatGTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCCTTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCtttt                                      <  3:169762‑M2/112‑1 (MQ=255)
                                                                                                              tcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCtttt                                      <  1:233843‑M2/112‑1 (MQ=255)
                                                                                                                      tttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTa                              <  6:6292‑M2/120‑1 (MQ=255)
                                                                                                                       ttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgagcgtga                                                   <  4:29260‑M2/99‑1 (MQ=255)
                                                                                                                       ttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgagcgtga                                                   >  3:29260‑M2/28‑126 (MQ=255)
                                                                                                                        tctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTAt                             <  2:301393‑M2/121‑1 (MQ=255)
                                                                                                                        tctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTAt                             <  5:50976‑M2/121‑1 (MQ=255)
                                                                                                                          tttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCg                           >  4:54569‑M2/25‑147 (MQ=255)
                                                                                                                           ttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGttt                                                                    <  1:37358‑M2/82‑1 (MQ=255)
                                                                                                                           ttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGttt                                                                    >  2:37358‑M2/24‑105 (MQ=255)
                                                                                                                              aagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGc                      >  3:164091‑M2/21‑148 (MQ=255)
                                                                                                                                    acgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgagcgtga                                                   >  2:211512‑M2/15‑113 (MQ=255)
                                                                                                                                    acgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgagcgtga                                                   <  1:211512‑M2/99‑1 (MQ=255)
                                                                                                                                        cggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAg            <  8:50749‑M2/138‑1 (MQ=255)
                                                                                                                                            ctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAg        <  7:132573‑M2/142‑1 (MQ=255)
                                                                                                                                              gcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTc      >  4:292743‑M2/5‑148 (MQ=255)
                                                                                                                                               catCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTc      <  3:54569‑M2/144‑1 (MQ=255)
                                                                                                                                                  cTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGATg  >  4:271561/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGATG  >  NC_000913/4001504‑4001797

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