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 57 [0] [0] 56 [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)56 (0.960)
+TG
54/490 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

CAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTAGAGCAGTGCCCAACCCAGGGCTTTCCAGCTAAAAAGAGGTTTAGCCTGGACTTCTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACG  >  NC_000913/4150008‑4150499
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
atggtctgagcaaacagaatgccctcatcaccaatcgctttttgcagacgcggcagcacttcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTgg                                                                                                                                                                                                                                                   >  2:121981‑M2/247‑251 (MQ=255)
     ctgagcaaacagaatgccctcatcaccaatcgctttttgcagacgcggcagcacttcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGc                                                                                                                                                                                                                                              >  2:178208‑M2/242‑251 (MQ=255)
          caaacagaatgccctcatcaccaatcgctttttgcagacgcggcagcacttcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGt                                                                                                                                                                                                                                         >  2:62690‑M2/237‑251 (MQ=255)
            aacagaatgccctcatcaccaatcgctttttgcagacgcggcagcacttcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTa                                                                                                                                                                                                                                        >  1:37049‑M2/235‑250 (MQ=255)
                gaatgccctcatcaccaatcgctttttgcagacgcggcagcacttcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAAc                                                                                                                                                                                                                                   >  1:329394‑M2/231‑251 (MQ=255)
                             cccaatcgctttttgcagacgcggcagcacttcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGaaaaa                                                                                                                                                                                                                      <  1:375649‑M2/34‑1 (MQ=255)
                              ccaatcgcttttggcagacgcggcagcacttcccatatggactccttgctggcagcgataatgctcgggttagatgtcaccccggcaatggggaatatgcgtgccagccgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGGGGAATGCATAGTAAGAACCTGTCTTGcaaaa                                                                                                                                                                                                                      >  2:338124‑M2/217‑250 (MQ=255)
                                aatcgctttttgcagacgcggcagcacttcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggggtccagatacagttacatgatgtttccttttacagttccagttcatgttgcagcagtctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTTTCTTGACAAAat                                                                                                                                                                                                                    >  2:271964‑M2/215‑250 (MQ=255)
                                aatcgctttttgcagacgaggcagcacttcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAata                                                                                                                                                                                                                   <  1:338124‑M2/37‑1 (MQ=255)
                                 atcgctttttgcagacgcggcagcacttcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGCAAAAatat                                                                                                                                                                                                                  >  1:64033‑M2/214‑251 (MQ=255)
                                      tttttgcagacgcggcagcacttcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCg                                                                                                                                                                                                             >  2:46407‑M2/209‑251 (MQ=255)
                                            cagacgcggcagcacttcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTa                                                                                                                                                                                                       <  1:446771‑M2/49‑1 (MQ=255)
                                                    ggcagcacttcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgtagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAAt                                                                                                                                                                                               >  1:59775‑M2/195‑251 (MQ=255)
                                                         acttcccatatggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCtt                                                                                                                                                                                          >  1:320455‑M2/190‑251 (MQ=255)
                                                            tcccatatggactccttgctggcagcgataatgcgcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTaa                                                                                                                                                                                        <  1:368681‑M2/64‑1 (MQ=255)
                                                                   tggactccttgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCt                                                                                                                                                                                >  1:154022‑M2/180‑251 (MQ=255)
                                                                            tgctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTTAGGATATCTGAAGGtat                                                                                                                                                                        <  1:84920‑M2/80‑1 (MQ=255)
                                                                              ctggcagcgataatgctcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATAtt                                                                                                                                                                     >  1:179465‑M2/169‑251 (MQ=255)
                                                                                         aatgctcgggttagttgtcaccccggcaatggggaatatgtgtgccagacgttcgacttctgcgccgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGACTGCAACGACAATTCCTTAAGGATATCTGAAGGAATATTCAGAATTTGa                                                                                                                                                           >  2:368158‑M2/158‑250 (MQ=255)
                                                                                          attccctggttaattgtcaccccggcaatggggaatattcgttccagacgttcggcttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAATGTATATTCAGAATGTGAAt                                                                                                                                                         <  2:382463‑M2/95‑1 (MQ=255)
                                                                                              tcgggttagttgtcaccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATaaa                                                                                                                                                      >  2:194055‑M2/153‑250 (MQ=255)
                                                                                                           caccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAat                                                                                                                                        <  2:164311‑M2/112‑1 (MQ=255)
                                                                                                            accccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAat                                                                                                                                        >  1:12635‑M2/139‑250 (MQ=255)
                                                                                                             ccccggcaatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAatat                                                                                                                                      <  1:185695‑M2/114‑1 (MQ=255)
                                                                                                                    aatggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGaaa                                                                                                                               <  1:429583‑M2/121‑1 (MQ=255)
                                                                                                                     atggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGaaa                                                                                                                               >  1:169298‑M2/130‑250 (MQ=255)
                                                                                                                       ggggaatatgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACg                                                                                                                            <  1:121981‑M2/124‑1 (MQ=255)
                                                                                                                              atgcgtgccagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGgtg                                                                                                                      <  1:96157‑M2/130‑1 (MQ=255)
                                                                                                                                      cagacgttcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAg                                                                                                             <  2:329394‑M2/139‑1 (MQ=255)
                                                                                                                                             tcgacttctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTAt                                                                                                       <  1:368158‑M2/145‑1 (MQ=255)
                                                                                                                                                   tctgcgacgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCATAAGGATATCTGAAGGTATATTCAGAAGTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTAGTTCAg                                                                                                  <  2:156726‑M2/150‑1 (MQ=255)
                                                                                                                                                          cgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTc                                                                                         <  1:46407‑M2/159‑1 (MQ=255)
                                                                                                                                                          cgttagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTc                                                                                         <  1:178208‑M2/159‑1 (MQ=255)
                                                                                                                                                             tagcggtgtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTcgcg                                                                                      >  2:448145‑M2/90‑251 (MQ=255)
                                                                                                                                                                    gtccagatacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGc                                                                               <  2:59775‑M2/169‑1 (MQ=255)
                                                                                                                                                                          atacagttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAgg                                                                         >  1:309846‑M2/77‑251 (MQ=255)
                                                                                                                                                                               gttccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTg                                                                    <  1:212462‑M2/180‑1 (MQ=255)
                                                                                                                                                                                 tccatgatgtttccttttacagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGc                                                                   <  1:396699‑M2/181‑1 (MQ=255)
                                                                                                                                                                                                    cagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCgg                                               >  2:15594‑M2/51‑251 (MQ=255)
                                                                                                                                                                                                    cagttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCgg                                               >  1:189466‑M2/51‑251 (MQ=255)
                                                                                                                                                                                                      gttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGAt                                                      <  2:345852‑M2/194‑1 (MQ=255)
                                                                                                                                                                                                      gttccagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGAt                                                      >  1:345852‑M2/49‑242 (MQ=255)
                                                                                                                                                                                                          cagttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTAc                                         <  1:24719‑M2/207‑1 (MQ=255)
                                                                                                                                                                                                             ttcatgttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGCCAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAACGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGc                                       >  2:438762‑M2/42‑250 (MQ=255)
                                                                                                                                                                                                                  gttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCgg                                               <  1:379491‑M2/201‑1 (MQ=255)
                                                                                                                                                                                                                  gttgcagcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCgg                                               >  2:379491‑M2/37‑237 (MQ=255)
                                                                                                                                                                                                                       agcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATCCCTGCCGc                            <  1:40167‑M2/220‑1 (MQ=255)
                                                                                                                                                                                                                       agcaggctggcgatagcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGc                            >  1:76876‑M2/32‑251 (MQ=255)
                                                                                                                                                                                                                                   atagcggcggcagagggtgGTTTGAAGGCATAGTAAGAACCTTTCTTTAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCTCAATCCCGGCAAt                <  1:370685‑M2/232‑1 (MQ=255)
                                                                                                                                                                                                                                     agcgtcggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGATTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAAt                >  2:87219‑M2/18‑249 (MQ=255)
                                                                                                                                                                                                                                     agcgtcggcagaggcggTGGGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTAAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAACCCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATCCCTGCCGCAATCCCGGCAATAg              <  2:385630‑M2/234‑1 (MQ=255)
                                                                                                                                                                                                                                          cggcagagtctgTGTGGAATGCAGATTAAGACCCTGTCTTGAAAAAATATCGCCGAATGTAACGCCAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGAAAAAGCTCGAGGGTTGGCAGAAGAGGAAGATTAGCCGGTATTCCGCATTCCTGCCGCAATCCCGGCAATAGTGAc          <  2:303107‑M2/238‑1 (MQ=255)
                                                                                                                                                                                                                                           ggcagaggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGAcc         >  1:397985‑M2/12‑250 (MQ=255)
                                                                                                                                                                                                                                            gcagaggctgTTTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAACCCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCAt       <  1:385903‑M2/241‑1 (MQ=255)
                                                                                                                                                                                                                                                aggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAAc   >  1:301153‑M2/7‑251 (MQ=255)
                                                                                                                                                                                                                                                 ggctgTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACg  >  1:98450‑M2/6‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
CAGCGCCGCCAGAGAGGCCGCCCATAGCACCACGCCGATTACTGCGGCAATAATTTTAATCCGCTTCGGAAAGAGGAATACCGGGATCAACCACAGCGAACTGAATAACAGCGAGTCGCGAATGCCGTTAGTGCCACTATAACCACTGATGTAAATAATGGCCTGTAGCAGAGTAGAGAAAAACCAAAAGTAGAGCAGTGCCCAACCCAGGGCTTTCCAGCTAAAAAGAGGTTTAGCCTGGACTTCTGTGGAATGCATAGTAAGAACCTGTCTTGAAAAAATATCGCCGAATGTAACGACAATTCCTTAAGGATATCTGAAGGTATATTCAGAATTTGAATAAAATGCAGACAGAAATATATTGAAAACGAGGGTGTTAGAACAGAAGTATTTCAGAAAACCCTCGCGCAAAAGCACGAGGGTTTGCAGAAGAGGAAGATTAGCCGGTATTACGCATACCTGCCGCAATCCCGGCAATAGTGACCATTAACG  >  NC_000913/4150008‑4150499

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