Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 4,139,771 10483 bp→TG [ptsA][eptC] [ptsA], frwC, frwB, pflD, pflC, frwD, yijO, [eptC]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 4139771 4150253 10483 51 [0] [0] 54 [ptsA]–[eptC] [ptsA],frwC,frwB,pflD,pflC,frwD,yijO,[eptC]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 41397700 (0.000)51 (0.670)
+TG
48/490 0.3 100% coding (2452/2502 nt) ptsA putative PTS multiphosphoryl transfer protein PtsA
?NC_000913 4150254 = 0 (0.000)coding (12/1734 nt) eptC phosphoethanolamine transferase EptC

TGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCGGCGTTAACCAGCGCGTTACGGAATTCTTCGTGCATGATGCGGCGAGCGAGGCGCGAGAAAATGCGCATATGCTGATCGCCCGCAGCGTGTTTGTTCAGGGTTAACATAATGATGAATTGCGCTTCATCATCGCCCCAGCGCACCGGCGCTTGCAGACGCGCCACGCTGATGGTGGATTGCTCAATGTGTTCTGATTTGCTGTGTGGAATGGCAAAACTAAAGCCCAGACCGGTAGAGAAAACGGCCT  >  NC_000913/4139528‑4140018
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                        
gggtCTGAGCAAACAGAATGCCCTCATCACCAATCGCTGTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtgga                                                                                                                                                                                                                                                  <  2:589093‑M1/250‑9 (MQ=255)
           acgatcctgccctcaTCACCAAGCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagta                                                                                                                                                                                                                                        <  2:203157‑M1/243‑19 (MQ=255)
                  tGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagta                                                                                                                                                                                                                                        >  1:203157‑M1/1‑225 (MQ=255)
                  ggcCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCAATATGGACGCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTAGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgggtggaatgcatagtaagaacctg                                                                                                                                                                                                                                <  1:225750‑M1/250‑27 (MQ=255)
                   gCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGTCTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTGCCATGATGTTTCCTTTGACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgt                                                                                                                                                                                                                               <  2:2737‑M1/251‑28 (MQ=255)
                    cccTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtc                                                                                                                                                                                                                              <  1:255238‑M1/251‑29 (MQ=255)
                     ccTCATCCACAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtct                                                                                                                                                                                                                             <  2:583302‑M1/251‑30 (MQ=255)
                     ccTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGAAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtct                                                                                                                                                                                                                             <  2:30855‑M1/251‑30 (MQ=255)
                        catcaCCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCGCTGGTTAGTTGTCACCCCGGCAATGTGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGTTGTCCAGATACAGTTCCATGATGTTTACTGTTGCAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGACATAGGCtgtgtggaattcatagtaagatcctgtctt                                                                                                                                                                                                                            <  2:175452‑M1/249‑31 (MQ=255)
                        aatcacCACTCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTTTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGCTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGATGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttg                                                                                                                                                                                                                           <  2:374323‑M1/249‑32 (MQ=255)
                              cAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaat                                                                                                                                                                                                                    <  1:500592‑M1/251‑39 (MQ=255)
                                aTCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatat                                                                                                                                                                                                                  >  1:124094‑M1/1‑211 (MQ=255)
                                           cAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgggtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgt                                                                                                                                                                                                        >  2:314486‑M1/1‑200 (MQ=255)
                                            aGACGCGGCCGCACTTCTCATATGGACTCATTGCTGGCCGCGATAATGCTCGTGTTAGTTGGAAACACGGCAATGGGGACTACGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGGGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCggtgtggaatgcatagtacgaacccgtcttgaaaaaatatcgccgaatgaa                                                                                                                                                                                                       >  2:576038‑M1/1‑199 (MQ=255)
                                                gcgGCAGCACTTCCCATTTGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgac                                                                                                                                                                                                  <  2:115157‑M1/251‑57 (MQ=255)
                                                      gcaCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattcc                                                                                                                                                                                            <  2:51198‑M1/251‑63 (MQ=255)
                                                      gcaCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattcc                                                                                                                                                                                            <  2:208664‑M1/251‑63 (MQ=255)
                                                            cccATATGGCCTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaagg                                                                                                                                                                                      >  1:40105‑M1/1‑183 (MQ=255)
                                                                   ggACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctg                                                                                                                                                                               <  1:334036‑M1/251‑76 (MQ=255)
                                                                               ggCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatatt                                                                                                                                                                     <  2:31651‑M1/249‑86 (MQ=255)
                                                                                     gATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaattt                                                                                                                                                             <  1:399781‑M1/251‑94 (MQ=255)
                                                                                     gATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaa                                                                                                                                                                             <  2:72185‑M1/235‑78 (MQ=255)
                                                                                     gATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaa                                                                                                                                                                             >  1:72185‑M1/1‑158 (MQ=255)
                                                                                       tAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttg                                                                                                                                                            >  1:522515‑M1/1‑156 (MQ=255)
                                                                                         aTGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaat                                                                                                                                                         <  1:398627‑M1/251‑98 (MQ=255)
                                                                                           gCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaata                                                                                                                                                        <  2:522515‑M1/250‑99 (MQ=255)
                                                                                                       tGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagac                                                                                                                                              >  2:367989‑M1/1‑140 (MQ=255)
                                                                                                            ccccGGCAACGGGGAAGATGCTTGCCAGACGTTCGACTTCTGCGACGTTCGCGGTGTCCAGATACAGTTCCAGGATTTTTCCTTTTACAGGTCCAGTTCATGTTGCAGCAGGATGGCGATTGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaagctgtcttgaaaaaatatcgccgaatgtaacgacaattctttaaggatatctgaaggtagattcagaattggaataaaattcagacagagatat                                                                                                                                      <  2:524273‑M1/251‑117 (MQ=255)
                                                                                                                ggCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattg                                                                                                                                  <  1:97896‑M1/251‑121 (MQ=255)
                                                                                                                    aTGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaa                                                                                                                              <  1:113335‑M1/251‑125 (MQ=255)
                                                                                                                               gCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgtt                                                                                                                    <  2:349384‑M1/250‑135 (MQ=255)
                                                                                                                                    ccAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaac                                                                                                               >  1:353307‑M1/1‑111 (MQ=255)
                                                                                                                                      ggaCGTTCTACTGCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACATTTCCAGTTCATGTTGCAGCAGGCTGGCTATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtctttaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacaga                                                                                                            <  1:587556‑M1/250‑143 (MQ=255)
                                                                                                                                          gttcgcattctgcgacGTTAGCGGTGTCAAGATACAGTGCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTGGCAGCAGGCTGACGATATCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaattaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagta                                                                                                        <  2:622552‑M1/244‑147 (MQ=255)
                                                                                                                                                       gACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatagattgaaaacgagggtgttagaacagaagtatttcagaaaaccc                                                                                           <  1:540495‑M1/251‑160 (MQ=255)
                                                                                                                                                         cGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctc                                                                                         <  1:293893‑M1/251‑162 (MQ=255)
                                                                                                                                                          gTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcg                                                                                        <  1:26884‑M1/251‑163 (MQ=255)
                                                                                                                                                               cGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCCGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgca                                                                                    >  2:540973‑M1/1‑84 (MQ=255)
                                                                                                                                                                     ccAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatataatgaaaacgagggtgttagaccagaagtatttcagaaaaccctcgcgcaaaagctc                                                                             <  2:392962‑M1/251‑174 (MQ=255)
                                                                                                                                                                            cAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtt                                                                      <  1:374224‑M1/251‑181 (MQ=255)
                                                                                                                                                                             aGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggttt                                                                     <  1:244525‑M1/251‑182 (MQ=255)
                                                                                                                                                                                 ccATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcag                                                                 <  2:32893‑M1/251‑186 (MQ=255)
                                                                                                                                                                                       tgTTTCCTTTTACAGTTCAAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCTGCAGAGGCtgtgtggaatgcatcgtcagaacctgtcttgaaaaaatatcgccgaatgtaccgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagaaagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagagg                                                           >  2:414151‑M1/1‑60 (MQ=255)
                                                                                                                                                                                        gTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagagg                                                           >  1:218990‑M1/1‑59 (MQ=255)
                                                                                                                                                                                         tttCCTTTTACAGTTCCAGTTCATGTTGCAGCAGCATGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcctagaaagaacctgtcttgaaaaaatatcgccgaatgaaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagtaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagagga                                                          >  1:542546‑M1/1‑58 (MQ=255)
                                                                                                                                                                                                   cagttccagttcATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagaggaagattagccgg                                               >  2:17170‑M1/1‑48 (MQ=255)
                                                                                                                                                                                                        ccagttcATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagaggaagattagccggtatt                                           <  1:455072‑M1/250‑208 (MQ=255)
                                                                                                                                                                                                              cATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagaggaagattagccggtattacgcat                                     <  1:169901‑M1/250‑214 (MQ=255)
                                                                                                                                                                                                               aTGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagaggaagattagccggtattacgcatac                                   >  1:558256‑M1/1‑36 (MQ=255)
                                                                                                                                                                                                                                cGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagaggaagattagccggtattacgcatacctgccgcaatcccggc                   >  1:169752‑M1/1‑19 (MQ=255)
                                                                                                                                                                                                                                                ggCtgtgtggaatgcatagtaagaaccagtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacaagggtttgcagaagaggaagattagccggtattacgcatacctgccgcaatcccggcaatagtgaccattaacg  <  1:388187‑M1/251‑249 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                        
TGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCGGCGTTAACCAGCGCGTTACGGAATTCTTCGTGCATGATGCGGCGAGCGAGGCGCGAGAAAATGCGCATATGCTGATCGCCCGCAGCGTGTTTGTTCAGGGTTAACATAATGATGAATTGCGCTTCATCATCGCCCCAGCGCACCGGCGCTTGCAGACGCGCCACGCTGATGGTGGATTGCTCAATGTGTTCTGATTTGCTGTGTGGAATGGCAAAACTAAAGCCCAGACCGGTAGAGAAAACGGCCT  >  NC_000913/4139528‑4140018

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