Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 2,498,768 9486 bp→CCC [pyrS]fryB [pyrS], pyrR, ypdC, fryA, ypdE, ypdF, fryC, fryB

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 2498768 2508253 9486 41 [0] [0] 41 [pyrS]–fryB [pyrS],pyrR,ypdC,fryA,ypdE,ypdF,fryC,fryB

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 24987670 (0.000)40 (0.620)
+CCC
37/484 0.4 100% coding (97/1698 nt) pyrS sensor histidine kinase PyrS
?NC_000913 2508254 = 0 (0.000)intergenic (‑12/+207) fryB/glk putative PTS enzyme IIB component FryB/glucokinase

AACCAGAAGTGCATTATCAAGAGATGCATTGCTAATACACCATTTATCTATCAATCCCCTTTTAAAATTGTGGTTTTTCCTCAAAAACGGCCTGCTGGTCATCAGCCTTCAGGTTACCTATCATAGAGGTTTAATCCTTATTCAGAGTCACCCGTGCACGAAATATTCAACATGCTGCTGGCGGTCTTCGATCGGGCAGCGTTAATGCTTATCTGCCTGTTCTTTCTCATCCGTATCCGCCTGTTTCGCGAACTGTTGCACAAGTCGGCGCACTCCCCAAAAGAATTGCTCGCCGTTACCGCCATTTTCTCGCTGTTCGCCCTGTTCAGCACCTGGTCCGGCGTTCCCGTAGAAGGCTCGCTGGTGAACGTACGTATTATCGCGGTGATGTCCGGCGGGATTCTGTTTGGCCCGTGGGTAGGCATCATTACCGGCGTGATTGCGGGTATTCACCGGTATTTAATTGATATCGGCGGCGTGACGGCGATCCCCT  >  NC_000913/2498518‑2499010
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                   
aaCCAGAAGTGCATTATCAAGAGATGCATTGCTAATACACCATTTATCTATCAATCCCCTTTTAAAATTGTGGTTTTTCCTCAAAAACGGCCTGCTGGTCATCAGCCTTCAGGTTACCTATCATAGAGGTTTAATCCTTATTCAGAGTCACCCGTGCACGAAATATTCAACATGCTGCTGGCGGTCTTCGATCGGGCAGCGTTAATGCTTATCTGCCTGTTCTTTCTCCTCCGTATCCGCCTGTTTCGCGc                                                                                                                                                                                                                                                    >  1:175201‑M1/1‑250 (MQ=255)
          gCATTATCAAGAGATGCATTGCTAATACACCATTTATCTATCAATCCCCTTTTAAAATTGTGGTTTTTCCTCAAAAACGGCCTGCTGGTCATCAGCCTTCAGGTTACCTATCATAGAGGTTTAATCCTTATTCAGAGTCACCCGTGCACGAAATATTCAACATGCTGCTGGCGGTCTTCGATCGGGCAGCGTTAATGCTTATCTGCCTGTTCTTTCTCATCCGTATCCGCCTGTTTCGCGccccgttcagg                                                                                                                                                                                                                                          >  1:428527‑M1/1‑240 (MQ=255)
               aTCAAGAGATGCATTGCTAATACACCATTTATCTATCAATCCCCTTTTAAAATTGTGGTTTTTCCTCAAAAACGGCCTGCTGGTCATCAGCCTTCAGGTTACCTATCATAGAGGTTTAATCCTTATTCAGAGTCACCCGTGCACGAAATATTCAACATGCTGCTGGCGGTCTTCGATCGGGCAGCGTTAATGCTTATCTTCCTGTTCTTTCTCATCCGTATCCGCCTGTTTCGCGccccggtcagtgat                                                                                                                                                                                                                                       >  2:531305‑M1/1‑235 (MQ=255)
                    gagaTGCATTGCTAATACACCATTTATCTAGCAATCCCCTTTTAAAATTGTGGTTTTTCCTCAAAAACGGCCTGCTGGTCATCAGCCTTCAGGTTACCTATCATAGAGGTTTAATCCTTATTCAGAGGCCCCCGTGCACGAAATATTCAACATGCTGTTGGCGGTCTTCGATCGGGCAGCGTTAATGATTATCTGCCTGTTCTTTCTCATCCGTATCCGCCTGTTTCGCGccccgttcagtgatagttgac                                                                                                                                                                                                                                <  2:80249‑M1/251‑22 (MQ=255)
                             tGCTAATGCCCCATTTATCTATCAATCCCCTTTTAAAATTGTGGTTTTTCCTCAGAAACGGCCTGCTGGTCAACAGCCTTCAGGTTCCCTATGATAGAGGTTTAATTCTTATGCAGAGTCGTCCTTTCCCGAAATATTCAACATGCTGCGGGCGGTCTTGGAGCGGGCAGCGTTAATGCTTATCTGCCTGTTCTTTCTCATCCGTTGCCGCCTGTTTCGCGccccgttcagtgatagttgacgagtgaaa                                                                                                                                                                                                                        <  2:326406‑M1/250‑30 (MQ=255)
                                         aTTTATCTATCAATCCCCTTTTAAAATTGTGGTTTTTCCTCAAAAACGGCCTGCTGGTCATCAGCCTTCAGGTTACCTATCATAGAGGTTTAATCCTTATTCAGAGTCACCCGTGCACGAAATATTCAACATGCTGCTGGCGGTCTTCGATCGGGAAGAGTTAATGCTTTTCTGCCTGTTCTTTCTCATCCGTATCCGCCTGTTTCTCGccccgttcagtgatagttgacgagtgaaattatgaacgtgc                                                                                                                                                                                                            >  1:367905‑M1/1‑209 (MQ=255)
                                            tatctatcAATCCCCTTTTAAAATTGTGGTTTTTCCTCAAAAACGGCCTGCTGGTCATCAGCCTTCAGGTTACCTATCATAGAGGTTTAATCCTTATTCAGAGTCACCCGTGCACGAAATATTCAACATGCTGCTGGCGGGCTTCGATCGGGCAGCGTTAATGCTTTACTGCCTGTTTTTTCTCATCCGTATCCCCCTTTTTCGCGccccgttcagtgataggtgacgagtgaaattatgaacgggcaaa                                                                                                                                                                                                         >  2:489284‑M1/1‑206 (MQ=255)
                                            taactatcaatCCCCTTTTAAAATTGTGGTTTTTCCTCAAAAACGGCCTGCTGGTCATCATCTTTCAGGTTCCCTATCATAGAGGTTTAATCCTTATTCAGAGTCACCCGTGCACGAAATATTCAACATGCTGCTGGCGGTCTTCGATCGGGCAGCGTTAATGCTTATCTGCCTGTTCTTTCTCAACCGTAGCCGCCTGGTTCGCGccccggtcagtgatagttgacgagtgaaattatgaacgtgcaaa                                                                                                                                                                                                         <  2:209267‑M1/247‑45 (MQ=255)
                                             atctatcAATCCCCTTTTAAAATTGTGGTTTTTCCTCAAAAACGGCCTGCTGGTCATCAGCCTTCAGGTTACCTATCATAGAGGTTTAATCCTTATTCAGAGTCACCCGTGCACGAAATATTCAACATGCTGCTGGCGGTCTTCGATCGGGCAGCGTTAATGCTTATCTGCCTGTTCTTTCTCATCCGTATCCGCCTGTTTCGCGccccgttcagtgatagtt                                                                                                                                                                                                                                   >  1:355493‑M1/1‑205 (MQ=255)
                                             atctatcAATCCCCTTTTAAAATTGTGGTTTTTCCTCAAAAACGGCCTGCTGGTCATCAGCCTTCAGGTTACCTATCATAGAGGTTTAATCCTTATTCAGAGTCACCCGGGCACGAAATATTCAACATGCTGCTGGCGGTCTTCGATCGGGCAGCGTTAATGCTTATCTGCCTGTTCTTTCTCATCCGTATCCGCCTGTTTCGCGccccgttcagtgatagtt                                                                                                                                                                                                                                   <  2:355493‑M1/223‑19 (MQ=255)
                                              tctatcAATCCCCTTTTAAAATTGTGGTTTTTCCTCAAAAACGGCCTGCTGGTCATCAGCCTTCAGGTTACCTATCATAGAGGTTTAATCCTTATTCAGAGTCACCCGTGCACGAAATATTCAACATGCTGCTGGCGGTCTTCGATCGGGCAGCGTTAATGCTTATCTGCCTGTTCTTTCTCATCCGTATCCGCCTGTTTCGCGccccgttcagtgatagttgacgagtgaaattatgaacgtgcaaaac                                                                                                                                                                                                       >  1:45868‑M1/1‑204 (MQ=255)
                                                            tttAAAATTGTGGTTTTTCCTCAAAAACGGCCTGCTGGTCATCAGCCTTCAGGTTACCTATCATAGAGGTTTAATCCTTATTCAGATTCACCCGTGCACGAAATATTCAACATGCTGCTGGCGGTCTTCGTTCGGGCAGCGTTAATGCTTATCTGCCTGTTCTTTCTCATCCGTATCCGCCTGTTTCGCGccccgttcagtgatagttgacgagtgaaattatgaacgtgcaaaacaaatcgccgttac                                                                                                                                                                                          <  1:228887‑M1/249‑60 (MQ=255)
                                                                 aatttgggtttttcCTCAAAAACGGCCTGCTGGTCATCAGCCTTCAGGTTACCTATCATAGAGGTTTAATCCTTATTCAGAGTCACCCGTGCACGAAATATTCAACATGCTGCTGGCGGTCTTCGATCGGGCAGCGTTAATGCTTATCTGCCTGTTCTTTCTCATCCGTATCCGCCTGTTTCGCGccccgttcagtgatagttgacgagtgaaattatgaacgtgcaaaacaaatcgccgttaccagacaa                                                                                                                                                                                   <  2:513882‑M1/245‑67 (MQ=255)
                                                                  aTTGTTGTTTTTCCTCAAAAACGGCCTGCTGGTCATCAGCCTTCAGGTTACCTATCATAGAGGTTTAATCCTTATTCAGAGTCACCCGTGCACGAAATATTCAACATGCTGCTGGCGGTCTTCGATCGGGCAGCGTTAATGCTTATCTGCCTGTTCTTTCTCATCCGTATCCGCCTGTTTCGCGccccgttcagtgatagttgacgagtgaaattatgaacgtgcaaaacaaatcgccgttaccagacaa                                                                                                                                                                                   >  1:432813‑M1/1‑184 (MQ=255)
                                                                              ccTCAAAAACGGCCTGCTGGTCATCAGCCTTCAGGTTACCTATCATAGAGGTTTAATCCTTATTCAGAGTCACCCGTGCACGAAATATTCAACATGCTGCTGGCGGTCTTCGATCGGGCAGCGTTAATGCTTATCTGCCTGTTCTTTCTCATCCGTATCCGCCTGTTTCGCGccccgttcagtgatagttgacgagtgaaattatgaacgtgcaaaacaaatcgccgttaccagacaaataatgcattta                                                                                                                                                                       <  1:427085‑M1/250‑79 (MQ=255)
                                                                                             gctgGTCATCAGCCTTCAGGTTACCTATCATAGAGGTTTAATCCTTATTCAGAGTCACCCGTGCACGAAATATTCAACATGCTGCTGGCGGTCTTCGATCGGGCAGCGTTAATGCTTATCTGCCTGTTCTTTCTCATCCGTATCCGCCTGTTTCGCGccccgttcagtgatagttgacgagtgaaattatgaacgtgcaaaacaaatcgccgttaccagacaaataatgcatttactggagaattaatc                                                                                                                                                         <  1:49022‑M1/249‑93 (MQ=255)
                                                                                                                            agagGTTTAATCCTTATTCAGAGTCACCCGTGCACGAAATATTCAACATGCTGCTGGCGGTCTTCGATCGGGCAGCGTTAATGCTTATCGGCCTGTTCTTTCTCATCCGTATCCGCCTGTTTCGCGccccgttcagtgatagttgacgagtgaaattatgaacgtgcaaaacaaatcgccgttaccagacaaataatgcatttactggagaattaatcccgctgtacggaagtgtgacgcaggtcgacctt                                                                                                                        <  2:419865‑M1/251‑126 (MQ=255)
                                                                                                                                          tATTCAGAGTCACCCGTGCACGAAATATTCAACATGCTGCTGGCGGTCTTCGATCGGGCAGCGTTAATGCTTATCTGCCTGTTCTTTCTCATCCGTATCCGCCTGTTTCGCGccccgttcagtgatagttgacgagtgaaatta                                                                                                                                                                                                                     <  1:327522‑M1/144‑33 (MQ=255)
                                                                                                                                          tATTCAGAGTCACCCGTGCACGAAATATTCAACATGCTGCTGGCGGTCTTCGATCGGGCAGCGGTAATGCTTATCTGCCTGTTCTTTATCATCCGTATCCGCCTGTTTCGCGccccgttcagtgatagttgccgagtgaaatta                                                                                                                                                                                                                     >  2:327522‑M1/1‑112 (MQ=255)
                                                                                                                                           aTTCAGAGTACCCCGTCCACGAAATATTCTACATTCTGCTGGCGTTCTTCGATCGGGCAGCGTTAATGCTTATCTGCCTGTTCTTTCTCATCCGTATCCGCCTGTTTCGCGccccgttcagtgatagttgacgagtgaaattatgaacgtgcaaaacaaatcgccgttaccagacaaataatgcatttactttagaattaatcccgctgtacggaagtgtgacgcatgtcgaccttgtgaggggacaaaa                                                                                                          <  2:508496‑M1/250‑140 (MQ=255)
                                                                                                                                             tCAGAGTCACCCGTGCACGAAATATTCAACATGCTGCTGGCGGTCTTCGATCGGGCAGCGTTAATGCTTATCTGCCTGTTCTTTCTCATCCGTATCCGCCTGTTTCGCGccccgttcagtgatagttgacgagtgaaattatgaacgtgcaaaacaaatcgccgttaccagacaaataatgcatttactggagaattaatcccgctgtacggaagtgtgacgcaggtcgaccttgtgaggggacacaaag                                                                                                        <  1:89541‑M1/250‑142 (MQ=255)
                                                                                                                                               agagCCACCCGTGCACGAAATATTCAACATGCTGCTGGCGGTCTTCGTTCGGGCAGCGTTAATGCTTATCTGCCTGTTCTTTCTCATCCGTATCCGCCTGTTTCGCGccccgttcagtgatagttgacgagtgaaattatgaacgtgcaaaacaaatcgccgttaccagacaaataatgcatttactggagaattaatcccgctgtacgtaagtgtgacgcaggtcgaccttgtgaggggacacaaagaaa                                                                                                     <  1:37865‑M1/251‑145 (MQ=255)
                                                                                                                                                gagTCACCCGTGCACGAAATATTCAACATGCTGCTGGCGGTCTTCGATCGGGCAGCGTTAATGCTTATCTGCCTGTTCTTTCTCATCCGTATCCGCCTGTTTCGCGccccgttcagtgctagttgacgagtgaaattatgaacgtgcaaaacaaatcgccgttaccagacaaataatgcatttactggagaattaatcccgctgtacggaagtgtgacgcaggtcgaccttgtgaggggacacaaagaaa                                                                                                     >  1:140997‑M1/1‑106 (MQ=255)
                                                                                                                                                 agTCACCCGTGCACGAAATATTCAACATGCTGCTGGCGGTCTTCGATCGGGCAGCGTTAATGCTTTTCTGCCTGTTATTTCTCATCCGTATCCGCCTGTTTCGCGacccgttcagcgatagtttacgagtgaaattatgaacgtgcaaaacaaatcgccgttaccagacaaataaagcatttactggagaattaatcccgctgtacggaagtgtgacgcaggtcggccttgtgaggggacacaaagaaaa                                                                                                    >  2:150054‑M1/1‑105 (MQ=255)
                                                                                                                                                                  atatTCAACATGCTGCTGGCGGTCTTCGATCGGGCAGCGTTAATGCTTATCTGCCTGTTCTTTCTCATCCGTATCCGCCTGTTTCGCGccccgttcagtgatagttgacgagtgaaattatgaacgtgcaaaacaaatcgccgttaccagacaaataatgcatttactggagaattaatcccgctgtacggaagtgtgacgcaggtcgaccttgtgaggggacacaaagaaaaattgaataaactgtatga                                                                                  <  1:509531‑M1/251‑164 (MQ=255)
                                                                                                                                                                       cAACATGCTGCTGGCGGTCTTCGATCGGGCAGCGTTAATGCTTATCTGCCTGTTCTTTCTCATCCGTATCCGCCTGTTTCGCGccccgttcagtgatagttgacgagtgaaattatgaacgtgcaaaacaaatcgccgttaccagacaaataatgcatttactggagaattaatcccgctgtacggaagtgtgacgcaggtcgaccttgtgaggggacacaaagaaaaattgaataaactgtatgattta                                                                              <  1:232235‑M1/250‑168 (MQ=255)
                                                                                                                                                                         aCATGCTTCTTGTGGTCTTCGATCGGGCAGCGTTAATGCTTATCGGCCTGTTCTTTCTCATCCGTATCCGCCTGGTTCGCGccccgttcagtgatagttgacgagtgaaattatgaacgtgcaaaacaaatcgccgttaccagacaaataatgcatttactggagaattaagcccgctgtacggaagtgtgacgcaggtcgaccttgtgaggggacacaaagaaaaattgaagaaactgtatgatttaaaa                                                                           <  2:460101‑M1/251‑171 (MQ=255)
                                                                                                                                                                            tgctgcTGGCGGTCTTCGATCGGGCAGCGTTAATGCTTATCTGCCTGTTCTTTCTCATCCGTATCCGCCTGTTTCGCGccccgttcagtgatagttgacgagtgaaattatgaacgtgcaaaacaaatcgccgttaccagacaaataatgcatttactggagaattaatcacgctgtacggaagtgtgacgcaggtcgaccttgtgaggggacacaaagaaaaattgaataaactgtatgctttaaa                                                                            >  1:33020‑M1/1‑78 (MQ=255)
                                                                                                                                                                             gctgctGTCGTTCTTCGATCGGTCAGCGTTAATGCTTATCTGCCTGTTCTTTCTCATCCGTATCCTCCTGTTTCGCGccccgttcagtgatagttgacgagtaaaattatgaacgtgcaaaacaaatcgccgttaccagacaaataatgcatttactggagaattaatcccgctgtacggaagtgtgacgcatgtcgaccttgtgaggggacacaaagaaaaattgaataaactgtatgatttaaaatat                                                                        <  1:302178‑M1/250‑174 (MQ=255)
                                                                                                                                                                                   ggCGGTCTTCGATCGGGCAGCGTTAATGCTTATCTGCCTGTTCTTTCTCATCCGTATCCGCCTGTTTCGCGccccgttcagtgatagttgacgagtgaaattatgaacgtgcaaaacaaatcgccgttaccagacaaataatgcatttactggagaattaatcccgctgtacggaagtgtgacgcaggtcgaccttgtgaggggacacaaagaaaaattgaataaactgtatgatttaaaagattatcg                                                                   <  1:9070‑M1/249‑179 (MQ=255)
                                                                                                                                                                                       gTCTTCGATCGGGCAGCGTTAATGCTTATCTGCCTGTTCTTTCTCATCCGTATCCGCCTGTTTCGCGccccgttcagtgatagttgacgagtgaaattatgaacgtgcaaaacaaatcgccgttaccagacaaataatgcatttactggagaattaatcccgctgtacggaagtgtgacgcaggtcgaccttgtgaggggacacaaagaaaaattgaataaactgtatgatttaaaagattaacgggagag                                                             <  1:379852‑M1/251‑185 (MQ=255)
                                                                                                                                                                                            cGATCGGGCAGCGTTAATGCTTATCTGCCTGTTCTTTCTCATCCGTATCCGCCTGTTTCGCGccccgttcagtgatagttgacgagtgaaattatgaacgtgcaaaacaaatcgccgttaccagacaaataatgcatttactggagaattaatcccgctgtacggaagtgtgacgcaggtcgaccttgtgaggggacacaaagaaaaattgaataaactgtatgatttaaaagattatcgggagagttacc                                                        >  1:97659‑M1/1‑62 (MQ=255)
                                                                                                                                                                                                     aGCGTTAATGCTTATCTCCCTGTTCTTTCTCATCCGTATCCGCCTGTTTCGCGccccgttcagtgatagttgacgagtgaaattatgaacgtgcaaaacaaatcgccgttaccagacaaataatgcatgtactggagaattaatcccgctgtacggaagtgtgacgcaggtcgaccttgtgaggggacacaaagaaaaatggaataaactgtattatttaaaagattatcgggagagttacctcccgatat                                               <  1:395820‑M1/251‑199 (MQ=255)
                                                                                                                                                                                                        gTTAATGCTTATCTGCCTTTTCTTTCTCATCCGTATCCGCCTGTTTCGCGccccgttcagtgatagttgacgagtgaaattatgaacgtgcaaaacaaatcgccgttaccagacaaataatgcatttactggagaattaagcccgctgtacggaagtgtgacgcaggtcgaccttgtgaggggacacaaagaaaaattgaataaactgtatgatttaaaagagtatcgggagaggtacctcccgatataaa                                            <  2:182085‑M1/251‑202 (MQ=255)
                                                                                                                                                                                                        gTTAATGCTTATCTGCCTGTTCTTTCTCATCCGTATCCGCCTGTTTCGCGccccgttcagtgatagttgacgagtgaaattatgaacgtgcaaaacaaatcgccgttaccagacaaataatgcatttactggagaattaatcccgctgtacggaagtgtgacgcaggtcgaccttgtgaggggacacaaagaaaaattgaataaactgtatgatttaaaagattatcgggagagttacctcccgatataa                                             >  1:517265‑M1/1‑50 (MQ=255)
                                                                                                                                                                                                        gTTAATGCTTATCTGCCTGTTCTTTCTCATCCGTATCCGCCTGTTTCGCGccccgttcagtgatagttgacgagtgaaattatgaacgtgcaaaacaaatcgccgttaccagacaaataatgcatttactggagaattaatcccgctgtacggaagtgtgacgcaggtcgacct                                                                                                                         <  2:235535‑M1/174‑125 (MQ=255)
                                                                                                                                                                                                        gTTAATGCTTATCTGCCTGTTCTTTCTCATCCGTATCCGCCTGTTTCGCGccccgttcagtgatagttgacgagtgaaattatgaacgtgcaaaacaaatcgccgttaccagacaaataatgcatttactggagaattaatcccgctgtacggaagtgtgacgcaggtcgacct                                                                                                                         >  1:235535‑M1/1‑50 (MQ=255)
                                                                                                                                                                                                                                   cATCCGTATCCGCCTGTTTCGCGccccgttcagtgatagttgacgagtgaaattatgaacgtgcaaaacaaatcgccgttaccagacaaataatgcatttactggagaattaatcccgctgtacggaagtgtgacgcaggtcgaccttgtgaggggacacaaagaaaaattgaataaactgtatgatttaaaagattatcgggagagttacctcccgatataaaaggaaggatttacagaatgtgaccta                  <  1:140420‑M1/250‑228 (MQ=255)
                                                                                                                                                                                                                                      ccGTATCCGCCTGTTTCGCGccccgttcagtgatagttgacgagtgaaattatgaacgtgcaaaacaaatcgccgttaccagacaaataatgcatttactggagaattaatcccgctgtacggaagtgtgacgcaggtcgaccttgtgaggggacacaaagaaaaattgaataaactgtatgatttaaaagattatcgggagagttacctcccgatataaaaggaaggatttacagaatgtgacctaaggt              <  1:468457‑M1/251‑232 (MQ=255)
                                                                                                                                                                                                                                      ccGTATCCGCCTGTTTCGCGccccgttcagtgatagttgacgagtgaaattatgaacgtgcaaaacaaatcgccgttaccagaaaaataatgcatttactggagaattaatcccgctgtacggaagtgtgacgcaggtcgaccttgtgaggggacacaaagaaaaattgaataaactgtatgatttaaaagattatcgggagagttacctcccgatataaaaggaaggatttacagaatgtgacctaaggt              >  1:485609‑M1/1‑20 (MQ=255)
                                                                                                                                                                                                                                                  gTTTCGCGccccgttcagtgatagttgacgagtgaaattatgaacgtgcaaaacaaatcgccgttaccagacaaataatgcatttactggagaattaatcccgctgtacggaagtgtgacgcaggtcgaccttgtgaggggacacaaagaaaaattgaataaactgtatgatttaaaagattatcgggcgagttacctcccgatataaaaggaaggatttacagaatgtgacctaaggtctggcgtaaatg  >  1:394759‑M1/1‑8 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                   
AACCAGAAGTGCATTATCAAGAGATGCATTGCTAATACACCATTTATCTATCAATCCCCTTTTAAAATTGTGGTTTTTCCTCAAAAACGGCCTGCTGGTCATCAGCCTTCAGGTTACCTATCATAGAGGTTTAATCCTTATTCAGAGTCACCCGTGCACGAAATATTCAACATGCTGCTGGCGGTCTTCGATCGGGCAGCGTTAATGCTTATCTGCCTGTTCTTTCTCATCCGTATCCGCCTGTTTCGCGAACTGTTGCACAAGTCGGCGCACTCCCCAAAAGAATTGCTCGCCGTTACCGCCATTTTCTCGCTGTTCGCCCTGTTCAGCACCTGGTCCGGCGTTCCCGTAGAAGGCTCGCTGGTGAACGTACGTATTATCGCGGTGATGTCCGGCGGGATTCTGTTTGGCCCGTGGGTAGGCATCATTACCGGCGTGATTGCGGGTATTCACCGGTATTTAATTGATATCGGCGGCGTGACGGCGATCCCCT  >  NC_000913/2498518‑2499010

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