Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 282,361 C→T G28G (GGC→GGT)  yagE → 2‑keto‑3‑deoxy gluconate (KDG) aldolase; CP4‑6 prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913282,3610CT100.0% 532.0 / NA 131G28G (GGC→GGT) yagE2‑keto‑3‑deoxy gluconate (KDG) aldolase; CP4‑6 prophage
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (67/64);  total (67/64)

Only 100 of 131 total aligned reads displayed.

ATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCG  >  NC_000913/282213‑282508
                                                                                                                                                    |                                                                                                                                                   
     cGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTAccgcc                                                                                                                                               <  1:1970018/150‑1 (MQ=255)
     cGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTAccgcc                                                                                                                                               <  2:2077158/150‑1 (MQ=255)
      gTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACcgccgc                                                                                                                                             >  2:1720958/1‑151 (MQ=255)
         aaataaaCGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCt                                                                                                                                          <  2:1606332/151‑1 (MQ=255)
         aaataaaCGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCt                                                                                                                                          >  1:1574578/1‑151 (MQ=255)
         aaataaaCGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCt                                                                                                                                          <  1:365720/151‑1 (MQ=255)
         aaataaaCGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCt                                                                                                                                          >  1:428705/1‑151 (MQ=255)
          aataaaCGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTg                                                                                                                                         >  1:1059833/1‑151 (MQ=255)
           ataaaCGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTg                                                                                                                                         >  1:2155902/1‑150 (MQ=255)
            taaaCGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGAt                                                                                                                                       >  1:446492/1‑151 (MQ=255)
              aacgaacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATcg                                                                                                                                     >  1:224813/1‑151 (MQ=255)
              aacgaacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATcg                                                                                                                                     <  1:1340779/151‑1 (MQ=255)
                cgaacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATcgac                                                                                                                                   >  2:2078748/1‑151 (MQ=255)
                cgaacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATcgac                                                                                                                                   >  2:1734395/1‑151 (MQ=255)
                  aacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATcgacga                                                                                                                                 <  1:816421/151‑1 (MQ=255)
                  aacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATcgacg                                                                                                                                  >  2:108723/1‑150 (MQ=255)
                  aacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATcgacg                                                                                                                                  >  2:1504588/1‑150 (MQ=255)
                   acgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGAt                                                                                                                                >  2:832784/1‑151 (MQ=255)
                     gCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCt                                                                                                                              >  1:580343/1‑151 (MQ=255)
                       tGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTg                                                                                                                             >  1:931896/1‑150 (MQ=255)
                        gTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTg                                                                                                                             <  1:784955/149‑1 (MQ=255)
                         ttCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATc                                                                                                                          >  2:622093/1‑151 (MQ=255)
                         ttCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATc                                                                                                                          >  1:1740717/1‑151 (MQ=255)
                           cTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCa                                                                                                                         >  1:627025/1‑150 (MQ=255)
                             ataATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCCCCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCaa                                                                                                                        <  1:71411/149‑1 (MQ=255)
                             ataATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCaa                                                                                                                        <  2:1915908/149‑1 (MQ=255)
                             ataATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCa                                                                                                                         <  2:495389/148‑1 (MQ=255)
                               aaTGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCCCCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGc                                                                                                                     <  2:1264226/150‑1 (MQ=255)
                                aTGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGc                                                                                                                     >  2:134886/1‑149 (MQ=255)
                                aTGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGc                                                                                                                     >  1:680873/1‑149 (MQ=255)
                                aTGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGc                                                                                                                     >  1:1003744/1‑149 (MQ=255)
                                aTGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGc                                                                                                                     >  1:2069800/1‑149 (MQ=255)
                                 tGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAgg                                                                                                                  >  1:1678055/1‑151 (MQ=255)
                                   tAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCg                                                                                                                >  1:2158558/1‑151 (MQ=255)
                                         aaaaTGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACg                                                                                                          >  1:1877597/1‑151 (MQ=255)
                                          aaaTGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACgg                                                                                                         >  2:1913469/1‑151 (MQ=255)
                                           aaTGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGc                                                                                                        <  2:1699538/151‑1 (MQ=255)
                                                  cAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCGCCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGtt                                                                                                   <  1:1734426/149‑1 (MQ=255)
                                                  cAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGttct                                                                                                 >  2:201113/1‑151 (MQ=255)
                                                   aGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGttctt                                                                                                <  2:2122326/151‑1 (MQ=255)
                                                    gCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGttcttc                                                                                               >  1:268783/1‑151 (MQ=255)
                                                     cAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTcc                                                                                              >  1:1462269/1‑151 (MQ=255)
                                                         gagaTCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTggg                                                                                          >  2:1635367/1‑151 (MQ=255)
                                                            aTCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGc                                                                                         >  1:487555/1‑149 (MQ=255)
                                                            aTCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGc                                                                                         >  1:1154924/1‑149 (MQ=255)
                                                               tcATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCgg                                                                                    >  1:1464380/1‑151 (MQ=255)
                                                                cATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGc                                                                                      <  2:340565/148‑1 (MQ=255)
                                                                cATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGt                                                                                   >  2:817298/1‑151 (MQ=255)
                                                                 aTGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTg                                                                                  <  1:559576/151‑1 (MQ=255)
                                                                         gTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTtc                                                                          <  1:1320388/151‑1 (MQ=255)
                                                                         gTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTtc                                                                          >  2:1153279/1‑151 (MQ=255)
                                                                           ccGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTtct                                                                         <  1:388685/150‑1 (MQ=255)
                                                                            cgcgTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTcc                                                                       >  1:1448439/1‑151 (MQ=255)
                                                                              cgTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCa                                                                     <  1:2124082/151‑1 (MQ=255)
                                                                                ttgttCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGc                                                                   >  1:1967387/1‑151 (MQ=255)
                                                                                   ttCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTc                                                                 <  2:480647/150‑1 (MQ=255)
                                                                                    tCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCgg                                                               >  1:190814/1‑151 (MQ=255)
                                                                                    tCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCg                                                                <  1:1408064/150‑1 (MQ=255)
                                                                                      aCGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGgcg                                                             >  1:1145571/1‑151 (MQ=255)
                                                                                          gAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCg                                                          >  2:178688/1‑150 (MQ=255)
                                                                                           aaTCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGaa                                                        <  2:1668810/151‑1 (MQ=255)
                                                                                           aaTCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGaa                                                        <  1:100509/151‑1 (MQ=255)
                                                                                           aaTCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGa                                                         >  2:695666/1‑150 (MQ=255)
                                                                                           aaTCATTCCCCCTGTCTCCACCATTTTTACCACCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGaa                                                        <  2:1243194/151‑1 (MQ=255)
                                                                                            aTCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAag                                                       <  2:590072/151‑1 (MQ=255)
                                                                                            aTCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAag                                                       <  1:1618261/151‑1 (MQ=255)
                                                                                            aTCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAag                                                       <  2:1936741/151‑1 (MQ=255)
                                                                                              cATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAagag                                                     <  1:1377979/151‑1 (MQ=255)
                                                                                              cATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAag                                                       >  1:1147666/1‑149 (MQ=255)
                                                                                               aTTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGc                                                    >  2:1173401/1‑151 (MQ=255)
                                                                                                ttCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGc                                                    <  2:543020/150‑1 (MQ=255)
                                                                                                ttCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCg                                                   >  2:659878/1‑151 (MQ=255)
                                                                                                    cccTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTaaa                                               <  1:1785287/151‑1 (MQ=255)
                                                                                                    cccTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTaaa                                               <  2:1462616/151‑1 (MQ=255)
                                                                                                    cccTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTaaa                                               <  2:1177924/151‑1 (MQ=255)
                                                                                                    cccTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTaa                                                >  2:2101855/1‑150 (MQ=255)
                                                                                                     ccTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAg                                              <  1:378409/151‑1 (MQ=255)
                                                                                                        gTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGc                                             <  1:89030/149‑1 (MQ=255)
                                                                                                          ctcCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCa                                           <  2:435607/149‑1 (MQ=255)
                                                                                                          ctcCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCAtt                                         >  2:359030/1‑151 (MQ=255)
                                                                                                             caccaTTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGcc                                      <  2:602304/151‑1 (MQ=255)
                                                                                                                   ttttACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCtt                                 <  2:1592861/150‑1 (MQ=255)
                                                                                                                    tttACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCttt                                <  1:1297701/150‑1 (MQ=255)
                                                                                                                    tttACCGCCGACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCttt                                <  1:1061411/150‑1 (MQ=255)
                                                                                                                         cgccgACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATc                          >  1:878704/1‑151 (MQ=255)
                                                                                                                           ccgACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCg                         >  1:2005750/1‑150 (MQ=255)
                                                                                                                           ccgACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGa                        <  1:2170223/151‑1 (MQ=255)
                                                                                                                            cgACGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGAt                       <  2:1440570/151‑1 (MQ=255)
                                                                                                                              aCGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCa                     <  1:1227835/151‑1 (MQ=255)
                                                                                                                               cGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCAt                    <  1:605431/151‑1 (MQ=255)
                                                                                                                               cGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCAt                    <  2:1406594/151‑1 (MQ=255)
                                                                                                                               cGGCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCAt                    >  2:1722574/1‑151 (MQ=255)
                                                                                                                                ggCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATg                   <  1:1866989/151‑1 (MQ=255)
                                                                                                                                 gCCAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGt                  <  2:507678/151‑1 (MQ=255)
                                                                                                                                  ccAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTc                 >  2:922999/1‑151 (MQ=255)
                                                                                                                                   cAGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCg                >  1:812432/1‑151 (MQ=255)
                                                                                                                                    aGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGa               <  1:680950/151‑1 (MQ=255)
                                                                                                                                    aGCTCGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGa               <  1:1517514/151‑1 (MQ=255)
                                                                                                                                        cGATAAGCCGGGTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGAtcg            <  2:2077439/150‑1 (MQ=255)
                                                                                                                                                 gggTACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCg  >  2:668111/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACCGCCGCGCTGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCG  >  NC_000913/282213‑282508

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