Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 565,054 44181 bp→TGA [intD]insL2 52 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 565054 609234–608008 42955–44181 70 [0] [1] 70 [intD]–insL2 [intD],renD,insE3,insF3,renD,emrE,ylcJ,ybcK,ybcL,ybcM,ylcH,ybcN,ninE,ybcO,rusA,ylcG,quuD,nmpC,insH2,nmpC,essD,rrrD,rzpD,rzoD,borD,ybcV,ybcW,ylcI,nohD,tfaD,tfaX,appY,ompT,envY,ybcH,nfrA,nfrB,cusS,cusR,cusC,cusF,cusB,cusA,pheP,ybdG,nfsB,ybdF,ybdJ,ybdK,sokE,hokE,insL2

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 5650530 (0.000)69 (0.900)
+TGA
62/488 0.1 99.3% coding (926/1164 nt) intD putative integrase
?NC_000913 609235 = 1 (0.010)noncoding (1229/1345 nt) IS186 repeat region

CAAATTGAACTGGCTTTTAAGCGGCTCAAAAGTTTGCTGCACCTGGATGCTTTGCGTGCAAAGGAACCTGAACTCGCGAAAGCGTGGATATTTGCTAATCTACTCGCCGCATTTTTAATTGACGACATAATCCAGCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGT  >  NC_000913/608992‑609479
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
gttacgccttctttatatcctccatcatttcagagtgggacatatttgggacattatcaccaaaaatgtcgtctattttcctcgcatgctctgtcaaatgattaggcgcaaggtgagcatacctacgaaccatttctatggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTAt                                                                                                                                                                                                                                               >  2:499977‑M2/244‑251 (MQ=255)
 ttacgccttctttatatcctccataatttcagagtgggacatatttgggacattatcaccaaaaatgtcgtctattttcctcgcatgctctgtcaaatgattaggcgcaaggtgagcatacctacgaaccatttctatggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATc                                                                                                                                                                                                                                              >  2:564324‑M2/243‑251 (MQ=255)
   ccgccttctttatatcctccataatttcagagtgggacatatttgggacattatcaccaaaaatgtcgtctattttcctcgcatgctctgtcaaatgattaggcgcaaggtgagcatacctacgaaccatttctatggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATccc                                                                                                                                                                                                                                            <  2:313585‑M2/11‑1 (MQ=255)
             tatatcctccataatttcagagtgggacatatttgggacattatcaccaaaaatgtcgtctattttcctcgcatgctctgtcaaatgattaggcgcaaggtgagcatacctacgaaccatttctatggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCt                                                                                                                                                                                                                                   >  1:380633‑M2/231‑250 (MQ=255)
              atatcccccataatttcagagtgggacatatttgggacattatcaccaaaaatgtcgtctattttcctcgcatgctctgtcaaatgattaggcgcaaggtgagcatacctacgaaccatttctatggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTAc                                                                                                                                                                                                                                 <  2:514515‑M2/22‑1 (MQ=255)
                atcctccataatttcagagtgggacatatttgggacattatcaccaaaaatgtcgtctattttcctcgcatgctctgtcaaatgattaggcgcaaggtgagcatacctacgaaccatttctatggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACa                                                                                                                                                                                                                                >  1:550249‑M2/228‑250 (MQ=255)
                atcctccataatttcagagtgggacatatttgggacattatcaccaaaaatgtcgtctattttcctcgcatgctctgtcaaatgattaggcgcaaggtgagcatacctacgaaccatttctatggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACa                                                                                                                                                                                                                                >  1:208958‑M2/228‑250 (MQ=255)
                  cctccagaaattcatagtggtacatatttgggacattatcaccaaaaatgtcgtctattttcctcgcatgctctgtcaaatgattaggcgcaaggtgatcatacctgcgaaccatttctatggactcccatccgcccatttccttaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACaaaa                                                                                                                                                                                                                             <  1:250449‑M2/26‑1 (MQ=255)
                  actccataatttcagagtgggacctatttggaacattatcaccaaaaatgcagtctattttcctcgcatgctctgtcaaatgattaggcgcaaggtgagcatacctacgaaccatttctatggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgccccggtgtgtctgaggccatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACaaa                                                                                                                                                                                                                              >  2:455628‑M2/226‑250 (MQ=255)
                    tccataatttcagagtgggacatatttgggacattatcaccaaaaatgtcgtctattttcctcgcatgctctgtcaaatgattaggcgcaaggtgagcatacctacgaaccatttctatggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACg                                                                                                                                                                                                                           >  1:287165‑M2/224‑251 (MQ=255)
                     ccataatttcagagtgggacatatttgggacattatcaccaaaaatgtcgtctattttcctcgcatgctctgtcaaatgattaggcgcaaggtgagcatacctacgaaccatttctatggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGc                                                                                                                                                                                                                          <  1:563532‑M2/29‑1 (MQ=255)
                         aatttcagagtgggacatatttgggacattatcaccaaaaatgtcgtctattttcctcgcatgctctgtcaaatgattaggcgcaaggtgagcatacctacgaaccatttctatggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAAt                                                                                                                                                                                                                       <  1:332684‑M2/32‑1 (MQ=255)
                            ttcagagtgggacatatttgggacattatcaccacaaatgtcgtctattttcctcgcatgctctgtcaaatgattaggcgcaaggtgcgcatacctacgaaccatttctatggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgcggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTga                                                                                                                                                                                                                    >  2:588144‑M2/216‑250 (MQ=255)
                                   tgggacatatttgggacattatcaccaaaaatgtcgtctattttcctcgcatgctctgtcaaatgattaggcgcaaggtgagcatacctacgaaccatttctatggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAAcgcg                                                                                                                                                                                                            >  2:89943‑M2/209‑251 (MQ=255)
                                      gacatatttgggacattatcaccaaaaatgtcgtctattttcctcgcatgctctgtcaaatgattaggcgcaaggtgagcatacctacgaaccatttctatggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGCCAGCCTCCAAAACGCAATTGAAGAACGCGAgg                                                                                                                                                                                                         >  1:386056‑M2/206‑251 (MQ=255)
                                      gacatatttgggacattatcaccaaaaatgtcgtctattttcctcgcatgctctgtcaaatgattaggcgcaaggtgagcatacctacgaaccatttctatggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAgg                                                                                                                                                                                                         <  1:136548‑M2/46‑1 (MQ=255)
                                             ttgggacattatcaccaaaaatgtcgtctattttcctcgcatgctctgtcaaatgattaggcgcaaggtgagcatacctacgaaccatttctatggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTc                                                                                                                                                                                                  >  1:339605‑M2/199‑251 (MQ=255)
                                               gggacattatcaccaaaaatgtcgtctattttcctcgcatgctctgtcaaatgattaggcgcaaggtgagcatacctacgatccatttctatggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatcgttgaGACAGGATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGGGGCATCGTCtt                                                                                                                                                                                                <  2:164352‑M2/55‑1 (MQ=255)
                                               gggacattatcaccaaaaatgtcgtctattttcctcgcatgctctgtcaaatgattaggcgcaaggtgagcatacctacgaaccatttctatggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCtt                                                                                                                                                                                                >  1:5479‑M2/197‑251 (MQ=255)
                                               gggacattatcaccaaaaatgtcgtctattttcctcgcatgctctgtcaaatgattaggcgcaaggtgagcatacctacgaaccatttctatggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCtt                                                                                                                                                                                                <  1:499977‑M2/55‑1 (MQ=255)
                                                 gacattatcaccaaaaatgtcgtctattttcctcgcatgctctgtcaaatgattaggcgcaaggtgagcatacctacgaaccatttctatggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTa                                                                                                                                                                                               >  2:231131‑M2/195‑250 (MQ=255)
                                                    attatcaccaaaaatgtcgtctattttcctcgcatgctctgtcaaatgattaggcgcaaggtgagcatacctacgaaccatttctatggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACg                                                                                                                                                                                            >  1:403433‑M2/192‑250 (MQ=255)
                                                     ttatcaccaaaaatgtcgtctattttcctcgcatgctctgtcaaatgattaggcgcaaggtgagcatacctacgaaccatttctatggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACg                                                                                                                                                                                            >  1:145496‑M2/191‑249 (MQ=255)
                                                                   gtcgtctattttcctcgcatgctctgtaaaatgattaggcgcaaggtgagcatacctacgaaccatttctatggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCg                                                                                                                                                                            >  2:109638‑M2/177‑251 (MQ=255)
                                                                   gtcgtctatttgcctcgcatgctctgtcaaatgattaggcgcaaggtgagcatacctacgaaccatttctatgggctcacatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaagctttgaGACATTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCg                                                                                                                                                                            <  2:257050‑M2/75‑1 (MQ=255)
                                                                     cgtctattttcctcgcatgctctgtcaaatgattaggcgcaaggtgagcatacctacgaaccatttctatggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGc                                                                                                                                                                           <  1:212720‑M2/76‑1 (MQ=255)
                                                                            tttcctcgcatgctctgtcaaatgattaggcgcaaggtgagcatacctacgaaccatttctatggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGGCGCAttcttc                                                                                                                                                                    >  2:11496‑M2/168‑250 (MQ=255)
                                                                               ccgcgcatgctctgtcaaatgattaggcgcaaggtgagcatacctacgaaccatttctatggactcccatccgcccatgtcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCTAGGCGTCGCATGCTTCAGAt                                                                                                                                                                <  2:548947‑M2/87‑1 (MQ=255)
                                                                                  cgcatgctctgtcaaatgattaggcgcaaggtgagcatacctacgaaccatttctatggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatgtaaccggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCTAGGCGTCGCATTCTTCAGATGGt                                                                                                                                                             <  2:613059‑M2/90‑1 (MQ=255)
                                                                                  cgcatgctctgtcaaatgattaggcgcaaggtgagcatacctacgaaccatttctatggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGt                                                                                                                                                             <  1:408463‑M2/90‑1 (MQ=255)
                                                                                     atgctctgtcaaatgattaggcgcaaggtgagcatacctacgaaccatttctatggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTc                                                                                                                                                           >  2:321585‑M2/159‑250 (MQ=255)
                                                                                     atgctctgtcaaatgattaggcgcaaggtgagcatacctacgaaccatttctatggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGAGTCGCATTCTTCAGATGGTTc                                                                                                                                                           >  2:179967‑M2/159‑250 (MQ=255)
                                                                                       gctctgtcaaatgattaggcgcaaggtgagcatacctacgaaccatttctatggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCAtt                                                                                                                                                                        <  1:260802‑M2/79‑1 (MQ=255)
                                                                                       gctctgtcaaatgattaggcgcaaggtgagcatacctacgaaccatttctatggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCAtt                                                                                                                                                                        >  2:260802‑M2/157‑235 (MQ=255)
                                                                                           tgtcaaatgattaggcgcaaggtgagcatacctacgaaccatttctatggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGAt                                                                                                                                                                <  1:325615‑M2/87‑1 (MQ=255)
                                                                                                 atgattaggcgcaaggtgagcatacctacgaaccatttctatggactcccatccgcccatttcctgaagcactgatactgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTAACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACAGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGtt                                                                                                                                               >  2:186520‑M2/147‑250 (MQ=255)
                                                                                                         gcgcaaggtgagcatacctacgaaccatttctatggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTAt                                                                                                                                      >  1:434200‑M2/139‑251 (MQ=255)
                                                                                                            caaggtgagcatacctacgaaccattgctatggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGGCAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATggg                                                                                                                                   <  2:71132‑M2/116‑1 (MQ=255)
                                                                                                                  gagcatacctacgaaccatttctatggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAAt                                                                                                                             >  2:381821‑M2/130‑251 (MQ=255)
                                                                                                                     catacctacgaaccatttctatggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATccc                                                                                                                          <  1:109638‑M2/125‑1 (MQ=255)
                                                                                                                      atacctacgaaccatttctatggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATcccc                                                                                                                         <  1:130058‑M2/126‑1 (MQ=255)
                                                                                                                        acctacgaaccatttctatggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGc                                                                                                                       >  2:629365‑M2/124‑251 (MQ=255)
                                                                                                                         cctacgatccttttctatggactcccatcctcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggttatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCa                                                                                                                      <  2:466595‑M2/129‑1 (MQ=255)
                                                                                                                               aaccatttctatggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCg                                                                                                                >  1:44682‑M2/117‑251 (MQ=255)
                                                                                                                                accatttctatggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCgg                                                                                                               <  1:53471‑M2/136‑1 (MQ=255)
                                                                                                                                 ccatttctatggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGt                                                                                                              <  2:208958‑M2/137‑1 (MQ=255)
                                                                                                                                         atggactcccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTgtcaggt                                                                                                      >  2:515538‑M2/107‑251 (MQ=255)
                                                                                                                                                cccatccgcccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctggggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATg                                                                                               >  2:193247‑M2/100‑251 (MQ=255)
                                                                                                                                                         ccatttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCAc                                                                                      <  1:514591‑M2/161‑1 (MQ=255)
                                                                                                                                                            tttcctgaagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTg                                                                                   >  1:35203‑M2/88‑251 (MQ=255)
                                                                                                                                                                   aagcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTaca                                                                            <  1:203562‑M2/171‑1 (MQ=255)
                                                                                                                                                                     gcactgataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACt                                                                          <  2:147576‑M2/173‑1 (MQ=255)
                                                                                                                                                                          gataatgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACTCTTTCGc                                                                     <  1:543538‑M2/178‑1 (MQ=255)
                                                                                                                                                                              atgggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGCGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGt                                                                 >  1:581889‑M2/70‑251 (MQ=255)
                                                                                                                                                                                gggacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTgg                                                               <  2:145496‑M2/184‑1 (MQ=255)
                                                                                                                                                                                  gacgcctgactgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCGACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTc                                                             <  2:259039‑M2/186‑1 (MQ=255)
                                                                                                                                                                                        tgactgaatcagccagcattcccaggtgcgactacggtcatggaaacggaaatctttgaGACAGTATCACTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGCGGCGTCGCATTCTTCAGATTGTTCAAACCTTAAGTTATCGCTTATGGGAGGAATCCGCGCATATCCGGATGCCAGGTCCGGCTGGTAAGGCACCTGCTTTACACTTTCGCTCGTGGTCAGTGAt                                                       >  2:509841‑M2/60‑251 (MQ=255)
                                                                                                                                                                                            tgaatcagccagcttgcccaggtgtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCt                                                   >  1:125044‑M2/56‑251 (MQ=255)
                                                                                                                                                                                                            cccaggtgagtctgaggacaggggaacggaaatctttgaTCCAGCATCCCTGACAGCCTACAACACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACAGAGTCGCCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTCCACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAAtcgtccc                                   >  2:523921‑M2/40‑248 (MQ=255)
                                                                                                                                                                                                                aTGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTAcc                                   >  1:381981/1‑247 (MQ=255)
                                                                                                                                                                                                                   gtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGAtgtt                            >  2:391757‑M2/33‑251 (MQ=255)
                                                                                                                                                                                                                   gtgtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGAtgt                             <  1:196011‑M2/218‑1 (MQ=255)
                                                                                                                                                                                                                     gtctgaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGAtgttg                           <  1:381821‑M2/220‑1 (MQ=255)
                                                                                                                                                                                                                         gaggtcatggaaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCa                        <  1:629365‑M2/223‑1 (MQ=255)
                                                                                                                                                                                                                              catggaaacggaaatctttgaGACAGGATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGGTGTGGGCGAATCGTACCAGATGTTGTCAATTAATc                 >  2:410841‑M2/22‑251 (MQ=255)
                                                                                                                                                                                                                                   aaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGtt            <  2:434200‑M2/235‑1 (MQ=255)
                                                                                                                                                                                                                                   aaacggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGtt            >  2:104329‑M2/17‑251 (MQ=255)
                                                                                                                                                                                                                                      cagaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGc         >  2:265185‑M2/14‑251 (MQ=255)
                                                                                                                                                                                                                                       ggaaatctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGc         >  2:143457‑M2/13‑250 (MQ=255)
                                                                                                                                                                                                                                             ctttgaGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGt  >  1:121220‑M2/7‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
CAAATTGAACTGGCTTTTAAGCGGCTCAAAAGTTTGCTGCACCTGGATGCTTTGCGTGCAAAGGAACCTGAACTCGCGAAAGCGTGGATATTTGCTAATCTACTCGCCGCATTTTTAATTGACGACATAATCCAGCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGT  >  NC_000913/608992‑609479

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