Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 499,082 Δ52,248 bp [aes][ybcF] 46 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 499082 551329 52248 45 [0] [1] 46 [aes]–[ybcF] [aes],gsk,ybaL,fsr,ushA,ybaK,chiX,ybaP,ybaQ,copA,copA,copA,glsA,ybaT,cueR,ybbJ,qmcA,fetA,fetB,cnoX,ybbO,tesA,ybbA,ybbP,rhsD,ybbC,ylbG,selU,allS,allA,allR,allZ,gcl,hyi,glxR,ybbW,allB,ybbY,glxK,allE,allC,allD,fdrA,ylbE,ylbF,[ybcF]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 4990810 (0.000)45 (0.740) 38/480 0.6 100% coding (893/960 nt) aes acetylesterase
?NC_000913 551330 = 0 (0.000)coding (892/894 nt) ybcF putative carbamate kinase

CTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACGCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGC  >  NC_000913/551099‑551578
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        
gcccgggctttgctgcggattgtctggagacgctggaagagattgccgagcaaaaccgtgaggtcttcctcggtgccggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcggcggcggttttcatcatccgtgTAGTCGTAGGCATTAGACAt                                                                                                                                                                                                                                       <  2:293354‑M2/20‑1 (MQ=255)
  ccgggctttgctgcggattgtctggagacgctggaagagattgccgagcaaaaccgtgaggtcttcctcggtgccggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcggcggcggttttcatcatccgtgTAGTCGTAGGCATTAGACAttt                                                                                                                                                                                                                                     <  2:119050‑M2/22‑1 (MQ=255)
     ggctttgctgcggattgtctggagacgctggaagagattgccgagcaaaaccgtgaggtcttcctcggtgccggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaaccttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcggcggcggttttcatcatccgtgTAGTCGTAGGCATTAGACATTtgtg                                                                                                                                                                                                                                  <  1:204170‑M2/25‑1 (MQ=255)
         ttgctgcggattgtctggagacgctggaagagattgccgagcaaaaccgtgaggtcttcctcggtgccggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgtcgtcgcgaagagcctcggcggcggttttcatcatccgtgTAGTCGTAGGCATTAGATATTTGTGCCt                                                                                                                                                                                                                               <  2:473578‑M2/28‑1 (MQ=255)
                         ggagacgctggaagagattgccgagcaaaaccgtgaggtcttcctcggtgccggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcggcggcggttttcatcatccgtgTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGAc                                                                                                                                                                                                                <  2:250504‑M2/43‑1 (MQ=255)
                         ggagacgctggaagagattgccgagcaaaaccgtgaggtcttcctcggtgccggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcggcggcggttttcatcatccgtgTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGAc                                                                                                                                                                                                                <  1:386013‑M2/43‑1 (MQ=255)
                                       agattgccgagcaaaaccgtgaggtcttcctcggtgccggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaaagatggctaatcttgttgccgcggatcgctaaagctgagcggtaaagaaccgagcgccgccgcgaagagcctcggcggcggttttcatcatccctgTAGTCGTTGGCATTAGACATTTTTGCCCGATGCGACGCCTGACGCCTCTTTTCAAGCCt                                                                                                                                                                                                >  2:197910‑M2/193‑251 (MQ=255)
                                               nagcaaaaccgtgaggtcttcctcggtgccggcgggaaaaaatatgaatatattccggcgcctaatgccacgccgggaaatatcgcaatgatggctaatcttgttgccgcgtctcgccaaagctgagcggtaaagaaatgagcgccgtcgcgaagagcctcggcggcgggttttatcatccgggTAGTCGTAGGCATTAGAAATTTGTGCCCTATGCGCCGCTTGACGCGTCCTATCAAGCCTACAACCgg                                                                                                                                                                                        >  2:271118‑M2/185‑251 (MQ=255)
                                                   aaaatcgtgaggtcttcctcggtgccggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcggcggcggttttcatcatccgtgTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGcc                                                                                                                                                                                    <  2:566314‑M2/71‑1 (MQ=255)
                                                     aaccgtgaggtcttcctcggtgccggcgggaaaaaatatgaatatattcctgcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaaactgggcgggaaagaactgagcgccgtcgcgaagagcctcggcggcgggtttcttcaaccgtgTAGTCGTAGGCATTTGACATTTGTGCCTTATGCGAAGCTTGACGCGTCTTATCAGGCCTACAACCGGGGCCGc                                                                                                                                                                                  >  1:589579‑M2/179‑251 (MQ=255)
                                                       ccgtgaggtcttcctcggtgccggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcggcggcggttttcatcatccgtgTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCAt                                                                                                                                                                                <  1:372193‑M2/75‑1 (MQ=255)
                                                         gtgaggtcttcctcggtgccggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcggcggcggttttcatcatccgtgTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATcc                                                                                                                                                                              <  2:576547‑M2/77‑1 (MQ=255)
                                                         gtgaggtcttcctcggtgccggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcggcggcggttttcagcatccgtgTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGTTGCCGCATcc                                                                                                                                                                              <  2:470692‑M2/77‑1 (MQ=255)
                                                                             ggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcggcggcggttttcatcatccgtgTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGcc                                                                                                                                                          >  1:425252‑M2/155‑251 (MQ=255)
                                                                                  gaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcggcggcggttttcatcatccgtgTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTgatg                                                                                                                                                     >  1:157170‑M2/150‑251 (MQ=255)
                                                                                     aaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcggcggcggttttcatcatccgtgTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCtgatgcgn                                                                                                                                                  <  2:269364‑M2/105‑2 (MQ=255)
                                                                                      aaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcggcggcggttttcatcatccgtgTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTgatgcgat                                                                                                                                                 <  1:59573‑M2/106‑1 (MQ=255)
                                                                                          atgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcggcggcggttttcatcatccgtgTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCtt                                                                                                                                             <  2:430515‑M2/110‑1 (MQ=255)
                                                                                          atgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcggcggcggttttcatcatccgtgTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCCGATGCGATGCtt                                                                                                                                             <  1:426835‑M2/110‑1 (MQ=255)
                                                                                                atattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcggcggcggttttcatcatccgtgTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGc                                                                                                                                        >  1:411275‑M2/136‑250 (MQ=255)
                                                                                                atattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcggcggcggttttcatcatccgtgTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGc                                                                                                                                        <  2:522891‑M2/115‑1 (MQ=255)
                                                                                                 atattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaagaactgagcgccgtcgcgaagagcctcggcggcggttttcatcatccgtgTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCAt                                                                                                                                      >  2:64921‑M2/135‑251 (MQ=255)
                                                                                                        gcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcggcggcggttttcatcatccgtgTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATc                                                                                                                                >  2:121896‑M2/128‑250 (MQ=255)
                                                                                                                        cggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcggcggcggttttcatcatccgtgTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTAc                                                                                                               >  1:286749‑M2/112‑251 (MQ=255)
                                                                                                                           aacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcggcggcggttttcatcatccgtgTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCgg                                                                                                            >  1:227687‑M2/109‑251 (MQ=255)
                                                                                                                                atcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcggcggcggttttcatcatccgtgTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCgg                                                                                                       >  1:146272‑M2/104‑251 (MQ=255)
                                                                                                                                   gaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcggcggcggttttcatcatccgtgTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGCTGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGt                                                                                                      >  2:474332‑M2/101‑249 (MQ=255)
                                                                                                                                    aaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcggcggcggttttcatcatccgtgTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAgg                                                                                                   >  1:38876‑M2/100‑251 (MQ=255)
                                                                                                                                         atggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcggcggcggttttcatcatccgtgTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGa                                                                                              <  1:271118‑M2/157‑1 (MQ=255)
                                                                                                                                              taatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcggcggcggttttcatcatccgtgTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTCGGCCGGATAAgg                                                                                         >  2:190849‑M2/90‑251 (MQ=255)
                                                                                                                                                  cttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcggcggcggttttcatcatccgtgTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGtt                                                                                     >  1:478853‑M2/86‑251 (MQ=255)
                                                                                                                                                         ccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcggcggcggttttcatcatccgtgTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTAc                                                                                  <  1:370032‑M2/169‑1 (MQ=255)
                                                                                                                                                         ccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcggcggcggttttcatcatccgtgTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGGAGGCCGGATAAGGCGTTTAc                                                                                  >  2:370032‑M2/79‑247 (MQ=255)
                                                                                                                                                            cgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcggcggcggttttcatcatccgtgTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCAt                                                                           >  1:48190‑M2/76‑251 (MQ=255)
                                                                                                                                                             gtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcggcggcggttttcatcatccgtgTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATc                                                                          >  1:314775‑M2/75‑251 (MQ=255)
                                                                                                                                                             gtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcggcggcggttttcatcatccgtgTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATc                                                                          >  2:552521‑M2/75‑251 (MQ=255)
                                                                                                                                                                                  aaagaactgagcgccgtcgcgaagagcctcggcggcggttttcatcatccgtgTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAg                                                      <  1:47991‑M2/197‑1 (MQ=255)
                                                                                                                                                                                   aagaactgagcgccgtcgggaagagcctcggcggcggttttcatcttccgtgTAGTCGTAGGCATTAGACATTGGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGtt                                                    <  2:425252‑M2/199‑1 (MQ=255)
                                                                                                                                                                                     gaactgagcgccgtcgcgaagagcctcggcggcggttttcatcatccgtgTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTa                                                   >  1:110543‑M2/51‑250 (MQ=255)
                                                                                                                                                                                      aactgagcgccgtcgcgaagagcctcggcggcggttttcatcatccgtgTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAAc                                                 <  1:458385‑M2/202‑1 (MQ=255)
                                                                                                                                                                                        ctgagcgccgtcgcgaagagcctcggcggcggttttcatcatccgtgTATTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCg                                               >  1:448325‑M2/48‑251 (MQ=255)
                                                                                                                                                                                                 gtcgcgaagagcctcggcggcggttttcatcatccgtgTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGt                                                     >  2:518636‑M2/39‑236 (MQ=255)
                                                                                                                                                                                                 gtcgcgaagagcctcggcggcggttttcatcatccgtgTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGt                                                     <  1:518636‑M2/198‑1 (MQ=255)
                                                                                                                                                                                                                          ttcatcatccgggTAGTCGTAGGCAGTCGACATTTGGGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGAGGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCGGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCgcccaaatcacgacgacggcat             >  2:150398‑M2/14‑241 (MQ=255)
                                                                                                                                                                                                                                    gtgTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCcggcgg   >  1:180768‑M2/4‑251 (MQ=255)
                                                                                                                                                                                                                                      gtaggcgtaggcATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGCACTTACTCTACGCCCAAATCACGCCGCTGGCATATTCCcgcggc  >  2:330310/1‑250 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        
CTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACGCCGCAGCAACGTGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAATAGAGCACCAGTTAACCGAACTTACTCTGCGCCCAAATCACGCCGCTGGCATATTCCGGCGGC  >  NC_000913/551099‑551578

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