Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 256,108 Δ115,274 bp [gpt][mhpC] 117 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 256108 371381 115274 39 [0] [1] 40 [gpt]–[mhpC] [gpt],frsA,crl,insB9,insA9,crl,phoE,proB,proA,thrW,ykfI,yafW,ykfH,ykfG,yafX,ykfF,ykfB,yafY,yafZ,ykfA,perR,insN,insI1,insN,eyeA,ykfC,insH1,mmuP,mmuM,afuC,insB2,insA2,yagB,yagA,yagE,yagF,yagG,yagH,xynR,argF,argL,ykgS,insB3,insA3,yagJ,yagK,yagL,yagM,ykgV,yagN,intF,yagP,paoD,paoC,paoB,paoA,yagU,ykgJ,ecpE,ecpD,ecpC,ecpB,ecpA,matA,ykgL,ykgO,ykgM,ykgR,insE1,insF1,rclC,rclB,rclA,rclR,ykgE,ykgF,ykgG,ykgH,betA,betB,betI,betT,yahV,pdeL,yahB,yahC,yahD,yahE,yahF,yahG,yahH,yahI,yahJ,yahK,yahL,yahM,yahN,yahO,prpR,prpB,prpC,prpD,prpE,codB,codA,cynR,cynT,cynS,cynX,lacA,lacY,lacZ,lacI,mhpR,mhpA,mhpB,[mhpC]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 2561070 (0.000)39 (1.070) 37/474 0.1 100% coding (131/459 nt) gpt xanthine‑guanine phsophoribosyltransferase
?NC_000913 371382 = 0 (0.000)coding (144/867 nt) mhpC 2‑hydroxy‑6‑ketonona‑2,4‑dienedioate hydrolase

GGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCTGTGGCGTTCACCC  >  NC_000913/371132‑371611
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                      
aagacgcggactcatgtgaaatactggtttttagtgcgccagatctctataatctcgcgcaacctattttcccctcgaacactttttaagccgtagataaacaggctgggacacttcacatgagcgaaaaatacatcgtcacctgggacatgttgcagatccatgcacgtaaactcgcaagccgactgatgccttctgaacaatggaaaggcattattgccgtaagccgtggcggtctggtaccgggtgct                                                                                                                                                                                                                                        <  2:2596‑M2/1‑1 (MQ=255)
      cggactcatgtgaaatactggtttttagtgcgccagatctctataatctcgcgcaacctattttcccctcgaacactttttaagccgtagataaacaggctgggacacttcacatgagcgaaaaatacatcgtcacctgggacatgttgcagatccatgcacgtaaactcgcaagccgactgatgccttctgaacaatggaaaggcattattgccgtaagccgtggcggtctggtaccgggttcTActggc                                                                                                                                                                                                                                  >  2:26258‑M2/245‑251 (MQ=255)
       ggactcatgtgaaatactggtttttagtgcgccagatctctataatctcgcgcaacctattttcccctcgaacactttttaagccgtagataaacaggctgggacacttcacatgagcgaaaaatacatcgtcacctgggacatgttgcagatccatgcacgtaaactcgcaagccgactgatgccttctgaacaatggaaaggcattattgccgtaagccgtggcggtctggtaccgggtgcTActggct                                                                                                                                                                                                                                 <  1:50487‑M2/8‑1 (MQ=255)
                 gaaatactggtttttagtgcgccagatctctataatctcgcgcaacctattttcccctcgaacactttttaagccgtagataaacaggctgggacacttcacatgagcgaaaaatacatcgtcacctgggacatgttgcagatccatgcacgtaaactcgcaagccgactgatgccttctgaacaatggaaaggcattattgccgtaagccgtggcggtctggtaccgggtgcTACTGGCTGGGCGAACtt                                                                                                                                                                                                                       <  1:190706‑M2/18‑1 (MQ=255)
                            ttttagtgcgccagatctctataatctcgcgcaacctattttcccctcgaacactttttaagccgtagataaacaggctgggacacttcacatgagcgaaaaatacatcgtcacctgggacatgttgcagatccatgcacgtaaactcgcaagccgactgatgccttctgaacaatggaaaggcattattgccgtaagccgtggcggtctggtaccgggtgcTACTGGCTGGGCGAACTTCAGCCGCAat                                                                                                                                                                                                             <  2:50506‑M2/28‑1 (MQ=255)
                             tttagtgcgccagatctctataatctcgcgcaacctattttcccctcgaacactttttaagccgtagataaacaggctgggacacttcacatgagcgaaaaatacatcgtcacctgggacatgttgcagatccatgcacgtaaactcgcaagccgactgatgccttctgaacaatggaaaggcattattgccgtaagccgtggcggtctggtaccgggtgcTACTGGCTGGGCGAACTTCAGCCGCAata                                                                                                                                                                                                            <  2:323723‑M2/29‑1 (MQ=255)
                                        agatctctataatctcgcgcaacctattttcccctcgaacactttttaagccgtagataaacaggctgggacacttcacatgagcgaaaaatacatcgtcacctgggacatgttgcagatccatgcacgtaaacgcgcaagccgactgatgccttctgaacaatggaaaggcatgattgccgtaagccgtggcgggctggtaccgggggcGACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTgg                                                                                                                                                                                                <  1:208618‑M2/41‑1 (MQ=255)
                                                ataatctcgcgcaacctattttcccctcgaacactttttaagccgtagataaacaggctgggacacttcacatgagcgaaaaatacatcgtcacctgggacatgttgcagatccatgcacgtaaactcgcaagccgactgatgccttctgaacaatggaaaggcattattgccgtaagccgtggcggtctggtaccgggtgcTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCg                                                                                                                                                                                        <  2:855‑M2/49‑1 (MQ=255)
                                                  aatctcgcgcaacctattttcccctcgaacactttttaagccgtagataaacaggctgggacacttcacatgagcgaaaaatacatcgtcacctgggacatgttgcagatccatgcacgtaaactcgcaagccgactgatgccttctgaacaatggaaaggcattattgccgtaagccgtggcggtctggtaccgggtgcTACTGGCTGGGCGAACTTCAGCCGCAATCTCGATCCGCTGGTAGAGGCggg                                                                                                                                                                                      >  2:167722‑M2/201‑251 (MQ=255)
                                                       cgcgcaacctattttcccctcgaacactttttaagccgtagataaacaggctgggacacttcacatgagcgaaaaatacatcgtcacctgggacatgttgcagatccatgcacgtaaactcgcaagccgactgatgccttctgaacaatggaaaggcattattgccgtaagccgtggcggtctggtaccgggtgcTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCt                                                                                                                                                                                    <  1:47881‑M2/53‑1 (MQ=255)
                                                       cgcgcaacctattttcccctcgaacactttttaagccgtagataaacaggctgggacacttcacatgagcgaaaaatacatcgtcacctgggacatgttgcagatccatgcacgtaaactcgcaagccgactgatgccttctgaacaatggaaaggcattattgccgtaagccgtggcggtctggtaccgggtgcTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATc                                                                                                                                                                                 <  1:113873‑M2/56‑1 (MQ=255)
                                                           caacctattttcccctcgaacactttttaagccgtagataaacaggctgggacacttcacatgagcgaaaaatacatcgtcacctgggacatgttgcagatccatgcacgtaaactcgcaagccgactgatgccttctgaacaatggaaaggcattattgccgtaagccgtggcggtctggtaccgggtgcTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTATAggcgggctcttgggg                                                                                                                                                                             >  2:331709‑M2/192‑244 (MQ=255)
                                                              cctattttcccctcgaacactttttaagccgtagataaacaggctgggacacttcacatgagcgaaaaatacatcgtcacctgggacatgttgcagatccatgcacgtaaactcgcaagccgactgatgccttctgaacaatggaaaggcattattgccgtaagccgtggcggtctggtaccgggtgcTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGAt                                                                                                                                                                          <  1:139000‑M2/63‑1 (MQ=255)
                                                                      cccctcgaacactttttaagccgtagataaacaggctgggacacttcacatgagcgaaaaatacatcgtcacctgggacatgttgcagatccatgcacgtaaactcgcaagccgactgatgccttctgaacaatggaaaggcattattgccgtaagccgtggcggtctggtaccgggtgcTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGGGATCCTGCTgg                                                                                                                                                                  >  1:187751‑M2/181‑251 (MQ=255)
                                                                      cccatagaacactttttaagccgtagataaacaggctgggacacttcacaggagcgaaaaatacatcgtcacctgggacatgttgcagatccatgcacgtaaactcgcaagccgactgatgccttctgaacaagggaaaggcattattgccgtaagccgtggcggtctggtaccgggtgcTACTGGCTGGGCGAACGTCAGCCGCAATATCGAGCCGCTGGTAGAGGCGTTCTATCGGGTGATCCGGCTgg                                                                                                                                                                  <  2:29827‑M2/71‑1 (MQ=255)
                                                                              acactttttaagccgtagataaacaggctgggacacttcacatgagcgaaaaatacatcgtcacctgggacatgttgcagatccatgcacgtaaactcgcaagccgactgatgccttctgaacaatggaaaggcattattgccgtaagccgtggcggtctggtaccgggggcTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTCTAGGGGGGTCCTGCTGGATTGTCCg                                                                                                                                                          >  2:137546‑M2/173‑251 (MQ=255)
                                                                                                  aaacaggctgggacacttcacatgagcgaaaaatacatcgtcacctgggacatgttgcagatccatgcacgtaaactcgcaagccgactgatgccttctgaacaatggaaaggcattattgccgtaagccgtggcggtctggtaccgggtgcTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCCATCGGGGGGTCCTGCTGGGTTGTCCGGGTTGGGGCAAAAGCGATTc                                                                                                                                      >  1:344989‑M2/153‑251 (MQ=255)
                                                                                                     caggctgggacacttcacatgagcgaaaaatacatcgtcacctgggacatgttgtagatccatgcacgtaaactcgcaagccgactgatgccttctgaacaatggaaaggcattattgccgtaagccgtggcggtctggtaccgggtgcTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGt                                                                                                                                   <  1:252751‑M2/102‑1 (MQ=255)
                                                                                                         ctgggacacttcacatgagcgaaaaatacatcgtcacctgggacatgttgcagatccatgcacgtaaactcgcaagccgactgatgccttctgaacaatggaaaggcattattgccgtaagccgtggcggtctggtaccgggtgct                                                                                                                                                                                                                                        <  1:210587‑M2/1‑1 (MQ=255)
                                                                                                         ctgggacacttcacatgagcgaaaaatacatcgtcacctgggacatgttgcagatccatgcacgtaaactcgcaagccgactgatgccttctgaacaatggaaaggcattattgccgtaagccgtggcggtctggaaccgggtgct                                                                                                                                                                                                                                        >  2:210587‑M2/146‑146 (MQ=255)
                                                                                                             gacacttcacatgagcgaaaaatacatcgtcacctgggacatgttgcagatccatgcacgtaaactcgcaagccgactgatgccttctgaacaatggaaaggcattattgccgtaagccgtggcggtctggtaccgggtgcTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGGGATCCTGCTGGGTTGGCCGGGGTGGGGCAAGAGCGATTCGGGCGTTAAt                                                                                                                            >  1:147617‑M2/142‑250 (MQ=255)
                                                                                                                            cgaaaaatacatcggcacctgggacatgttgcagatccatgcacgtaaactcgcaagccgactgatgccttctgaacaatggaaaggcattattgccgtaagccgtggcggtctggtaccgggtgcTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCgg                                                                                                            <  2:15718‑M2/125‑1 (MQ=255)
                                                                                                                               aaaatacatcgtcacctgggacatgttgcagatccatgcacgtaaactcgcaagccgactgatgccttctgaacaatggaaaggcattattgccgtaagccgtggcggtctggtaccgggtgcTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATc                                                                                                         <  1:76714‑M2/128‑1 (MQ=255)
                                                                                                                                         gtcacctgggacatgttgcagatccatgcacgtaaactcgcaagccgactgatgccttctgaacaatggaaaggcattattgccgtaagccgtggcggtctggtaccgggtgcTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTTTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCa                                                                                                 >  1:273162‑M2/114‑249 (MQ=255)
                                                                                                                                         gtcacctgggacatgttgcagatccatgcacgtaaactcgcaagccgactgatgccttctgaacaatggaaaggcattattgccgtaagccgtggcggtctggtaccgggtgcTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGGTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGGTCGCGATCGGGTCTTAatgcccg                                                                                               >  1:20689‑M2/114‑248 (MQ=255)
                                                                                                                                             cctgggacatgttgcagatccatgcacgtaaactcgcaagccgactgatgccttctgaacaatggaaaggcattattgccgtaagccgtggcggtctggtaccgggtgcTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGGGGTCCTGCTGGGTTTTCCGGGGTGGGGGAAGAGCGGTTTGGGCGTTTATAGGGGGTTGCGGTCGGATCTTTATTCACGAATc                                                                                           >  2:103688‑M2/110‑251 (MQ=255)
                                                                                                                                                  gacatgttgcagctccatgcacgtaaactcgcaagccgactgatgccttctgaacaatggaaaggcattattgccgtaagccgtggcggtctggtaccgggtgcTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGGTCGCGATCGAATCTTAATGCACGAATCCTGa                                                                                       >  1:266293‑M2/105‑250 (MQ=255)
                                                                                                                                                  gacatgttgcagatccatgcacgtaaactcgcaagccgactgatgccttctgaacaatggaaaggcattattgccgtaagccgtggcggtctggtaccgggtgcTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGaa                                                                                      <  2:355334‑M2/147‑1 (MQ=255)
                                                                                                                                                                 catgcacgtaaactcgcaagccgactgatgccttctgaacaatggaaaggcattattgccgtaagccgtggcggtctggtaccgggtgcTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATc                                                                        >  1:38325‑M2/90‑250 (MQ=255)
                                                                                                                                                                            actcgcaagccgactgatgccttctgaacaatggaaaggcattattgccgtaagccgtggcggtctggtaccgggtgcTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTTATAGTGGTTCGCGGTCGGGTCTTAATGCACGAATCCTGGAAAACGTGGTGGATCAACTGGGTAtcc                                                            >  1:52662‑M2/79‑250 (MQ=255)
                                                                                                                                                                            actcgcaagccgactgatgccttctgaacaagggaaaggcattatggccgtaagccgtggcggtctggtaccgggtgcTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCg                                                            <  2:196661‑M2/173‑1 (MQ=255)
                                                                                                                                                                                           gatgccttctgaacaatggaaaggcattattgccgtaagccgtggcggtctggtaccgggtgcTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGAATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACctg                                             <  2:266293‑M2/187‑1 (MQ=255)
                                                                                                                                                                                              gccttctgaacaatggaaaggcattattgccgtaagccgtggcggtctggtaccgggtgcTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTgg                                         <  2:366204‑M2/191‑1 (MQ=255)
                                                                                                                                                                                                      aacaatggaaaggcattattgccgtaagccgtggcggtctggtaccgggtgcTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCg                                 <  1:194553‑M2/199‑1 (MQ=255)
                                                                                                                                                                                                               aaggcattattgccgtaagccgtggcggtctggtaccgggtgcTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTTGGGCACGAGCGCTTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTTATGCACGAATCCTTAAAAGCGTGGTGGGTCAACTg                                                                   >  1:245428‑M2/44‑209 (MQ=255)
                                                                                                                                                                                                               aaggcattattgccgtaagccgtggcggtctggtaccgggtgcTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTg                                                                   <  2:245428‑M2/166‑1 (MQ=255)
                                                                                                                                                                                                                  gcattattgccgtaagccgtggcggtctggtaccgggtgcTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGTGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATTCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCAt                     <  1:3195‑M2/211‑1 (MQ=255)
                                                                                                                                                                                                                           ccgtaagccgtggcggtctggtaccgggtgcTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCtgtg            <  1:300067‑M2/220‑1 (MQ=255)
                                                                                                                                                                                                                               aagccgtggcggtctggtaccgggtgcTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGAGCCTGCTGGGTTGTCCGGGGTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCTGTGGCGt        >  1:50680‑M2/28‑251 (MQ=255)
                                                                                                                                                                                                                                     tggcggtctggtaccGGGAGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGGTATCGGGGGATCCCGCTGGAGTTTCCGGGTTGGGGGAAAAGCGATTCGGGCGGTTATAGTGGTTCGCCGTCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCCACTGGGTATCGCCAAAATCCCACCGCTGGGAAAATCGGGGGGGGGCCATAGTCCTGTGGCGTTCAccc  >  2:248154/14‑251 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                      
GGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCTGCGCATCCATTTTAATGACTGCGGACAAGGCGACGAAACCGTTGTCCTGCTGCATGGTTCCGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCTGTGGCGTTCACCC  >  NC_000913/371132‑371611

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