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 45 [0] [0] 44 [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)44 (0.860)
+TG
42/488 0.1 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

ACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTAGAGCAGTGCCCAACCCAGGGCTTTCCAGCTAAAAAGAGGTTTAGCCTGGACTTCTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAA  >  NC_000913/4150007‑4150497
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   
catggtctgagcaaacagaatgccctcatcaccaatcgctttttgcagacgcggcagcacttcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgtgtg                                                                                                                                                                                                                                                  <  1:392585‑M2/4‑1 (MQ=255)
  tggtctgagcaaacagaatgccctcatcaccaatcgctttttgcagacgcggcagcacttcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTgg                                                                                                                                                                                                                                                 >  1:60403‑M2/246‑250 (MQ=255)
        gagcaaacagaatgccctcatcaccaatcgctttttgcagacgcggcagcacttcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCAt                                                                                                                                                                                                                                          <  1:439456‑M2/12‑1 (MQ=255)
         agcaaacagaatgccctcatcaccaatcgctttttgcagacgcggcagcacttcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCAt                                                                                                                                                                                                                                          >  1:79858‑M2/239‑250 (MQ=255)
           caaacagaatgccctcatcaccaatcgcttgttgcagacgcggcagcacttcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatattcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccgtttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGt                                                                                                                                                                                                                                       <  2:357858‑M2/15‑1 (MQ=255)
                                        ttttgcagacgcgacagcacttcccatatggactccttgctggcagcgatgatgctcgggttagttgtcgccccggcaatggttaatatgcgttccagacgttcgacttctgcgtcgctagcggtgtccatatacagttccatgatgtttccttttacagttccagttcatgctgcagcaggctggcgatagcgtctgcataggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGc                                                                                                                                                                                                             <  2:329709‑M2/41‑1 (MQ=255)
                                           tgcagacgcggcagcacttcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAAt                                                                                                                                                                                                        <  1:431610‑M2/46‑1 (MQ=255)
                                                       agcacttcccatatggactccttgctggcagcgaaaatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgtcgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTc                                                                                                                                                                                           <  1:259355‑M2/59‑1 (MQ=255)
                                                                catatggactccttgctggcagcgagtatgctcgggttagttgtgaccccggcaatggggaatatgcgtgccgggcgttcgacttctgcgacgttagcggtgttcagatacagttccatgatgtgtccttttacagttccagttcatgttgcatcaggctggcgatagcgtcgtcagaggcggTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGat                                                                                                                                                                                  <  2:242147‑M2/68‑1 (MQ=255)
                                                                 atatggactccttgctggcagcgataatgctcgggttagttgtcaccacggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGCTTGTAACGACAATTCCTTAAGGat                                                                                                                                                                                  >  1:413912‑M2/183‑250 (MQ=255)
                                                                   atggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATc                                                                                                                                                                               >  2:441876‑M2/181‑251 (MQ=255)
                                                                          ccttgctggcagcgataatgctcggcttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagctccatgatgtttcctgttacagtgccagtttatgttgcagcaggctgtccctagcttcggcagaggctgTGTGGAATGCATCGTCAGAGCCTGTCTTGAAAAAATATCGCCGAATGGAACGACATTTCCTTAAGGATCTCTGAAGGt                                                                                                                                                                        <  2:382994‑M2/78‑1 (MQ=255)
                                                                              gctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTatat                                                                                                                                                                    <  2:22769‑M2/82‑1 (MQ=255)
                                                                                 ggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCTCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATAtt                                                                                                                                                                   <  2:342501‑M2/83‑1 (MQ=255)
                                                                                 ggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTc                                                                                                                                                                  <  1:256822‑M2/84‑1 (MQ=255)
                                                                                   cagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGa                                                                                                                                                               <  1:354425‑M2/87‑1 (MQ=255)
                                                                                    agcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGaa                                                                                                                                                              <  1:189867‑M2/88‑1 (MQ=255)
                                                                                       gataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAAttt                                                                                                                                                           >  1:360637‑M2/161‑251 (MQ=255)
                                                                                             gctcgggttaggtgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTGAAGGATATCTGAAGGTATATTCAGAATTTGAATa                                                                                                                                                      <  2:198302‑M2/96‑1 (MQ=255)
                                                                                                 gggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATg                                                                                                                                                 <  1:157007‑M2/101‑1 (MQ=255)
                                                                                                      cgttgtcaccccggcactggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCATAATTTGCATAAAATGcagac                                                                                                                                            >  2:128635‑M2/146‑251 (MQ=255)
                                                                                                      agttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGcagac                                                                                                                                            <  1:324780‑M2/106‑1 (MQ=255)
                                                                                                                          ggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAAtt                                                                                                                                                                                            >  2:285865‑M2/126‑183 (MQ=255)
                                                                                                                          ggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAAtt                                                                                                                                                                                            <  1:285865‑M2/58‑1 (MQ=255)
                                                                                                                                   gtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGtt                                                                                                                  >  1:33888‑M2/117‑248 (MQ=255)
                                                                                                                                           cgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGt                                                                                                       >  2:12857‑M2/109‑251 (MQ=255)
                                                                                                                                            gttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGt                                                                                                       >  1:98700‑M2/108‑250 (MQ=255)
                                                                                                                                                    tctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAg                                                                                                <  1:395862‑M2/150‑1 (MQ=255)
                                                                                                                                                       gcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAAc                                                                                           <  1:358066‑M2/155‑1 (MQ=255)
                                                                                                                                                                 cggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCa                                                                                  >  1:73646‑M2/87‑250 (MQ=255)
                                                                                                                                                                        cagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACg                                                                          <  1:318214‑M2/172‑1 (MQ=255)
                                                                                                                                                                                  tccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGc                                                                 >  1:109006‑M2/70‑250 (MQ=255)
                                                                                                                                                                                       gaagtttccttttactgttccagttaatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGGACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTACTTAAGGATATCTGAAGGTATATTCAGAATTTGAGTAAAATGCAGACAGAACTATATGGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCATCGCGCAAAAGCACGAGGGTTTacagagg                                                            >  2:202702‑M2/65‑248 (MQ=255)
                                                                                                                                                                                        atgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAagag                                                          <  2:33888‑M2/188‑1 (MQ=255)
                                                                                                                                                                                                ttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGAtt                                                   >  1:62385‑M2/56‑250 (MQ=255)
                                                                                                                                                                                                   tacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGcc                                               >  1:195007‑M2/53‑251 (MQ=255)
                                                                                                                                                                                                             gttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGc                                     <  1:113773‑M2/209‑1 (MQ=255)
                                                                                                                                                                                                                         gcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGc                          >  1:32671‑M2/31‑250 (MQ=255)
                                                                                                                                                                                                                          caggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCa                         >  1:132090‑M2/30‑250 (MQ=255)
                                                                                                                                                                                                                                   gatagcgtcggaagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAAAAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCaa               <  2:98700‑M2/231‑1 (MQ=255)
                                                                                                                                                                                                                                      agcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAg            <  1:12857‑M2/234‑1 (MQ=255)
                                                                                                                                                                                                                                        cgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAATATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTg          <  2:124956‑M2/236‑1 (MQ=255)
                                                                                                                                                                                                                                           cggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGAcc       >  1:382834‑M2/13‑251 (MQ=255)
                                                                                                                                                                                                                                                gaggctgTGTGAAATGCATAGTAAGAACCGGTCTTGAAAAAATATCGCCGAATGTAGCGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTGGAATAAAATGCAGACAGAAATATATGGAAAAGGAGGGTGTTAGAACAGAAGTATTTCAGAAACCCCTCGCGCAAAAGCACGGGGGTTGGGGGAAGAGGAAGATTAGCCGGTAGTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTaa  <  2:79858‑M2/244‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   
ACAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTAGAGCAGTGCCCAACCCAGGGCTTTCCAGCTAAAAAGAGGTTTAGCCTGGACTTCTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAA  >  NC_000913/4150007‑4150497

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