Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,173,076 C→A 22.8% pseudogene (287/442 nt) gatC ← pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,173,0760CA22.8% 267.3 / 78.2 136pseudogene (287/442 nt)gatCpseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC
Reads supporting (aligned to +/- strand):  ref base C (52/53);  new base A (15/16);  total (67/69)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.18e-01

Only 100 of 138 total aligned reads displayed.

CATGAAAATACCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAG  >  NC_000913/2172978‑2173175
                                                                                                  |                                                                                                   
 aTGAAAATACCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTTCGGGGAGAAAACCTGAATCAGTAACCAGTTAATGGGAGAACCGc                                                                                                   >  2:1844780/1‑100 (MQ=255)
  tGAAAATACCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGcc                                                                                                  >  2:1792709/1‑100 (MQ=255)
  tGAAAATACCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGcc                                                                                                  <  1:1970041/100‑1 (MQ=255)
  tGAAAATACCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGcc                                                                                                  >  1:3729181/1‑100 (MQ=255)
  tGAAAAGACCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAAACAGTAACCAGGTAATGGGGGAACCGcc                                                                                                  >  2:2978677/1‑100 (MQ=255)
  ggAAAATACCGGTCAGATAAATTGCACGGATAATAATGGAACGGGGAAGATTTGGGGGGGAGAAAAGCTGAATCAGGAACCAGGTAATGGGAGAACCg                                                                                                    >  2:3508214/2‑98 (MQ=255)
   gAAAATACCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGccc                                                                                                 >  1:2573516/1‑100 (MQ=255)
    aaaaTACCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCt                                                                                                >  1:3353316/1‑100 (MQ=255)
       aTACCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGAt                                                                                             >  2:3493116/1‑100 (MQ=255)
       aTACCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGAt                                                                                             <  1:2826802/100‑1 (MQ=255)
        tACCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATc                                                                                            <  1:4093583/100‑1 (MQ=255)
         aCCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATcc                                                                                           <  2:3768917/100‑1 (MQ=255)
          ncGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCa                                                                                          >  2:482146/2‑100 (MQ=255)
             gTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTg                                                                                       >  2:1062042/1‑100 (MQ=255)
             gTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTg                                                                                       <  1:3117105/100‑1 (MQ=255)
               cAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGaa                                                                                     >  1:2477660/1‑100 (MQ=255)
               cAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCCGATCCAttca                                                                                      >  2:2121512/1‑97 (MQ=255)
                aGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGaa                                                                                     >  1:1132248/1‑99 (MQ=255)
                aGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAg                                                                                    <  2:3176535/100‑1 (MQ=255)
                 gATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGc                                                                                   >  2:4056196/1‑100 (MQ=255)
                 gATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGc                                                                                   >  1:3524051/1‑100 (MQ=255)
                   tAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGaca                                                                                 <  2:3066398/100‑1 (MQ=255)
                     aaTTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGccacc                                                                               <  2:1336171/100‑1 (MQ=255)
                     aaTTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGCAATGGGAGAACCGCCCTGATCCATTGAAGccacc                                                                               >  2:180280/1‑100 (MQ=255)
                       ttGCACCGATAATAATGAAACCGGGCATATCTTGCGGCCAGAAAACCTGCATCAGTCACCAGGTACTGGGACAACCGCCCTGATCCATCGAACCCACCAt                                                                             >  2:3665843/1‑100 (MQ=255)
                        tGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATa                                                                            >  1:125002/1‑100 (MQ=255)
                        tGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATa                                                                            >  1:92790/1‑100 (MQ=255)
                         gCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATAc                                                                           >  2:352254/1‑100 (MQ=255)
                          cACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATAcc                                                                          >  1:2315948/1‑100 (MQ=255)
                          cACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATAcc                                                                          <  1:3371275/100‑1 (MQ=255)
                            ccGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATAc                                                                           >  2:1438404/1‑97 (MQ=255)
                                 aataatGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCt                                                                   <  2:1814108/100‑1 (MQ=255)
                                  ataatGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCtt                                                                  <  1:1545507/100‑1 (MQ=255)
                                  ataatGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCtt                                                                  >  1:106119/1‑100 (MQ=255)
                                   taatGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCttt                                                                 <  1:465529/100‑1 (MQ=255)
                                   taatGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCt                                                                   >  1:1710226/1‑98 (MQ=255)
                                    aatGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCtttt                                                                <  1:554153/100‑1 (MQ=255)
                                     atGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTa                                                               >  1:3974630/1‑100 (MQ=255)
                                      tGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTaa                                                              >  2:3772949/1‑100 (MQ=255)
                                      tGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTaa                                                              <  2:781199/100‑1 (MQ=255)
                                      tGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTaa                                                              <  2:2613543/100‑1 (MQ=255)
                                       gAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAAc                                                             <  2:863392/100‑1 (MQ=255)
                                           ccGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGcc                                                          >  2:946558/1‑99 (MQ=255)
                                            cGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGcccc                                                        <  2:3340889/100‑1 (MQ=255)
                                              ggAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAg                                                      <  1:508833/100‑1 (MQ=255)
                                              ggAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAg                                                      >  2:248694/1‑100 (MQ=255)
                                               gAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGc                                                     <  1:2018572/100‑1 (MQ=255)
                                                 atatTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAAGGGGAGAACAGCCCTGATCCATTTAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCAt                                                   <  2:2741606/100‑1 (MQ=255)
                                                  tatTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCAtt                                                  >  1:3798285/1‑100 (MQ=255)
                                                    ttttGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAg                                                <  1:1780453/100‑1 (MQ=255)
                                                    ttttGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAg                                                >  1:1141363/1‑100 (MQ=255)
                                                      ttGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAgcc                                              <  1:790784/100‑1 (MQ=255)
                                                      ttGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAgcc                                              <  2:3822300/100‑1 (MQ=255)
                                                      ttGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAgcc                                              <  2:3692148/100‑1 (MQ=255)
                                                         cGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAgccgcc                                           <  1:921606/100‑1 (MQ=255)
                                                         cGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAgccgcc                                           >  1:3222979/1‑100 (MQ=255)
                                                         cGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAgccgcc                                           >  1:3830273/1‑100 (MQ=255)
                                                         cGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAgccgcc                                           <  1:2558653/100‑1 (MQ=255)
                                                         cGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTCAAGCCACCATACCCCCGGCCTTTACCGCCCCAGCATTAgccgcc                                           >  2:3429670/1‑100 (MQ=255)
                                                              agaAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTg                                      <  2:1217371/100‑1 (MQ=255)
                                                               gacaaCCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCAGGCTCTTAACGCCCCAGCATTAGCCGCCAgctag                                     >  2:629316/4‑96 (MQ=255)
                                                                aaaaCCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTggg                                    <  1:1376169/100‑1 (MQ=255)
                                                                  aaCCTGAATCAGTAACCAGGTAATGGGAGAACCACCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGgtg                                  <  2:387747/100‑1 (MQ=255)
                                                                  aaCCTGAATCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCATTTGGgtg                                  <  2:3798261/100‑1 (MQ=255)
                                                                    ccTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGtgtg                                >  1:3600772/1‑100 (MQ=255)
                                                                    ccTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGtgtg                                <  2:1048390/100‑1 (MQ=255)
                                                                    ccTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGtgtg                                >  1:3817234/1‑100 (MQ=255)
                                                                      tGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACACCCCAGCATTAGCCGCCAGTTGGGTGTGCa                              >  1:602466/1‑100 (MQ=255)
                                                                        aaTCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCaaa                            <  2:592148/100‑1 (MQ=255)
                                                                         aTCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAAc                           <  1:1653782/100‑1 (MQ=255)
                                                                          tCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAAcc                          <  1:3262976/100‑1 (MQ=255)
                                                                          tCAGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAAcc                          <  1:3420555/100‑1 (MQ=255)
                                                                            aGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCaa                        >  2:3729976/1‑100 (MQ=255)
                                                                            aGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCaa                        >  1:3481585/1‑100 (MQ=255)
                                                                            aGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCaa                        >  1:2970206/1‑100 (MQ=255)
                                                                            aGTAACCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATGCCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCaa                        >  2:1978080/1‑100 (MQ=255)
                                                                              tAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATa                      <  2:687484/100‑1 (MQ=255)
                                                                               aaCCAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAg                     <  2:3935129/100‑1 (MQ=255)
                                                                                aCCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGt                    >  2:1648941/1‑100 (MQ=255)
                                                                                  cAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGttt                  >  2:4006729/1‑100 (MQ=255)
                                                                                  cAGGTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGttt                  <  1:2004007/100‑1 (MQ=255)
                                                                                   aGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTg                 <  1:2412189/100‑1 (MQ=255)
                                                                                   aGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTg                 <  1:2855179/100‑1 (MQ=255)
                                                                                    ggTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGc                <  2:949563/100‑1 (MQ=255)
                                                                                     gTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCg               <  1:774117/100‑1 (MQ=255)
                                                                                     gTAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCg               <  2:3729157/100‑1 (MQ=255)
                                                                                      tAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTCACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGCTTGCg               >  2:207790/1‑99 (MQ=255)
                                                                                      tAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACTATTGGTTTGcac              <  1:2032598/100‑3 (MQ=255)
                                                                                      tAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACTATTGGTTTGcac              <  1:1438414/100‑3 (MQ=255)
                                                                                      tAATGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGt              >  1:4013724/1‑100 (MQ=255)
                                                                                       aaTGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTACCCGCCAGTTGGGTGTGCAAACCCATAGTTTGCGTc             >  2:3909116/1‑100 (MQ=255)
                                                                                        aTGGGAGACCCGCCCTGATCCATTGAGGCCTCCGAACCCCCGGCTTTGTACGGCTCAGCCTTTGCCGCCAGCGGGGTGTGCAGACAATTGGTTTGcacac            <  2:1189620/100‑5 (MQ=255)
                                                                                         tGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTAGGTGTGCAAACCAATAGTTTGCGTcgc           <  1:352257/100‑1 (MQ=255)
                                                                                         tGGGAGAACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCCGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTcgc           >  2:2351441/1‑100 (MQ=255)
                                                                                         tGGGAGAACAGCCCTGATCCATTGAAGCCACCAGACCCCCGGCTTTTAACGCCCCAGCACTAGCCGCCACATGGGTGTGCAAACCACTAGTTTACGTc             >  2:2482333/1‑98 (MQ=255)
                                                                                           ggAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGa         >  1:682643/1‑100 (MQ=255)
                                                                                            gagaACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGAt        <  2:322776/100‑1 (MQ=255)
                                                                                              gaACAGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATcc      >  1:995098/1‑100 (MQ=255)
                                                                                               aaCCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCca     <  1:3715611/100‑1 (MQ=255)
                                                                                                  cGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAg  <  1:3857770/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
CATGAAAATACCGGTCAGATAAATTGCACCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAG  >  NC_000913/2172978‑2173175

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