Predicted mutation
evidence seq id position mutation freq annotation gene description
MC JC NC_000913 4,188,959 Δ9 bp 100% coding (3611‑3619/4224 nt) rpoC → RNA polymerase, beta prime subunit

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 4188959 4188967 9 94 [0] [0] 94 rpoC RNA polymerase, beta prime subunit

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 41889580 (0.000)81 (1.110) 63/262 NT 100% coding (3610/4224 nt) rpoC RNA polymerase, beta prime subunit
?NC_000913 4188968 = 0 (0.000)coding (3620/4224 nt) rpoC RNA polymerase, beta prime subunit

GAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATC  >  NC_000913/4188820‑4189116
                                                                                                                                                    |                                                                                                                                                    
atcctggctgaaatcagcggtatcgtttccttcggtaaagaaaccaaaggtaaacgtcgtctggttatcaccccggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgg                                                                                                                                                      >  1:201071‑M2/149‑149 (MQ=255)
atcctggctgaaatcagcggtatcgtttccttcggtaaagaaaccaaaggtaaacgtcgtctggttatcaccccggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgg                                                                                                                                                      <  1:885888‑M2/1‑1 (MQ=255)
   ctggctgaaatcagcggtatcgtttccttcggtaaagaaaccaaaggtaaacgtcgtctggttatcaccccggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGa                                                                                                                                                   <  1:1022464‑M2/4‑1 (MQ=255)
   ctggctgaaatcagcggtatcgtttccttcggtaaagaaaccaaaggtaaacgtcgtctggttatcaccccggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGa                                                                                                                                                   <  1:432889‑M2/4‑1 (MQ=255)
   ctggctgaaatcagcggtatcgtttccttcggtaaagaaaccaaaggtaaacgtcgtctggttatcaccccggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGa                                                                                                                                                   >  1:858458‑M2/146‑149 (MQ=255)
    tggctgaaatcagcggtatcgtttccttcggtaaagaaaccaaaggtaaacgtcgtctggttatcaccccggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgatggaagctgaacgtgCTGAc                                                                                                                                                  >  2:1220227‑M2/145‑149 (MQ=255)
          aaatcagcggtatcgtttccttcggtaaagaaaccaaaggtaaacgtcgtctggttatcaccccggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAAtt                                                                                                                                            <  2:933869‑M2/11‑1 (MQ=255)
               agcggtatctgttccttcggtaaagaaaccaaaggtaagcgtcgtctggttatccccccggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagttcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGa                                                                                                                                       <  2:1229657‑M2/16‑1 (MQ=255)
                    gatcgtttccttctgtaaagaaaccaaaggtaaacgtcgtctggttatcaccccggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTc                                                                                                                                  <  2:1219331‑M2/21‑1 (MQ=255)
                          ttccctcggtaaagaaaccaaaggtaaacgtcgtctggttatcaccccggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcaactcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAg                                                                                                                            <  1:200160‑M2/27‑1 (MQ=255)
                           tccttcggtaaagaaaccaaaggtaaacgtcgtctggttatcaccccggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAgc                                                                                                                           <  2:1175624‑M2/28‑1 (MQ=255)
                           accttcggtaaagaacccaaaggtaagcgtcgtctggtgatgcccccgctagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAgc                                                                                                                           <  2:1312700‑M2/28‑1 (MQ=255)
                            ccttcggtaaagaaaccaaaggtcaacgtcgtctgaatatcaccccggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAgcg                                                                                                                          <  2:478550‑M2/29‑1 (MQ=255)
                             cttcggtaaagaaaccaaaggtaaacgtcgtcgggttatcaccccggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAgcgc                                                                                                                         <  2:297632‑M2/30‑1 (MQ=255)
                               tcggtaaagaaaccaaaggtaaacgtcgtctggttatcaccccggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGcgccg                                                                                                                       <  2:635381‑M2/32‑1 (MQ=255)
                                cggtaaagaaaccaaaggtaaacgtcgtctggttatcaccccggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGcgccgc                                                                                                                      <  1:1126465‑M2/33‑1 (MQ=255)
                                  gtaaagaaaccaaaggtaaacgtcgtctggttatcaccccggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcggaggtgaacgtgGGGACGTAACTTCCGACGGTCCGGACGCGCAGCAc                                                                                                                    >  2:135441‑M2/115‑149 (MQ=255)
                                   taaagaaaccaaaggtaaacgtcgtctggttatcaccccggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCaca                                                                                                                   >  1:277997‑M2/114‑148 (MQ=255)
                                        aaaccaaaggtaaacgtcgtctggttatcaccccggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACAtt                                                                                                              <  1:434659‑M2/41‑1 (MQ=255)
                                        aaaccaaaggtaaacgtcgtctggttatcaccccggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgcgGTGACGTACTTTCCGACGGTCCGGAAGCGCCGCACGACAtt                                                                                                              >  1:49624‑M2/109‑149 (MQ=255)
                                          accaaaggtaaacgtcgtctggttatcaccccggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATtct                                                                                                            >  1:536684‑M2/107‑149 (MQ=255)
                                           ccaaaggtaaacgtcgtctggttatcaccccggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATtctg                                                                                                           <  2:201071‑M2/44‑1 (MQ=255)
                                           ccaaaggtaaacgtcgtctggttatcaccccggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTCCCGACGGTCCGGAAGCGCCGCACGACATtctg                                                                                                           >  2:1205319‑M2/106‑149 (MQ=255)
                                            caaatgtaaacgtcgtctgtttatcaccccggtagacggtagcgatccctacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATtctgc                                                                                                          <  1:1008407‑M2/45‑1 (MQ=255)
                                            caaaggtaaacgtcgtctggttatcaccccggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATtctgc                                                                                                          >  2:1117444‑M2/105‑149 (MQ=255)
                                             aaaggtaaacgtcgtctggttatcaccccggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATtctgcg                                                                                                         <  2:270016‑M2/46‑1 (MQ=255)
                                             aaaggtaaacgtcgtctggttatcaccccggtagacggtagcgatccgaccgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATtctgcg                                                                                                         >  1:932638‑M2/104‑149 (MQ=255)
                                               aggtaaacgtcgtctggttatcaccccggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgggGTGACGTAATTTCCGACCGTGCGGAAGCGCCGACCGACATtctgcgtc                                                                                                       >  2:564756‑M2/102‑149 (MQ=255)
                                                 gtaaacgtcgtctggttatcaccccggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCGGACGGTCCGGCAGCGCTGCACGAGATtgtgcgtctg                                                                                                     >  1:1073686‑M2/100‑149 (MQ=255)
                                                    aacgtcgtctggttatcaccccggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTctgcgtctgcgt                                                                                                  <  1:355865‑M2/53‑1 (MQ=255)
                                                    aacgtcgtctggttatcaccccggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTctgcgtctgcgt                                                                                                  <  1:963517‑M2/53‑1 (MQ=255)
                                                            ctggttatcaccccggtagacggtagcgttccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTTATGTTCa                                                                                          <  2:568056‑M2/61‑1 (MQ=255)
                                                            ccggttatcacccaggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggggaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCa                                                                                          <  1:688127‑M2/61‑1 (MQ=255)
                                                               gttatcaccccggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagcgcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGc                                                                                       >  2:956451‑M2/86‑149 (MQ=255)
                                                                  atccccccggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgggttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGt                                                                                    <  1:153587‑M2/67‑1 (MQ=255)
                                                                  atcaccccggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGt                                                                                    <  1:1218038‑M2/67‑1 (MQ=255)
                                                                    caccccggaagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTTACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTa                                                                                  <  1:478572‑M2/69‑1 (MQ=255)
                                                                     accccggtagacggtagcgctctgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTAc                                                                                 >  2:1308391‑M2/80‑149 (MQ=255)
                                                                     accccggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTGATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGGGttcatgccggtgc                                                                                 >  2:360839‑M2/80‑143 (MQ=255)
                                                                     accccggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTAc                                                                                 >  2:839452‑M2/80‑149 (MQ=255)
                                                                     accccggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTAc                                                                                 >  1:1062745‑M2/80‑149 (MQ=255)
                                                                      ccccggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCATGACATTCTGCGTCGGCGTGGTGTTGATGCTGTTACt                                                                                >  2:551628‑M2/79‑149 (MQ=255)
                                                                      ccccggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACt                                                                                >  1:548515‑M2/79‑149 (MQ=255)
                                                                      ccccggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGACTGCGTGGTGTCCATGCTGTTACt                                                                                >  2:579684‑M2/79‑149 (MQ=255)
                                                                       cccggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTc                                                                               <  1:1239711‑M2/72‑1 (MQ=255)
                                                                       cccggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGGGGTGTTCATGCTGTTACTc                                                                               >  1:1164942‑M2/78‑149 (MQ=255)
                                                                        ccggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCg                                                                              <  1:832281‑M2/73‑1 (MQ=255)
                                                                         cggtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGt                                                                             <  1:1301044‑M2/74‑1 (MQ=255)
                                                                           gtagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTCCTCGTTa                                                                           >  1:838357‑M2/74‑149 (MQ=255)
                                                                            tagacggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTAc                                                                          >  2:662129‑M2/73‑149 (MQ=255)
                                                                                 ggtagcgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGATCCGGAAGCGCCGCACGACATTCTGCGTCTGCTTGGTGTTCATGCTGTTACTCGTTACATCGt                                                                     >  2:978512‑M2/68‑149 (MQ=255)
                                                                                     gcgattcgtcctacgagatgattccgcaatggggtcagctcaacgtgtgtgaagctgaacgtgGTGACGTAATTTCCGACGGTCCGGACTCGCCGCACGACATTATGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAAc                                                                 <  2:412632‑M2/86‑1 (MQ=255)
                                                                                      cgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGCAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACg                                                                <  2:976956‑M2/87‑1 (MQ=255)
                                                                                      cgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACg                                                                >  2:175812‑M2/63‑149 (MQ=255)
                                                                                      cgatccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggggaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGGTTCCATGCTGTTACTCGTTCCATGGCTAACg                                                                >  1:1158728‑M2/63‑149 (MQ=255)
                                                                                        atccgtacgaagagatgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGaa                                                              <  1:1250092‑M2/89‑1 (MQ=255)
                                                                                         tccgtacgaagagagtattccgaaagggcgtcagcccaacgtgttcgaaggggaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCCGTTACTCGTTACATCGTTAACGAAg                                                             <  2:686260‑M2/90‑1 (MQ=255)
                                                                                           cgtacgaagagctgattccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTa                                                           <  2:56648‑M2/92‑1 (MQ=255)
                                                                                                   gagatgattccgaaatggcgtcagctcaacgtgctcgaaggtgaacgtgATGACGAAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACGCGTTACATCGTTAACGACGTACCGGACGt                                                   >  2:864929‑M2/50‑149 (MQ=255)
                                                                                                   gagataattccggaatggcgtcagctcaacgtgttcggaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGt                                                   <  2:206957‑M2/100‑1 (MQ=255)
                                                                                                      atgattccgaaatggcgccagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGtata                                                <  1:1172285‑M2/103‑1 (MQ=255)
                                                                                                         attccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCg                                             >  2:1027591‑M2/44‑149 (MQ=255)
                                                                                                          ttccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGt                                            >  1:798766‑M2/43‑149 (MQ=255)
                                                                                                          ttccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGt                                            >  1:1078314‑M2/43‑149 (MQ=255)
                                                                                                          ttccgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGGCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGGGTTCATGCTGTTACTCGTTACATCGTTGACGAAGTACCGGACGTATGCCGt                                            >  2:538027‑M2/43‑149 (MQ=255)
                                                                                                             cgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTg                                         <  1:1078062‑M2/110‑1 (MQ=255)
                                                                                                             cgaaatggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACAGTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTg                                         <  1:1144285‑M2/110‑1 (MQ=255)
                                                                                                                  tggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAggg                                    >  1:611498‑M2/35‑149 (MQ=255)
                                                                                                                  tggcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGGCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAggg                                    <  2:76685‑M2/115‑1 (MQ=255)
                                                                                                                   ggcgtcagcccaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGc                                   <  1:202816‑M2/116‑1 (MQ=255)
                                                                                                                    gcgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCCGCAGcgcg                                  >  2:63248‑M2/33‑149 (MQ=255)
                                                                                                                    gcgtcagctcaacgtgttcgaaggggaacgtgGCGACGTACTTACAATCGGGCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCg                                  <  1:104930‑M2/117‑1 (MQ=255)
                                                                                                                     cgtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGt                                 >  1:9422‑M2/32‑149 (MQ=255)
                                                                                                                      gtcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGtt                                >  2:1343746‑M2/31‑149 (MQ=255)
                                                                                                                       tcagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTa                               >  1:179830‑M2/30‑149 (MQ=255)
                                                                                                                        cagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTaa                              <  1:1117444‑M2/121‑1 (MQ=255)
                                                                                                                        cagctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACACCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTaa                              <  1:729157‑M2/121‑1 (MQ=255)
                                                                                                                          gctcaacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGa                            <  1:1187198‑M2/123‑1 (MQ=255)
                                                                                                                              aacgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTaa                        <  2:1168134‑M2/127‑1 (MQ=255)
                                                                                                                                cgtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACg                      <  1:126309‑M2/129‑1 (MQ=255)
                                                                                                                                 gttgtcgaaggtgaccgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGa                     <  1:502602‑M2/130‑1 (MQ=255)
                                                                                                                                 gtgttcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGa                     <  1:662129‑M2/130‑1 (MQ=255)
                                                                                                                                     tcgaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATaaa                 <  2:242633‑M2/134‑1 (MQ=255)
                                                                                                                                     gcgaaggtgaacgtgGTGACGGAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATaaa                 <  2:49624‑M2/134‑1 (MQ=255)
                                                                                                                                       gagggtgaccgtgGTGACGTAATCTCCGACGGTCCGCACGCGCCGCCAGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAATGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAaca               <  1:94123‑M2/136‑1 (MQ=255)
                                                                                                                                       gaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATTGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAaacc               >  2:252493‑M2/14‑148 (MQ=255)
                                                                                                                                       gaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATCAaca               >  1:438720‑M2/14‑149 (MQ=255)
                                                                                                                                       gaaggtgaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAaca               >  2:1174980‑M2/14‑149 (MQ=255)
                                                                                                                                           gtgaacgtgGTGACGTAATTTCCGAAGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACAGTCTGCAGGGCGTTAAGATTAACGATAAACACATc           >  1:755102‑M2/10‑149 (MQ=255)
                                                                                                                                             gaacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTCCAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGa         >  2:640979‑M2/8‑149 (MQ=255)
                                                                                                                                              aacgtgGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGaa        <  1:929812‑M2/143‑1 (MQ=255)
                                                                                                                                              aacgggGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGaa        <  2:1255929‑M2/143‑1 (MQ=255)
                                                                                                                                                   gGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAcgttgt   >  2:771244‑M2/2‑147 (MQ=255)
                                                                                                                                                    gtgACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTAtc  <  2:179830/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGGTAGACGGTAGCGATCCGTACGAAGAGATGATTCCGAAATGGCGTCAGCTCAACGTGTTCGAAGGTGAACGTGTAGAACGTGGTGACGTAATTTCCGACGGTCCGGAAGCGCCGCACGACATTCTGCGTCTGCGTGGTGTTCATGCTGTTACTCGTTACATCGTTAACGAAGTACAGGACGTATACCGTCTGCAGGGCGTTAAGATTAACGATAAACACATCGAAGTTATC  >  NC_000913/4188820‑4189116

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