Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 4,235,905 16376 bp→AAGTAAAGGAAAACATCAACGTTCGCAAAAGAAAATGCCGATTTGCGGCCGC yjbEmalM yjbE, yjbF, yjbG, yjbH, yjbT, psiE, xylE, malH, malG, malF, malE, malK, lamB, malM

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 4235905 4252280 16376 47 [0] [0] 49 yjbE–malM yjbE,yjbF,yjbG,yjbH,yjbT,psiE,xylE,malH,malG,malF,malE,malK,lamB,malM

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 42359040 (0.000)39 (0.760)
+AAGTAAAGGAAAACATCAACGTTCGCAAAAGAAAATGCCGATTTGCGGCCGC
35/374 0.1 100% intergenic (+497/‑2) pgi/yjbE glucose‑6‑phosphate isomerase/uncharacterized protein YjbE
?NC_000913 4252281 = 0 (0.000)intergenic (+1807/‑225) malM/ubiC maltose regulon periplasmic protein/chorismate lyase

TTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAAAAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAA  >  NC_000913/4235681‑4236151
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                       
ttCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgc                                                                                                                                                                                                                               >  1:57029‑M1/1‑224 (MQ=255)
     aaaaTACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaa                                                                                                                                                                                                                            >  1:228083‑M1/1‑219 (MQ=255)
        aTACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaa                                                                                                                                                                                                                       >  1:27386‑M1/1‑216 (MQ=255)
                   aaTTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgc                                                                                                                                                                                                           <  1:474474‑M1/251‑47 (MQ=255)
                                   attttaAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaagg                                                                                                                                                                                           <  1:309250‑M1/251‑63 (MQ=255)
                                       taAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgc                                                                                                                                                                                                                   >  2:205609‑M1/1‑185 (MQ=255)
                                       taAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgc                                                                                                                                                                                                                   <  1:205609‑M1/223‑39 (MQ=255)
                                           tttATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCACGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacat                                                                                                                                                                                    >  2:265578‑M1/1‑181 (MQ=255)
                                                       aTTGTTATTTAAATATAAGATGAAACTTATATGTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaccatcaccgttcctggc                                                                                                                                                                       >  1:354984‑M1/1‑169 (MQ=255)
                                                          gTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatc                                                                                                                                                                    <  1:319106‑M1/251‑86 (MQ=255)
                                                          gTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAAGTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatc                                                                                                                                                                    <  2:273469‑M1/251‑86 (MQ=255)
                                                                tAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatcctggac                                                                                                                                                              >  1:400372‑M1/1‑160 (MQ=255)
                                                                                cTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatcctggacggtgatgcccctacgg                                                                                                                                              >  1:158692‑M1/1‑144 (MQ=255)
                                                                                                   ccAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaatgaaaacatcaccgttcctggcatcctggacggtgatgcccctacggttgccctcgccagcacggg                                                                                                                           <  1:58712‑M1/251‑127 (MQ=255)
                                                                                                   ccAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatcctggacggtgatgcccctacggttgccctcgccagcacggg                                                                                                                           >  2:458522‑M1/1‑125 (MQ=255)
                                                                                                           ccTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatcctggacggtgatgcccctacggttgccctcgccagcacgggcatcggt                                                                                                                    >  1:235636‑M1/1‑117 (MQ=255)
                                                                                                                aaCGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatcctggacggtgatgcccctacggntgccctcgccagcacgggcatcggtaaagcg                                                                                                              >  2:244256‑M1/1‑112 (MQ=255)
                                                                                                                          cTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatcctggacggtgatgcccctacggttgccctcgccagcacgggcatcggtaaagcgtaaggttca                                                                                                     >  1:307147‑M1/1‑102 (MQ=255)
                                                                                                                               tCCCCGCGCCGCAATTCGCATGCTTTAGTTGTGTATACTCGATCCTGCCCGAAATGTTTTTGGGTAAATCTCCCTTCATTCAATGAAGGGAAATTGTcagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgagttgcggccgctacggaaaggaaaacatcaccgttcctggcatcctggacggtgatgcccctacggttgccctcgccagcacgggcatcggaaaagcgtaaggttccacctct                                                                                               >  2:20610‑M1/1‑97 (MQ=255)
                                                                                                                                    gcgcCGTAATTCTCATGCTTTAGTTGTGTATACTCGAGCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatcctggacggtgatgcccctacggttgccctcgccagcacgggcatcggtaaagcgtaaggttcaacatcgtttta                                                                                          <  2:286093‑M1/251‑160 (MQ=255)
                                                                                                                                    gcgcCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaaga                                                                                                                                                                                                                         <  2:223046‑M1/124‑33 (MQ=255)
                                                                                                                                    gcgcCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaaga                                                                                                                                                                                                                         >  1:223046‑M1/1‑92 (MQ=255)
                                                                                                                                       ccGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatcctggacggtgatgcccctacggttgccctcgccagcacgggcatcggtaaagcgtaaggttcaacatcgttttacca                                                                                       <  2:285267‑M1/251‑163 (MQ=255)
                                                                                                                                       ccGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatcctggacggtgatgcccctacggttgccctcgccagcacgggcatcggtaaagcgtaaggttcaacatcgttttacc                                                                                        >  2:399090‑M1/1‑89 (MQ=255)
                                                                                                                                                  aTGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggcagctaagtaaaggaaaacatcaccgttcctggcatcctggacggtgatgcccctacggttgccctcgccagcacgggcatcggtaaagcgtaaggttaaccatcgttttaccacttcttccgat                                                                            >  2:497495‑M1/1‑78 (MQ=255)
                                                                                                                                                          gTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatcctggacggtgatgcccctacggttgccctcgccagcacgggcatcggtaaagcgtaaggttcaacatcgttttaccacttcatgcgattgttgcgt                                                                    <  1:463489‑M1/251‑182 (MQ=255)
                                                                                                                                                          gTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatcctggacggtgatgcccctacggttgccctcgccagcacgggcatcggtaaagcggaaggttcaacatcgttttaccacttcatgcgattgttgcgt                                                                    >  1:349137‑M1/1‑70 (MQ=255)
                                                                                                                                                                    cTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatcctggacggtgatgcccctacggttgccctcgccagcacgggcatcggtaaagcgtaaggttcaacatcgttttaccacttcatgcgattgttgcgtttttgttgcg                                                          >  1:97262‑M1/1‑60 (MQ=255)
                                                                                                                                                                    cTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatcctggacggtgatgcccctacggttgccctcgccagcacgggcatcggtaaagcgt                                                                                                             <  1:376218‑M1/200‑141 (MQ=255)
                                                                                                                                                                    cTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatcctggacggtgatgcccctacggttgccctcgccagcacgggcatcggtaaagcga                                                                                                             >  2:376218‑M1/1‑60 (MQ=255)
                                                                                                                                                                       gATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatcctggacggtgatgcccctacggttgccctcgccagcacgggcatcgg                                                                                                                     <  2:447238‑M1/189‑133 (MQ=255)
                                                                                                                                                                       gATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatcctggacggtgatgcccctacggttgccctcgccagcacgggcatcgg                                                                                                                     >  1:447238‑M1/1‑57 (MQ=255)
                                                                                                                                                                        aTCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatcctggacggtgatgcccctacggttgccctcgccagcacgggcatcggtaaagcgtaaggttcaacatcgttttaccacttcatgcgattgttgcgtttttgttgcgtatt                                                      <  2:155857‑M1/251‑196 (MQ=255)
                                                                                                                                                                          cccgcccgAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatcctggacggtgatgcccctacggttgccctcgccagcacgggcatcggtaaagcgtaaggttcaacatcgttttaccacttcatgcgattgttgcgtttttgttgcgtagtag                                                    <  2:106752‑M1/251‑198 (MQ=255)
                                                                                                                                                                               ccgAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatcctggacggtgatgcccctacggttgccctcgccagcacgggcatcggtaaagcgtaaggttcaacatcgttttaccacttcatgcgattgttgcgtttttgttgcgtattagatcac                                               <  1:319993‑M1/251‑203 (MQ=255)
                                                                                                                                                                                      gTTTTTGGGTAAATCTCCATTCATTCACTGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatcctggacggtgatgcccctacggttgccctcgccagcacgggcatcggtaaagcgtaaggttcaacatcgttttaccacttcatgcgattgttgcgtttttgttgcgtattagatcacttaattt                                        >  2:102498‑M1/1‑42 (MQ=255)
                                                                                                                                                                                                   tctcCATTCATTCAATGAAGGGAAATTGTaagtaaaggaaaacctcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatcctggacggtgatgcccctacggttgccctcgccagcacgggcatcggtaaagcgtaaggttcaacatcgttttaccacttcatgcgattgttgcgtttttgttgcgtattagatcacttaatttgctttacatctcc                           >  2:189168‑M1/1‑29 (MQ=255)
                                                                                                                                                                                                       cattcattCAATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatcctggacggtgatgcccctacggttgccctcgccagcacgggcatcggtaaagcgtaaggttcaacatcgttttaccacttcatgcgattgttgcgtttttgttgcgtattagatcacttaatttgctttacatctcccgt                        >  1:209535‑M1/1‑25 (MQ=255)
                                                                                                                                                                                                        attcattcAATGAAGGGACATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatcctggacggtgatgcccctacggttgccctcgccagcacgggcatcggtaaagcgtaaggttcaacatcgttttaccacttcatgcgattgttgcgtttttgttgcgtattagatcacttaatttgctttacatctcccgta                       >  2:409488‑M1/1‑24 (MQ=255)
                                                                                                                                                                                                           cattcaATGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatcctggacggtgatgcccctacggttgccctcgccagcacgggcatcggtaaagcgtaaggttcaacatcgttttaccacttcatgcgattgttgcgtttttgttgcgtattagatcacttaatttgctttacatctcccgtaaac                    >  2:82030‑M1/1‑21 (MQ=255)
                                                                                                                                                                                                                aaTGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatcctggacggtgatgcccctacggttgccctcgccagcacgggcatcggtaaagcgtaaggttcaacatcgttttaccactttatgcgattgttgcgtttttgttgcgtattagatcacttaatttgctttacatcgcccgtaaacactttt              <  2:57029‑M1/251‑236 (MQ=255)
                                                                                                                                                                                                                aaTGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatcctggacggtgatgcccctacggttgccctcgccagcacgggcatcggtaaagcgtaaggttcaacatcgttttaccacttcatgcgattgttgcgtttttgttgcgtattagatcacttaatttgctttacatctcccgtaaacactttt              <  1:458522‑M1/251‑236 (MQ=255)
                                                                                                                                                                                                                 aTGAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgaggccgctaagtaacggcaaaaatcaccgttcctggcatcctggacggtgatgcccctacggttgccctcgccagcacgagcatcggtaaagcgtaaggttcaacatcgttttaccacttcatgcgcttgttgcgtttttgttgcgtattagataacttaatttgctttacatctcccgtaaacacttttc             >  2:156890‑M1/1‑15 (MQ=255)
                                                                                                                                                                                                                   gAAGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgagttgcggccgctaagtaaatgaaaacatcaccgttcctggcatcctggacggtgatgcctctacggttgccctcgccagcacgggcatcggtaaagcgtaaggttcaacatcgttttaccacttcatgcgattgttgcgtttttgttgcgtattagatcccttaatttgctttacatctcccgtaaacacttttcg            <  2:341505‑M1/250‑238 (MQ=255)
                                                                                                                                                                                                                    aaGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatcctggacggtgatgcccctacggttgccctcgccagcacgggcatcggtaaagcgtaaggttcaacatcgttttaccacttcatgcgattgttgcgtttttgttgcgtattagatcacttaatttgctttacatctcccgtaaacacttttctgc          <  1:111804‑M1/251‑240 (MQ=255)
                                                                                                                                                                                                                     aGGGAAATTGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgctaagtaaaggaaaacatcaccgttcctggcatcctggacggtgatgcccctacggttgccctcgccagcacgggcatcggtaaagcgtaaggttcaacatcgttttaccacttcatgcgattgttgcgtttttgttgcgtattagatcacttaatttgctttacatctcccgtaaacacttttctgcg         <  2:97262‑M1/251‑241 (MQ=255)
                                                                                                                                                                                                                            ttGTaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggacgctaagtaaaggaaaacatcaccgttcctggcatcctggacggagatgcccctacggttgccctcgccagcacgggcatcggtaaagcgtaaggttcaacatcgttttaccacttcatgcgattgttgcgtttttgttgcgtattagatcacttaatttgctttacatctcccgtaaacacttttctgcgatacaat  >  2:264168‑M1/1‑4 (MQ=255)
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                       
TTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAATTGTTATGAAAAAAGTTCTGTATGGCATTTTTGCCATATCTGCGCTTGCGGCGACTTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAA  >  NC_000913/4235681‑4236151

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