Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 2,286,372 Δ22,969 bp yejO[alkB] 24 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 2286372 2309340 22969 65 [0] [0] 65 yejO–[alkB] yejO,insH8,yejO,narP,ccmH,ccmG,ccmF,ccmE,ccmD,ccmC,ccmB,ccmA,napC,napB,napH,napG,napA,napD,napF,yojO,eco,mqo,yojI,[alkB]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 22863710 (0.000)65 (0.910) 60/490 0.1 100% intergenic (+84/+19) proL/yejO tRNA‑Pro/adhesin‑like autotransporter YejO
?NC_000913 2309341 = 0 (0.000)coding (1/651 nt)
coding (1065/1065 nt)
alkB
ada
DNA oxidative demethylase
DNA‑binding transcriptional dual regulator/DNA repair protein Ada

GCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGCTCTTGCCACGGTTCAGCATCGGCAAACAGATCCAACATTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGATGACAGGGTATTATGATAGCCAGCTTGTTGGCGGCACAGGCGCTGGCAACGGCCCGTACCGCTTTCGGTTTGCCGATGGCGTTAGCCAGTTGCTGATAACTGACGGTTTCACCGCAAGGTATCGTGCGCAGTGCCTGCCAGACT  >  NC_000913/2309104‑2309584
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                     
gattaattataaatcagttagcgaaatatcttacttgcaatcggtgtggaaaacggtagtattagcagccacgagtcggcacgtagcgcagcctggtagcgcaccgtcatggggtgtcgggggtcggaggttcaaatcctctcgtgccgaccaaaaatcccaagaaaaaaccaacccttacggttggtttttttatatctgcaattaattcgataaacagaccgtgacacatcacatTTACCTCTCCTCAt                                                                                                                                                                                                                                         >  1:344786‑M2/238‑251 (MQ=255)
 attaattataaatcagttagcgaaatatcttacttgcaatcggtgtggaaaacggtagtattagcagccacgagtcggcacgtagcgcagcctggtagcgcaccgtcatggggtgtcgggggtcggaggttcaaatcctctcgtgccgaccaaaaatcccaagaaaaaaccaacccttacggttggtttttttatatctgcaattaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCAtt                                                                                                                                                                                                                                        >  1:486171‑M2/237‑251 (MQ=255)
       tataaatcagttagcgaaatatcttacttgcaatcggtgtggaaaacggtagtattagcagccacgagtcggcacgtagcgcagcctggtagcgcaccgtcatggggtgtcgggggtcggaggttcaaatcctctcgtgccgaccaaaaatcccaagaaaaaaccaacccttacggttggtttttttatatctgcaattaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGc                                                                                                                                                                                                                                  >  2:171911‑M2/231‑251 (MQ=255)
        ataaatcagttagcgaaatatcttacttgcaatcggtgtggaaaacggtagtattagcagccacgagtcggcacgtagcgcagcctggtagcgcaccgtcatggggtgtcgggggtcggaggttcaaatcctctcgtgccgaccaaaaatcccaagaaaaaaccaacccttacggttggtttttttatatctgcaattaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCt                                                                                                                                                                                                                                 <  1:148760‑M2/22‑1 (MQ=255)
              cagttagcgaaatatcttacttgcaatcggtgtggaaaacggtagtattagcagccacgagtcggcacgtagcgcagcctggtagcgcaccgtcatggggtgtcgggggtcggaggttcaaatcctctcgtgccgaccaaaaatcccaagaaaaaaccaacccttacggttggtttttttatatctgcaattaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCgcgg                                                                                                                                                                                                                           <  2:161505‑M2/28‑1 (MQ=255)
                gttagcgaaatatcttacttgcaatcggtgtggaaaacggtagtattagcagccacgagtcggcacgtagcgcagcctggtagcgcaccgtcatggggtgtcgggggtcggaggttcaaatcctctcgtgccgaccaaaaatcccaagaaaaaaccaacccttacggttggtttttttatatctgcaattaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCgcggcg                                                                                                                                                                                                                         >  1:407634‑M2/222‑251 (MQ=255)
                   agcgaaatatcttacttgcaatcggtgtggaaaacggtagtattagaagccacgagtcggcacgtagcgcagcctggtagcgaaccgtcatggggtgtcgggggtcggaggttcaaatcctctcgtgccgaccaaaaatcccaagcaaaaaccaacccttacggttggtttttttatatctgcaattaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGgcgc                                                                                                                                                                                                                        >  1:551059‑M2/219‑249 (MQ=255)
                    gcgaaatatcttacttgcaatcggtgtggaaaacggtagtattagcagccacgagtcggcacgtagcgcagcctggtagcgcaccgtcatggggtgtcgggggtcggaggttcaaatcctctcgtgccgaccaaaaatcccaagaaaaaaccaacccttacggttggtttttttatatctgcaattaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCgcagc                                                                                                                                                                                                                     <  1:325807‑M2/34‑1 (MQ=255)
                                       atcggtgtggaaaacggtagtattagcagccacgagtcggcacgtagcgcagcctggtagcgcaccgtcatggggtgtcgggggtcggaggttcaaatcctctcgtgccgaccaaaaatcccaagaaaaaaccaacccttacggttggtttttttatatctgcaattaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGac                                                                                                                                                                                                   >  2:132642‑M2/199‑250 (MQ=255)
                                             gtggaaaacggtagtattagcagccacgagtcggcacgtagcgcagcctggtagcgcaccgtcatggggtgtcgggggtcggaggttcaaatcctctcgtgccgaccaaaaatcccaagaaaaaaccaacccttacggttggtttttttatatctgcaattaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGccc                                                                                                                                                                                             <  1:170842‑M2/58‑1 (MQ=255)
                                                     cggtagtattagcagccacgagtcggcacgtagcgcagcctggtagcgcaccgtcatggggtgtcgggggtcggaggttcaaatcctctcgtgccgaccaaaaatcccaagaaaaaaccaacccttacggttggtttttttatatctgcaattaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTa                                                                                                                                                                                     >  1:1347‑M2/185‑250 (MQ=255)
                                                      ggtagtattagcagccacgagtcggcacgtagcgcagcctggtagcgcaccgtcatggggtgtcgggggtcggaggttcaaatcctctcgtgccgaccaaaaatcccaagaaaaaaccaacccttacggttggtttttttatatctgcaattaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAAc                                                                                                                                                                                   <  2:196974‑M2/68‑1 (MQ=255)
                                                       gtagtattagcagccacgagtcggcacgtagcgcagcctggtagcgcaccgtcatggggtgtcgggggtcggaggttcaaatcctctcgtgccgaccaaaaatcccaagaaaaaaccaacccttacggttggtttttttatatctgcaattaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAAcc                                                                                                                                                                                  >  1:159185‑M2/183‑251 (MQ=255)
                                                                          gtcggcacgtagcgcagcctggtagcgcaccgtcatggggagtcgggggtcggcggttcaaatcctctcgtgccgaccaaaaatcccaagaaaaaaccaacccttacggttggtttttttatatctgcaattaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGATTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCCGCGGTAACCGGAAAGTGTGCCATcacc                                                                                                                                                                >  1:48646‑M2/164‑250 (MQ=255)
                                                                             ggcacgtagcgcagcctggaagcgcaccgtcatggggcgtcgggggtcggaggttcaaatcctctcgcgccgaccaaaactcccaagaaaaaaccaacccttacggttggtttttttatatctgcaattaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCCCCACgg                                                                                                                                                            >  2:299304‑M2/161‑251 (MQ=255)
                                                                              gcacgtagcgcagcctgggagcgcaccgtcatggggtgacgggggtcggaggttcaaatcctctcgtgccgaccaaaaatcccaagaacaaaccaacccttacggttggtttttttatatctgcaattaattcgataaacagacagtaacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACgg                                                                                                                                                            >  1:374891‑M2/160‑250 (MQ=255)
                                                                                  gtagcgcagcatggtagcgcaccgtcatggggtgtcgggggtcggaggttcaaatcctctcgtgccgacccaaaatcccaagaaaaaaccaacccttacggttggtttttttatatctgcaattaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGaccac                                                                                                                                                       >  2:220121‑M2/156‑251 (MQ=255)
                                                                                                    gcaccgtcatggggtgtcgggggtcggaggttcaaatcctctcgtgccgaccaaaaatcccaagaaaaaaccaacccttacggttggtttttttatatctgcaattaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGATGACAGGGtattat                                                                                                                                     >  1:156140‑M2/138‑251 (MQ=255)
                                                                                                     caccgtcatggggtgtcgggggtcggaggttcaaatcctctcgtgccgaccaaaaatcccaagaaaaaaccaacccttacggttggtttttttatatctgcaattaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGATGACAGGGTATTATg                                                                                                                                    >  2:6675‑M2/137‑251 (MQ=255)
                                                                                                                ggtgtcgggggtcggaggttcaaatcctctcgtgccgaccaaaaatcccaagaaaaaaccaacccttacggttggtttttttatatctgcaattaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGATGACAGGGTATTATGATAGCCCGCtt                                                                                                                         >  2:179836‑M2/126‑251 (MQ=255)
                                                                                                                   gtcgggggtcggaggttcaaatcctctcgtgccgaccaaaaatcccaagaaaaaaccaacccttacggttggtttttttatatctgcaattaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCATTTGCGCTTTACGCGGCACTCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACTTACCACCCGATGACAGGGTATTATGATAGCCAGCttgtt                                                                                                                      <  1:254908‑M2/129‑1 (MQ=255)
                                                                                                                            cggaggttcaaatcctctcgtgccgaccaaaaatcccaagaaaaaaccaacccttacggttggtttttttatatctgcaattaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGATGACAGGGTATTATGATAGCCAGCTTGTTGGCGGcaca                                                                                                             <  2:168174‑M2/138‑1 (MQ=255)
                                                                                                                                           tctcgtgccgaccaaatctctcaacaaaacaccaacccttacgcttgcattttttatatctgcaattaattcgatacacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCCGTTGCGCTTTACGCGACACCCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGATGACAGGGTATTATGATAGCCAGCTTGATGGCGGACCAGGcgctggccaaggc                                                                                              >  2:308220‑M2/99‑245 (MQ=255)
                                                                                                                                                   cgaccaaaaatcccaagaaaaaaccaacccttacggttggtttttttatatctgcaattaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGATGACAGGGTATTATGATAGCCAGCTTGTTGGCGGCACAGGCGCTGGCAACGGCCCGTACCg                                                                                      >  2:341001‑M2/91‑251 (MQ=255)
                                                                                                                                                      ccaaaaatcccaagaaaaaaccaacccttacggttggtttttttatatctgcaattaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGATGACAGGGTATTATGATAGCCAGCTTGTTGGCGGCACAGGCGCTGGCAACGGCCCGTACCGCtt                                                                                   <  2:486171‑M2/164‑1 (MQ=255)
                                                                                                                                                           aatcccaagaaaaaaccaacccttacggttggtttttttatatctgcaattaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGATGAAAGGGTATTATGATAGCCAGCTTGTTGGCGGCACAGGCGCTGGCAACGGCCCGTACCGCTTTCGGt                                                                              >  1:314103‑M2/83‑251 (MQ=255)
                                                                                                                                                           aatcccaagaaaaaaccaacccttacggttggtttttttatatctgcaattaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGAGTGCCATCCCCACGGACCACCCGATGACAGGGTATTATGATAGCCAGCTTGTTGGCGGCACAGGCGCTGGCAACGGCCCGGACCGCTTTCGGt                                                                              >  2:27472‑M2/83‑251 (MQ=255)
                                                                                                                                                                caagaaaaaaccaaccctttcgtttggtttttttatagctgcaattaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTGCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGATGACAGGGTATTATGATAGCCGGCATGTTGGCGGCACAGGCGCTGGCAACGGCCCGTACCTCTTTCGGTTTGcc                                                                        <  2:596628‑M2/174‑1 (MQ=255)
                                                                                                                                                                   gaaaaaaccaacccttcaggttggtttttttatatctgcaattaattcgatacacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGCTGACAGGGTATTATGATAGCCAGCTTGTTGGCGGACCAGGCGCTACCAACGGCCAGTACCGCTTTCGGTTTGCCGa                                                                      >  2:501118‑M2/75‑250 (MQ=255)
                                                                                                                                                                   gaaaaaaccaacccttacggttggtttttttatatctgcaattaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGATGACAGGGTATTATGATAGCCAGCTTGTTGGCGGCACAGGCGCTGGCAACGGCCCGTACCGCTTTCGGTTTGCCGAt                                                                     <  1:341001‑M2/177‑1 (MQ=255)
                                                                                                                                                                     aaaaaccaacccttacggttggtttttttatatctgcaattaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGATGACAGGGTATTATGATAGCCAGCTTGTTGGCGGCACAGGCGCTGGCAACGGCCCGTACCGCTTTCGGTTTGCCGATgg                                                                   >  1:393624‑M2/73‑251 (MQ=255)
                                                                                                                                                                      aaaaccaacccttacggttggtttttttatatctgcaattaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGATGACAGGGTATTATGATAGCCAGCTTGTTGGCGGCACAGGCGCTGGCAACGGCCCGTACCGCTTTCGGTTTGCCGATGGc                                                                  >  1:40411‑M2/72‑251 (MQ=255)
                                                                                                                                                                         accaacccttacggttggtttttttatatctgcacttaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGATGACAGGGTATTATGATAGCCAGCTTGTTGGCGGCACAGGCGCTGGCAAAGGCCCGTAACGCCTTCGGTGTGGCcgatgcggt                                                                >  2:437887‑M2/69‑247 (MQ=255)
                                                                                                                                                                          ccaacccttacggttggtttttttatatctgcaattaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGATGACAGGGTATTATGATAGCCAGCTTGTTGGCGGCACAGGCGCTGGCAACGGCCCGTACCGCTTTCGGTTTGCCGATGGCGtt                                                               >  1:196137‑M2/68‑250 (MQ=255)
                                                                                                                                                                               ccttacggttggttttgttatatctgcaattaattcgataaacagaccgtgacacatcacagTGTCCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTCTCGCGACACGCCCCAGCGGTACCCGGAAAGTGGGCCATCACCACGGGCCACCCGATGAGAGGGTATTATGATAGCCAGCTTGGTGGCGGCAGAGGCGCTGGCAACGGCCCGTACCGCTTACGGTTTGCCGATGGCGTTAGCCAg                                                         <  2:356055‑M2/189‑1 (MQ=255)
                                                                                                                                                                                   acggttggtttttttatatctgcaattaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGATGACAGGGTATTATGATAGCCAGCTTGTTGGCGGCACAGGCGCTGGCAACGGCCCGTACCGCTTTCGGTTTGCCGATGGCGTTAGCCAGTTGc                                                     <  1:159979‑M2/193‑1 (MQ=255)
                                                                                                                                                                                   acggttggtttttttatatctgcaattaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGATGACAGGGTATTATGATAGCCAGCTTGTTGGCGGCACAGGCGCTGGCAACGGCCCGTACCGCTTTCGGTTTGCCGATGGCGTTAGCCAGTTGc                                                     >  2:224085‑M2/59‑251 (MQ=255)
                                                                                                                                                                                    cggttggtttttttatatctgcaattaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGATGACAGGGTATTATGATAGCCAGCTTGTTGGCGGCACAGGCGCTGGCAACGGCCCGTACCGCTTTCGGTTTGCCGATGGCGTTAGCCAGTTGCt                                                    >  1:63261‑M2/58‑251 (MQ=255)
                                                                                                                                                                                    cggttggtttttttatatctgcaataaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGATGACAGGGTATTATGATAGCCAGCTTGTTGGCGGCACAGGCGCTGGCAACGGCCCGTCCCGCTTTCGGTTTGCCGATGGCGTTAGCCAGTTgcg                                                    >  1:501568‑M2/58‑250 (MQ=255)
                                                                                                                                                                                      gttggtttttttatatctgcaattaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGATGACAGGGTATTATGATAGCCAGCTTGTTGGCGGCACAGGCGCTGGCAACGGCCCGTACCGCTTTCGGTTTGCCGATGGCGTTCGCCAGTTGCTg                                                   >  2:6855‑M2/56‑250 (MQ=255)
                                                                                                                                                                                         ggtttttttatatctgcaattaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGATGACAGGGTATTATGATAGCCAGCTTGTTGGCGGCACAGGCGCTGGCAACGGCCCGTACCGCTTTCGGTTTGCCGATGGCGTTAGCCAGTTGCTGATa                                                >  1:209812‑M2/53‑250 (MQ=255)
                                                                                                                                                                                            ttttttatatctgcaattaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGATGACAGGGTATTATGATAGCCAGCTTGTTGGCGGCACAGGCGCTGGCAACGGCCCGTACCGCTTTCGGTTTGCCGATGGCGTTAGCCAGTTGCTGATAACt                                             <  2:159185‑M2/201‑1 (MQ=255)
                                                                                                                                                                                                ttatatctgcaattaattcgataaacagaccgtgccacatcatcgTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTCGCGCTTTACGCGACACGCCCCTGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGATGACAGTGTATTATGATAGCCAGCTTGTTGGCGGCACAGGCGCTGGCAACGGCCCGTACCGTTTTCGGTTTGCCGATGGCGTTAGCCAGTTGCTGATACCTGACg                                         <  2:344786‑M2/205‑1 (MQ=255)
                                                                                                                                                                                                ttatatctgcaattaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGATGACAGGGTATTATGATAGCCAGCTTGTTGGCGGCACAGGCGCTGGCAACGGCCCGTACCGCTTTCGGTTTGCCGATGGCGTTAGCCAGTTGCTGATAACTGACgg                                        >  2:431133‑M2/46‑251 (MQ=255)
                                                                                                                                                                                                 tatatctgcaattaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGATGACAGGGTATTATGATAGCCAGCTTGTTGGCGGCACAGGCGCTGGCAACGGCCCGTACCGCTTTCGGTTTGCCGATGGCGTTAGCCAGTTGCTGATAACTGACGGt                                       >  1:258260‑M2/45‑251 (MQ=255)
                                                                                                                                                                                                        gcaaataattcgataaacagaccgtgacacatcacagTTACCTCTCCGCATTGTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCCTCACCACGGACCACCCGATGACAGTGTATTAGGATAGCCAGCTTGTTGGCGGCACAGGCGCTGGCAACGGCCCGTACCGCTGTCGGTTTGCCGATTGCGTTAGCCAGTTTCGGATAACTGACGGTGTCAcc                                 <  2:551059‑M2/213‑1 (MQ=255)
                                                                                                                                                                                                         caattaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCCGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGCTGACAGGGTATTATGATAGCCAGCTTGTTGGCGGCACAGGCGCTGGCAACGGCCCGTACCGCTTTCGGTTTGCCGATGGCGTTAGCCAGTTGCTGATAACCGACGGTTTCACcgg                               >  2:389929‑M2/37‑250 (MQ=255)
                                                                                                                                                                                                         caattaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGATGACAGGGTATTATGATAGCCAGCTTGTTGGCGGCACAGGCGCTGGCAACGGCCCGTACCGCTTTCGGTTTGCCGATGGCGTTAGCCAGTTGCTGATAACTGACGGTTTCACCGc                               >  2:188816‑M2/37‑251 (MQ=255)
                                                                                                                                                                                                           attaattcgataaacaaaccgtgacacatcacatTTACATATCCTCATTTTCAGCTTCGCGGCGCAGCAGTTTCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGCTGACAGGGTATTATGATAGCCAGCTTGTTGGCGGCACAGGCGCTGGCAACGGCCCGTACCGCTTTCGGGTTGCCGATGGCGTTAGCCAGTTGCTGATAACTGCCGGATTCACCGCa                              >  1:447796‑M2/35‑250 (MQ=255)
                                                                                                                                                                                                            ttaattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGATGACAGGGTATTATGATAGCCCGCTTGTTGGCGGCACAGGCGCTGGCAACGGCCCGTACCGCTTTCGGTTTGCCGATGGCGTTAGCCCGTTGCTGATAACTGACGGTTTCACCGCa                              >  1:276237‑M2/34‑249 (MQ=255)
                                                                                                                                                                                                             taattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGATGACAGGGTATTATGATAGCCAGCTTGTTGGCGGCACAGGCGCTGGCAACGGCCCGTACCGCTTTCGGTTTGCCGATGGCGTTAGCCAGTTGCTGATAACTGACGGTTTCACCGCAAgg                           >  1:279684‑M2/33‑251 (MQ=255)
                                                                                                                                                                                                              aattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGATGACAGGGTATTATGATAGCCAGCTTGTTGGCGGCACAGGCGCTGGCAACGGCCCGTACCGCTTTCGGTTTGCCGATGGCGTTAGCCAGTTGCTGATAACTGACGGTTTCACCGCAAGGt                          >  1:51692‑M2/32‑251 (MQ=255)
                                                                                                                                                                                                              aattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGATGACAGGGTATTATGATAGCCAGCTTGTTGGCGGCACAGGCGCTGGCAACGGCCCGTACCGCTTTCGGTTTGCCGATGGCGTTAGCCAGTTGCTGATAACTGACGGTTTCACCGCAAGGt                          >  1:380670‑M2/32‑251 (MQ=255)
                                                                                                                                                                                                              aattcgataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGATGACAGGGTATTATGATAGCCAGCTTGTTGGCGGCACAGGCGCTGGCAAAGGCCCGTAACGCTTTTGGTTTGTCGCTGGAGTTAGCCAGTTTCTGataactgacgcttttacacccatgg                          >  2:19333‑M2/32‑240 (MQ=255)
                                                                                                                                                                                                                   gataaacagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGACAGTGTGCCATCACCACGGACCACCCGATGACAGGGTATTATGATAGCCAGCTTGTTGGCGGCACAGGCGCTGGCAACGGCCCGTACCGCTTTCGGTTTGCCGATGGAGTTAGCCAGTTGCTGATAACTGACGGTTTCAACGCAAGGTATCGt                     >  1:82798‑M2/27‑251 (MQ=255)
                                                                                                                                                                                                                   cataaaacagccgtgacaccacacagTTACCCTCCCCCATTTTCAGCCTCGCGGCGCCGCAGTTGCGTCTTACCCGACCCGGCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGATGACAGGGTATTATGATAGCCAGCTTGTTGGCGGCACAGGCGCTGGCAACGGCCCGTACCGCTTTCGGTTTGCCGATGGCGTTAGCCAGTTGCTGAGAACTGACGGTTTCACCGCAAGGTATCGt                     <  1:196661‑M2/225‑1 (MQ=255)
                                                                                                                                                                                                                        acagaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGATGACAGGGTATTATGATAGCCAGCTTGTTGGCGGCCCAGGCGCTGGCAACGGCCCGTACCGCTTTCGGTTTGCCGATGGCGTTAGCCAGTTGCTGATAACTGACGGTTTCACCGCAAGGTATCGTgcgc                 >  1:309127‑M2/22‑250 (MQ=255)
                                                                                                                                                                                                                          agaccgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGATGACAGGGTATTATGATAGCCAGCTTGTTGGCGGCACAGGCGCTGGCAACGGCCCGTACCGCTTTCGGTTTGCCGATGGCGTTAGCCAGTTGCTGATAACTGACGGTTTCACCGCAAGGTATCGTGCGCAGt              >  1:248791‑M2/20‑251 (MQ=255)
                                                                                                                                                                                                                            accgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGATGACAGGGTATTATGATAGCCAGCTTGTTGGCGGCACAGGCGCTGGCAACGGCCCGTACCGCTTTCGGTTTGCCGATGGCGTTAGCCAGTTGCTGATAACTGACGGTTTCACCGCAAGGGATCGTGCGCAGtgc            >  2:466557‑M2/18‑251 (MQ=255)
                                                                                                                                                                                                                              cgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGATGACAGGGTATTATGATAGCCAGCTTGTTGGCGGCACAGGCGCTGGCAACGGCCCGTACCGCTTTCGGTTTGCCGATGGCGTTAGCCAGTTGCTGATAACTGACGGTTTCACCGCAAGGTATCGTGCGCAGtgcct          <  1:9576‑M2/236‑1 (MQ=255)
                                                                                                                                                                                                                              cgtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAAACGGAAAGTGTGCCATCACCACGGACCACCCGATGACAGGGTATTATGATAGCCAGCTTGTTGGCGGCACAGGCGCTGGCAACGGCCCGTACCGCTTTCGGTTTGCCGATGGCGTTAGCCAGTTGCTGATAACTGACGGTTTCACCGCAAGGTATCGTGCGCAGtgcct          >  1:262697‑M2/16‑251 (MQ=255)
                                                                                                                                                                                                                               gtgacacatcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGATGACAGGGTATTATGATAGCCAGCTTGTTGGCGGCACAGGCGCTGGCAACGGCCCGTACCGCTTTCGGTTTGCCGATGGCGTTAGCCAGTTGCTGATAACTGACGGTTTCACCGCAAGGTATCGTGCGCAGtgcctg         >  1:430191‑M2/15‑251 (MQ=255)
                                                                                                                                                                                                                                    acatcacagTTGCCTCTCCTCCTTTTCAGCTTAGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGATGACAGGGTATTATGATAGCCAGCTTGTTGGCGGCACCGGCGCTGGCAACGGCCCGTACCGCTTTCGGTTTGCCGCTGGCGTTAGCCAGTTGCTGATAACTGACGGTTTCACCGCCAGGTATCGTGCGCAgtgcctgca       >  2:600769‑M2/10‑247 (MQ=255)
                                                                                                                                                                                                                                     catcaccgTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGATGACAGGGTATTATGATAGCCAGCTTGTTGGCGGCACAGGCGCTGGCAACGGCCCGTACCGCTTTCGGTTTGCCGATGGCGTTAGCCAGTTGCTGATAACTGCCGGTTTCACCGCAAGGGTTCGTGCGCAGTGCCTGCCAGAc   >  2:68446‑M2/9‑251 (MQ=255)
                                                                                                                                                                                                                                      atcacagTTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGATGACAGGGTATTATGATAGCCAGCTTGTTGGCGGCACAGGCGCTGGCAACGGCCCGTACCGCTTTCGGTTTGCCGATGGCGTTAGCCAGTTGCTGATAACTGACGGTTTCACCGCAAGGTATCGTGCGCAGTGCCTGCCAGACt  >  1:359125‑M2/8‑251 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                     
GCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGCTCTTGCCACGGTTCAGCATCGGCAAACAGATCCAACATTACCTCTCCTCATTTTCAGCTTCGCGGCGCAGCAGTTGCGCTTTACGCGACACGCCCCAGCGGTAACCGGAAAGTGTGCCATCACCACGGACCACCCGATGACAGGGTATTATGATAGCCAGCTTGTTGGCGGCACAGGCGCTGGCAACGGCCCGTACCGCTTTCGGTTTGCCGATGGCGTTAGCCAGTTGCTGATAACTGACGGTTTCACCGCAAGGTATCGTGCGCAGTGCCTGCCAGACT  >  NC_000913/2309104‑2309584

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