Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,619,018 A→G T264A (ACG→GCG)  deoA → thymidine phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,619,0180AG100.0% 314.8 / NA 83T264A (ACG→GCG) deoAthymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (35/48);  total (35/48)

GAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAA  >  NC_000913/4618871‑4619161
                                                                                                                                                   |                                                                                                                                               
gAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCgg                                                                                                                                              >  1:1105851/1‑151 (MQ=255)
 aaCTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCggg                                                                                                                                             <  2:113960/151‑1 (MQ=255)
    tctctGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTg                                                                                                                                           >  2:597965/1‑150 (MQ=255)
     ctctGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGa                                                                                                                                          >  1:159653/1‑150 (MQ=255)
       ctGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGa                                                                                                                                          <  2:14867/148‑1 (MQ=255)
        tGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAatat                                                                                                                                      >  2:485988/1‑151 (MQ=255)
         gAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAat                                                                                                                                        <  2:593219/148‑1 (MQ=255)
          aaGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGGGAata                                                                                                                                       >  1:634303/1‑148 (MQ=255)
            gCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTa                                                                                                                                  <  2:345620/151‑1 (MQ=255)
            gCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTa                                                                                                                                  <  1:740214/151‑1 (MQ=255)
                 tGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAAccc                                                                                                                              <  2:85437/150‑1 (MQ=255)
                  gCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCcgc                                                                                                                            <  1:52264/151‑1 (MQ=255)
                    cGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCcgcg                                                                                                                           >  1:20439/1‑150 (MQ=255)
                    cGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCcgcg                                                                                                                           >  2:36138/1‑150 (MQ=255)
                    cGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGt                                                                                                                          >  2:80821/1‑151 (MQ=255)
                     gAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCcgcg                                                                                                                           >  2:1178596/1‑149 (MQ=255)
                     gAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTc                                                                                                                         <  1:764157/151‑1 (MQ=255)
                     gAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGAGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTc                                                                                                                         <  1:826949/151‑1 (MQ=255)
                       aGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGGGAATATCGTAACCCGCGTCt                                                                                                                        >  1:1046221/1‑150 (MQ=255)
                         cGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGtt                                                                                                                     >  2:613794/1‑151 (MQ=255)
                           aTTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTg                                                                                                                   <  2:15864/151‑1 (MQ=255)
                              gttgGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGAt                                                                                                                 <  2:1030698/150‑1 (MQ=255)
                                tgGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGt                                                                                                               <  2:110229/150‑1 (MQ=255)
                                  gCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAAACCGCGGCTGTTTTATGTc                                                                                                              >  1:929313/1‑149 (MQ=255)
                                    gTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACg                                                                                                           <  1:980504/150‑1 (MQ=255)
                                        cTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATgg                                                                                                       >  2:607806/1‑150 (MQ=255)
                                         tAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGgcg                                                                                                     >  1:215998/1‑151 (MQ=255)
                                           aCGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCt                                                                                                   <  1:460601/151‑1 (MQ=255)
                                           aCGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAAAATAGTAACCCGGGTCTGTTTGATGTCGCGCTgg                                                                                                       >  1:848782/1‑147 (MQ=255)
                                                gcTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGGGAATATTGTAACCCGCGTCTGTTTGATGTCACGATGgcgc                                                                                                    >  1:449905/1‑145 (MQ=255)
                                                gcTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGGGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCtgt                                                                                                 >  2:664682/1‑148 (MQ=255)
                                                 cTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGGGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCg                                                                                              >  1:887111/1‑150 (MQ=255)
                                                    gCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAAACCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGt                                                                                             >  1:652423/1‑148 (MQ=255)
                                                      gTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGgaga                                                                                        <  1:848875/151‑1 (MQ=255)
                                                        gcgcACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGAt                                                                                       <  2:1009337/150‑1 (MQ=255)
                                                             ccaccGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGAt                                                                                 <  1:576892/151‑1 (MQ=255)
                                                              caccGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGAtc                                                                                >  2:1171181/1‑151 (MQ=255)
                                                                  gcgcTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTcc                                                                             >  2:739001/1‑150 (MQ=255)
                                                                   cgcTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCgg                                                                           <  2:240530/151‑1 (MQ=255)
                                                                   cgcTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCgg                                                                           <  1:570287/151‑1 (MQ=255)
                                                                    gctgctCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGc                                                                          <  2:656093/151‑1 (MQ=255)
                                                                     ctgctCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGc                                                                          >  1:557429/1‑150 (MQ=255)
                                                                         tCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACt                                                                     >  2:236173/1‑151 (MQ=255)
                                                                          cACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTg                                                                    <  2:185291/151‑1 (MQ=255)
                                                                          cACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTg                                                                    <  1:1064848/151‑1 (MQ=255)
                                                                           aCCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTgg                                                                   <  2:159584/151‑1 (MQ=255)
                                                                            ccGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGc                                                                  <  2:544813/151‑1 (MQ=255)
                                                                             cGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCg                                                                 <  1:103596/151‑1 (MQ=255)
                                                                                cATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGaaa                                                              <  1:1178334/151‑1 (MQ=255)
                                                                                   gAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGAt                                                           <  2:766946/151‑1 (MQ=255)
                                                                                   gAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGAt                                                           <  1:1159258/151‑1 (MQ=255)
                                                                                    aaTCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATg                                                          <  1:645788/151‑1 (MQ=255)
                                                                                         ggTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGcc                                                     <  2:18423/151‑1 (MQ=255)
                                                                                           tACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGGGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGGGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGcc                                                     >  1:1126908/1‑149 (MQ=255)
                                                                                            aCTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCg                                                    >  1:1185712/1‑149 (MQ=255)
                                                                                              tGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAgc                                                >  2:525255/1‑151 (MQ=255)
                                                                                               ggCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAgcg                                               <  1:923271/151‑1 (MQ=255)
                                                                                               ggCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAg                                                 <  2:1011996/149‑1 (MQ=255)
                                                                                               agCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGAAGGCGCTGTGCGTGGGGATGCTGATCTCCGGCAAAATGGCGAAAGATGACGCCGAAgcg                                               >  1:1062879/2‑151 (MQ=255)
                                                                                                gCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAgcgc                                              <  2:1008248/151‑1 (MQ=255)
                                                                                                gCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAgcgc                                              <  1:693919/151‑1 (MQ=255)
                                                                                                    ccAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAgcgcgcg                                           >  1:20059/1‑150 (MQ=255)
                                                                                                    ccAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAgcgcgc                                            <  1:1022314/149‑1 (MQ=255)
                                                                                                     cAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGcgcgcgcg                                         <  2:550518/151‑1 (MQ=255)
                                                                                                       gTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCTAAGCGCGCGCGaa                                       <  1:571004/151‑1 (MQ=255)
                                                                                                         gCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAAt                                     <  1:1110491/151‑1 (MQ=255)
                                                                                                           aGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGGGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGGGCTGTGCGTGGAGATGCTGATATCCGGCAAACTGGCGCAAGAAGGCGCCGAAGCGCGCGCGCAAt                                     >  1:647726/1‑149 (MQ=255)
                                                                                                             gTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCa                                 <  2:466659/151‑1 (MQ=255)
                                                                                                               aaCGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAg                                <  1:485870/150‑1 (MQ=255)
                                                                                                                        gAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCTGAAGCGCGCGCGAAATTGCAGGCGGTGCTgg                      <  2:457311/151‑1 (MQ=255)
                                                                                                                        gAAGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTgg                      <  1:1183088/151‑1 (MQ=255)
                                                                                                                         aaGTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGa                     <  2:159687/151‑1 (MQ=255)
                                                                                                                           gTTCGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTgg                      >  1:585114/1‑148 (MQ=255)
                                                                                                                              cGTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACg                <  1:36130/151‑1 (MQ=255)
                                                                                                                               gTGAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGGGCTGGGCaa                  >  1:719531/1‑148 (MQ=255)
                                                                                                                                 gAAGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTa             <  2:329243/151‑1 (MQ=255)
                                                                                                                                  aaGCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTaa            <  1:597822/151‑1 (MQ=255)
                                                                                                                                    gCGGTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAAAGGGaaa           >  2:928584/1‑150 (MQ=255)
                                                                                                                                       gTGCAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGGAAAACgg       >  1:1170066/1‑151 (MQ=255)
                                                                                                                                        tGCAGTTCCTGGCGGGGGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGTTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAACGCGCGCGAAATTGCCGGCGGCGGGGGAGAACAGGAAAACgg       >  1:573658/1‑150 (MQ=255)
                                                                                                                                          cAGTTCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAg    <  2:1106102/151‑1 (MQ=255)
                                                                                                                                            gTTCCTGGGGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAATCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGaa  <  2:982861/151‑1 (MQ=255)
                                                                                                                                              tCCTGGCGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGGCAACGGTAAAGCg        >  1:996497/1‑143 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
GAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAA  >  NC_000913/4618871‑4619161

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