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 67 [0] [0] 67 [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)67 (1.060) 63/482 0.1 100% coding (893/960 nt) aes acetylesterase
?NC_000913 551330 = 0 (0.000)coding (892/894 nt) ybcF putative carbamate kinase

ATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTGGTAAAGCAGACGGCTGTCATCCAGCAGCGGATCGAACTCCGCCCCGGCAATAAAACAGGGCGGAACTTCGCGAGTGAGATCATTATTAAACAGACAGTAATACGGCGACTCGCGGTCCGCGTCGTTGCTTAAATATGCCTCTTCGTACATCTGCAAATCCTGTTGCG  >  NC_000913/498838‑499324
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   
atTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgt                                                                                                                                                                                                                                              <  1:481653‑M1/251‑8 (MQ=255)
    aGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTTGAAGAGATTGCCGAGCAAAACCGTTAGGTCTTCCTCGGTGCCGGCGATAAAAAATATGAATAGATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTGGTTGCCGCGTAGCACTAAAGCTGAGCGGTGAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTGTCGTCATCCGTTtattcgtaggc                                                                                                                                                                                                                                          <  2:337681‑M1/251‑12 (MQ=255)
           gtgCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagac                                                                                                                                                                                                                                   <  1:252805‑M1/251‑19 (MQ=255)
             gCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAAAGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGGAAAGAAATGAGCGCCGTCGCGAAGAGCCCCGGCGGCGGGTTTCATCATCCGTGtagtcgtaggccttagacat                                                                                                                                                                                                                                 >  1:92036‑M1/1‑231 (MQ=255)
              cccGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatt                                                                                                                                                                                                                                <  1:261975‑M1/251‑22 (MQ=255)
                        gctgcGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgat                                                                                                                                                                                                                      <  2:483815‑M1/251‑32 (MQ=255)
                            cGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcga                                                                                                                                                                                                                  <  1:25044‑M1/251‑36 (MQ=255)
                             ggATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGGTTTCATTATCCGGGgagtcggaggcattaaacaattttgcctgatgcgac                                                                                                                                                                                                                 >  1:406922‑M1/1‑215 (MQ=255)
                               aTTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCGCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgc                                                                                                                                                                                                               <  2:493443‑M1/251‑39 (MQ=255)
                                        agaCGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgt                                                                                                                                                                                                      <  2:395502‑M1/251‑48 (MQ=255)
                                               ggAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatc                                                                                                                                                                                                <  2:406922‑M1/250‑54 (MQ=255)
                                                                aaaaCCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcggcttatcaggcctacaaccggtgcc                                                                                                                                                                              >  1:138683‑M1/1‑180 (MQ=255)
                                                                 aaaCCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTTCGGCGCTTAATGCCACGCCGGAACCTATCGAAATGATGGCTAATCTTGTTGCCGCGCCTCGCTAAAGCTTAGCGGGAAAGAACTGAGCGCCGTCGCGGAGAGCCTCGGCGGCGGGTTTCATCCTCCGTGtagtcggtgggattagacacttttgcctgatgcgacgcttgacgcgcattattaggcctacaaccggtgccg                                                                                                                                                                             >  1:400101‑M1/1‑179 (MQ=255)
                                                                   aCCGTGAGGTCCTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAAGGCCACTCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCAGAGCGGTAAAGAACTGATCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcatgcctacaaccggtgccgca                                                                                                                                                                           <  1:208717‑M1/251‑75 (MQ=255)
                                                                     cGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTTCGGCGCTTTATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGGGGTTTTTATCATCCGTGtaggtgtaggccttagacatttttgcctgatgggaggcttggcgcgtctttttaggcctccaaccggtgcagcctc                                                                                                                                                                         >  1:17578‑M1/1‑175 (MQ=255)
                                                                                ccTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattgg                                                                                                                                                               <  2:166152‑M1/250‑87 (MQ=255)
                                                                                               gAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatg                                                                                                                                               >  1:294117‑M1/1‑149 (MQ=255)
                                                                                                   aaaTATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgat                                                                                                                                           >  1:92147‑M1/1‑145 (MQ=255)
                                                                                                   aaaTATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgat                                                                                                                                           >  1:57706‑M1/1‑145 (MQ=255)
                                                                                                   aaaTATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTGGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggggccgcatccggcaattggtgcacaatgcctgatgcgat                                                                                                                                           <  1:162594‑M1/251‑107 (MQ=255)
                                                                                                   aaaTATGAATAGATGCCGGCGCTTAAGGACACGCCGAAACATATCGAAATGATGGCGAATCTTGTGGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCTTCGCGAAGAGCCTCGGCGGCGGTTGTCATCATCCGTGtagtcgtagtcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctactaccggtgccgcatccggcaattggtgcacaatgcctgatgcgat                                                                                                                                           <  2:457905‑M1/251‑107 (MQ=255)
                                                                                                           atatatTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGCTGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcc                                                                                                                                                                                                                          >  1:504022‑M1/1‑137 (MQ=255)
                                                                                                           atatatTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCTGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatgtgtgcc                                                                                                                                                                                                                          <  2:504022‑M1/164‑28 (MQ=255)
                                                                                                             atatTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgc                                                                                                                                  >  1:334895‑M1/1‑135 (MQ=255)
                                                                                                                   cGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttat                                                                                                                           >  1:196489‑M1/1‑129 (MQ=255)
                                                                                                                   cGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcccgatgcgatgcttgacgcgtcttat                                                                                                                           <  1:187217‑M1/251‑123 (MQ=255)
                                                                                                                            aTGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctac                                                                                                                  >  1:402261‑M1/1‑120 (MQ=255)
                                                                                                                              gCCACGCCGGAACATATCGAAATGATGGCGAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACGGAGCGCCGTCGCGAAGAGCCGCGGCGGCGGTTTTCAGCAGCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcgtgacgcgtcttatcaggccgacaaccggtgcctcatccgtcaattgttgcacaatgcctgatgcgatgcttgacgcatcgtagcagtcctaca                                                                                                                 <  2:276071‑M1/250‑133 (MQ=255)
                                                                                                                               ccACGCCGGCAAATATCGAACTGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGCACTGAGCGCCGTCGCGCAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagccatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgctgcttgacgcatcttatcaggcctacaat                                                                                                               >  2:553071‑M1/1‑117 (MQ=255)
                                                                                                                                      ggAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTAGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtacc                                                                                                        >  1:398879‑M1/1‑110 (MQ=255)
                                                                                                                                        aaCATATCGAACTGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccgg                                                                                                      >  2:484307‑M1/1‑108 (MQ=255)
                                                                                                                                               cGAAATGATGGCGAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggta                                                                                               <  2:402261‑M1/251‑151 (MQ=255)
                                                                                                                                                gAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAACGAAATGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggt                                                                                                >  2:507651‑M1/1‑100 (MQ=255)
                                                                                                                                                 aaaTGATGTCTAATCTTGTTGCCGCGGATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcggcttatcaggcctacaaccggagccgcatccggcaattagtgcaccatgcctgatgcgatgcttgacgcatcttatcaggcctacaatggggaccggatcggtagg                                                                                             >  2:447721‑M1/1‑99 (MQ=255)
                                                                                                                                                 aaaTGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtagg                                                                                             >  1:425313‑M1/1‑99 (MQ=255)
                                                                                                                                                      atgGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggg                                                                                        >  2:180151‑M1/1‑94 (MQ=255)
                                                                                                                                                        ggCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatctgatcaggcctacaatgggtaccggatcggtaggccggat                                                                                       <  1:383433‑M1/250‑159 (MQ=255)
                                                                                                                                                            aaTCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGCGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggc                                                                                  >  1:268797‑M1/1‑88 (MQ=255)
                                                                                                                                                            aaTCTTGTTGCCGCGTAGCGCTAAAGCGGAGCGGTAAAGAACTGGGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATAATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaaggggtaccggatcggtaggccggataagg                                                                                   <  2:14191‑M1/250‑163 (MQ=255)
                                                                                                                                                                ttgttgCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggaaaaggcgttt                                                                              >  2:93296‑M1/1‑84 (MQ=255)
                                                                                                                                                                        gcgTATCGCTAAAGCGGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggatagggcgtttacgccgca                                                                      <  2:48448‑M1/251‑176 (MQ=255)
                                                                                                                                                                              cGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggc                                                                <  1:565743‑M1/251‑182 (MQ=255)
                                                                                                                                                                                 aaagtgttggcggtaaagaGATGAGCGCCGGCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGgagtcgtaggcattaggcatgtgggcctgatgcgccgcttggcgcgtctgatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaattggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaag                                                             <  2:102932‑M1/242‑185 (MQ=255)
                                                                                                                                                                                      cTGAGCGGTCAAGAACTGAGCGCCGTAGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaataat                                                          >  1:461741‑M1/1‑62 (MQ=255)
                                                                                                                                                                                        gAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaataga                                                       <  1:71639‑M1/250‑191 (MQ=255)
                                                                                                                                                                                           cGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGCTTTCATAATCCGTGcagtcgtaggcattagacatttgtgcctgatgcgacgcttgccgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcctcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttccgccgcatccggcaagaatagagcac                                                   >  2:266470‑M1/1‑57 (MQ=255)
                                                                                                                                                                                            ggTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgctgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcacc                                                  >  1:108727‑M1/1‑56 (MQ=255)
                                                                                                                                                                                            ggTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcacc                                                  >  1:270918‑M1/1‑56 (MQ=255)
                                                                                                                                                                                               aaaGAACGGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcaccagt                                               <  2:92036‑M1/251‑199 (MQ=255)
                                                                                                                                                                                                aagaccggagcgccgTCGCGAAGAGCCTCGGCGGCGGTGTTCAGCATCCGTTgagtcgtaggctttaggcatgtgtgcctgatgcgacgcttcacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggttcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcaccagtt                                              <  2:398879‑M1/244‑200 (MQ=255)
                                                                                                                                                                                                aaGAACTGCGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcctccggcaagcatagagcacccgtt                                              >  2:265082‑M1/1‑52 (MQ=255)
                                                                                                                                                                                                    aCTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtaggcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcagtcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatatagcaccagttaacc                                          <  2:400101‑M1/251‑204 (MQ=255)
                                                                                                                                                                                                       gAGCGCCGTCGCGAAGCGCCTCCGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggcaggccggataaggcgtttacgccgcatccggcaagaatagagcaccagttaaccga                                        >  1:115549‑M1/1‑45 (MQ=255)
                                                                                                                                                                                                       gAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcaccagttaaccgaa                                       <  1:451059‑M1/251‑207 (MQ=255)
                                                                                                                                                                                                        aGCGCCGTCGCGAAGCGACTCGGAGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgccgcgtcttatcaggactacaaccggtgacgcctccggcaattggtgcacaatgcctgatgcgccgcttgacgcatcttatcaggcataccatgggtaacggatcggtaggccggataaggcgtttacgccgcctccggcaagaatagagaacccgttaacctaac                                      >  2:337039‑M1/1‑44 (MQ=255)
                                                                                                                                                                                                         gcgcCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtacc                                                                                                        <  1:526387‑M1/184‑142 (MQ=255)
                                                                                                                                                                                                         gcgcCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtacc                                                                                                        >  2:526387‑M1/1‑43 (MQ=255)
                                                                                                                                                                                                           gcCGTCGCGAAGAGCCTCGGCGGCGGTTTTCCTAATCCGTGcagtcgtaggtattcgacatttgtgcatgatgcgacgcttgccgcctcttataaggcctacaaccggtgccgcatccggccattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccgggacaggcggttacgccgcatccggcaagactagagcaccagttaaccgcaattc                                   >  2:511985‑M1/1‑41 (MQ=255)
                                                                                                                                                                                                                    aaGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtatgccggataaggcgtttacgccgcatccggcaagaatagagcaccagttaaccgaacttactctgcgcc                          <  1:252452‑M1/251‑220 (MQ=255)
                                                                                                                                                                                                                      gagCCTCGGCGGCGGTTTTCATCATCCGTAtcgtcgtaggccttagacatttgtgcctgatgccacgcttgacgcgtcttaacaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcacaagttaaccgaacttactctgcgccc                         >  2:289441‑M1/1‑30 (MQ=255)
                                                                                                                                                                                                                         ccTCGGCGGCGGTTCTCATCATCCGTGtagtcgtaggcattcgacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaacaggcgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcaccagttaaccgaacttactctgcgcccaaat                     >  2:240585‑M1/1‑27 (MQ=255)
                                                                                                                                                                                                                                ggcggTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaacggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatccggcctacaatgggtaccggatcggtaggccggataaggagtttacgccgcatccggcaagaatagagcacccgttaaccgaacttactctgcgcccaaatcacgccc              >  2:394372‑M1/1‑20 (MQ=255)
                                                                                                                                                                                                                                  gggtTTTCATCATCCGTTgaggcgtaggcattagacattttttcctgaggcgacgctggacgcgtcttgtcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcttacaatgggtaccggagcggtaggccggataaggcgtttgcgccgcatccggcaagaaggcagcaccagtgaaccgaacttacgcttcgcccaaaccacgccgct            <  2:268797‑M1/250‑234 (MQ=255)
                                                                                                                                                                                                                                      tttCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagctccagttaaccgaacttactctgcgcccaaatcacgccgctggc         <  1:72304‑M1/250‑237 (MQ=255)
                                                                                                                                                                                                                                        tcatcaTCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcaccagttaaccgaacttactctgcgcccaaatcacgccgctggcat       <  1:182965‑M1/250‑239 (MQ=255)
                                                                                                                                                                                                                                         catcatCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttctccggcctacaaccggtgccgcctccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcaccagttaaacgaacttactctgcgcccaaatcacgccgctggcatat     >  1:568788‑M1/1‑11 (MQ=255)
                                                                                                                                                                                                                                            catcCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcaccagttaaccgaacttactctgcgcccaaatcacgccgctggcatattcc  <  1:138551‑M1/251‑244 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   
ATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTGGTAAAGCAGACGGCTGTCATCCAGCAGCGGATCGAACTCCGCCCCGGCAATAAAACAGGGCGGAACTTCGCGAGTGAGATCATTATTAAACAGACAGTAATACGGCGACTCGCGGTCCGCGTCGTTGCTTAAATATGCCTCTTCGTACATCTGCAAATCCTGTTGCG  >  NC_000913/498838‑499324

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