Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,078,186 C→T 15.3% G225D (GGC→GAC)  cmtA ← putative mannitol‑specific PTS IIB and IIC components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,078,1860CT15.3% 415.8 / 51.5 167G225D (GGC→GAC) cmtAputative mannitol‑specific PTS IIB and IIC components
Reads supporting (aligned to +/- strand):  ref base C (70/70);  new base T (17/9);  total (87/80)
Fisher's exact test for biased strand distribution p-value = 2.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

Only 100 of 174 total aligned reads displayed.

TCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCAT  >  NC_000913/3078087‑3078285
                                                                                                   |                                                                                                   
tCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGc                                                                                                     <  2:631459/100‑1 (MQ=255)
  gTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGccc                                                                                                   >  2:1598369/1‑100 (MQ=255)
   tGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGccc                                                                                                   >  1:1912863/1‑99 (MQ=255)
   tGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCg                                                                                                  <  1:1518120/100‑1 (MQ=255)
   tGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCg                                                                                                  <  1:1570959/100‑1 (MQ=255)
    ggATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCgg                                                                                                 <  1:3816897/100‑1 (MQ=255)
    ggATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCg                                                                                                  >  1:2499928/1‑99 (MQ=255)
     gATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCAGCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGa                                                                                                <  1:2115443/100‑1 (MQ=255)
       tCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGTCCGGGcc                                                                                              <  2:3448508/100‑1 (MQ=255)
         ccACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGAccag                                                                                            >  1:664103/1‑98 (MQ=255)
         ccACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTg                                                                                            >  2:3089903/1‑100 (MQ=255)
         ccACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTg                                                                                            <  1:1225780/100‑1 (MQ=255)
         ccACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTg                                                                                            >  1:533763/1‑100 (MQ=255)
         ccACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTg                                                                                            >  2:3695414/1‑100 (MQ=255)
         ccACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGCAAACGCCAGCAGCAGCCCCAGGCCCGGACCt                                                                                             >  1:301457/1‑99 (MQ=255)
           aCCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTggg                                                                                          <  2:3334055/100‑1 (MQ=255)
            ccGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCACCAGCCCCAGGCCCGGAc                                                                                               >  1:2966899/1‑94 (MQ=255)
             cGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGTCCGGGCCTggggc                                                                                        <  2:128258/100‑3 (MQ=255)
             cGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGtt                                                                                        <  2:2697497/100‑1 (MQ=255)
               aGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGTCCGGGCCTggggctg                                                                                      <  2:683538/100‑5 (MQ=255)
                ggAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGTCCGGGCCTggggctgc                                                                                     <  2:780687/100‑6 (MQ=255)
                ggAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTaga                                                                                     <  1:3208440/100‑1 (MQ=255)
                 gAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTagag                                                                                    <  1:260072/100‑1 (MQ=255)
                  atctGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAgg                                                                                   <  2:24440/97‑1 (MQ=255)
                  aaGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAgg                                                                                   <  2:1714022/100‑1 (MQ=255)
                   aGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGc                                                                                  <  1:3534184/100‑1 (MQ=255)
                   aGTGAATAATCAGCGCCCCGGGCGCAAAACGTTTACTCATCCCTTTACCAAGCAAGGTAAACGCCAGCAGCAGCACCAGGCCCGGACCTGGGTTAGAgg                                                                                   <  1:1446632/99‑1 (MQ=255)
                    gTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAACCGCCAGCAGCAGCCGCAGGCCCGGACCTGGGTTAGAGGcc                                                                                 <  1:1393877/100‑1 (MQ=255)
                    gTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGcc                                                                                 <  2:2065647/100‑1 (MQ=255)
                    gTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGcc                                                                                 >  2:1665121/1‑100 (MQ=255)
                    gTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGcc                                                                                 >  1:1925046/1‑100 (MQ=255)
                      gAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGccac                                                                               >  1:1249728/1‑100 (MQ=255)
                      gAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGccac                                                                               <  2:2238614/100‑1 (MQ=255)
                      gAATAATCATCGCACCGGGCGCAGAACGTTTACTCATCACTTTACCAAACAAGGTAAACGTCAGCAGCAGCCCCAGCCCCGGACCTGGGTTAGAgg                                                                                   >  1:24439/1‑96 (MQ=255)
                       natcaccatcGCCCCGGGAGCAGAACGTTTACTCAAACCTTAACCAAACACGGCAAACGCCAGCAGCAGCCCCAGGCCCGCACCTGGGTTAGAGGccacc                                                                              >  1:1660612/7‑100 (MQ=255)
                        ataatCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGccacca                                                                             <  2:3656626/100‑1 (MQ=255)
                        ataatCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGccacca                                                                             >  1:1001169/1‑100 (MQ=255)
                          aatCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATa                                                                           <  2:3417776/100‑1 (MQ=255)
                          aatCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATa                                                                           <  2:526162/100‑1 (MQ=255)
                           atcatcGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATaa                                                                          <  2:3064990/100‑1 (MQ=255)
                             catcGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCAGACCTAGGTTAGAGGCCACCATa                                                                           >  1:1168488/1‑97 (MQ=255)
                              atcGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAaaga                                                                       <  2:1816968/100‑1 (MQ=255)
                              atcGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAaaga                                                                       >  2:1595439/1‑100 (MQ=255)
                                      gggCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGAt                                                               <  2:2902477/100‑1 (MQ=255)
                                       ggCGCAGAACGTTTACTCGACCCTTTGCCAGGCACGGTAAGCGCCATCAGCAGGCCGTGGTCTGGACGTGGGTGAGAGGCCACGATTAAGAAGATGGAtt                                                              <  1:2618967/100‑1 (MQ=255)
                                         cgcAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTa                                                            <  2:24519/100‑1 (MQ=255)
                                            aGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCg                                                         >  2:3381713/1‑100 (MQ=255)
                                            aGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCg                                                         >  2:2644267/1‑100 (MQ=255)
                                                cGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTaa                                                     >  2:792738/1‑100 (MQ=255)
                                                 gTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAAc                                                    >  1:812714/1‑100 (MQ=255)
                                                 gTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAAc                                                    >  2:3791310/1‑100 (MQ=255)
                                                    tACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGa                                                 <  2:2171614/100‑1 (MQ=255)
                                                       tCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGc                                              >  2:1423633/1‑100 (MQ=255)
                                                        cATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGcc                                             >  2:2774958/1‑100 (MQ=255)
                                                        cATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGcc                                             >  1:3452197/1‑100 (MQ=255)
                                                        cATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGcc                                             <  1:1630227/100‑1 (MQ=255)
                                                         aTCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCt                                            >  2:3596186/1‑100 (MQ=255)
                                                            ccTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGtt                                         <  1:522964/100‑1 (MQ=255)
                                                            ccTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGCCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTgat                                         <  2:2578232/100‑3 (MQ=255)
                                                            ccTTTACCAAACAAGGAAAAAGCCAGCAGCAGCCCCAGGCCCGGACCTGGGATAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGcat                                            >  1:2859052/1‑95 (MQ=255)
                                                             cTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGACGATGGATTTACCGTTAACCGAAGCCTGTTg                                        <  1:715815/100‑1 (MQ=255)
                                                             cTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTg                                        <  2:812722/100‑1 (MQ=255)
                                                              tttACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTATAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTg                                        <  1:1003852/99‑1 (MQ=255)
                                                              tttACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGc                                       <  2:3542579/100‑1 (MQ=255)
                                                                tACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCAt                                     >  2:214093/1‑100 (MQ=255)
                                                                  ccAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGccagatcggaag                                   >  2:1922271/1‑90 (MQ=255)
                                                                   cAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTcc                                  >  1:3136629/1‑100 (MQ=255)
                                                                    aaaCAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTccc                                 >  2:1692829/1‑100 (MQ=255)
                                                                     aaCAAGGGAAACGCCAGCAGCAGCCCCAGGACCGGACCTGGGTGAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTcca                                 <  1:3355485/99‑2 (MQ=255)
                                                                          ggTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGa                           >  1:3367995/1‑100 (MQ=255)
                                                                           gTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGan                          <  2:3142548/100‑2 (MQ=255)
                                                                            tAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGata                         <  2:221732/100‑1 (MQ=255)
                                                                               aCGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGataata                      <  1:3647518/100‑1 (MQ=255)
                                                                                cGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGataat                       <  1:570916/98‑1 (MQ=255)
                                                                                 gCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGa                    <  1:851180/100‑1 (MQ=255)
                                                                                  ccAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGAc                   >  2:558092/1‑100 (MQ=255)
                                                                                     gcagcagcCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGc                 >  1:3054037/1‑99 (MQ=255)
                                                                                      cctctaacccAGGTCCGGGCCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGccc               >  2:1734340/8‑100 (MQ=255)
                                                                                      cctctaacccAGGTCCGGGCCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGc                 >  1:1855749/8‑98 (MQ=39)
                                                                                      cctctaacccAGGTCCGGGCCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCAAGCGGATAATAGACGc                 >  1:1517950/8‑98 (MQ=38)
                                                                                      cagcagcCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGccc               <  2:1522538/100‑1 (MQ=255)
                                                                                       ctctaacccAGGTCCGGGCCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCt              >  2:2092499/7‑100 (MQ=255)
                                                                                        tctaacccAGGTCCGGGCCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTg             >  2:3019881/6‑100 (MQ=255)
                                                                                         ctaacccAGGTCCGGGCCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTga            >  2:2948922/5‑100 (MQ=255)
                                                                                         ctaacccAGGTCCGGGCCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTga            >  2:2872833/5‑100 (MQ=255)
                                                                                         ctaacccAGGTCCGGGCCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCACTga            >  2:591162/5‑100 (MQ=255)
                                                                                          taacccAGGTCCGGGCCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTga            >  1:3694814/4‑99 (MQ=255)
                                                                                          gcccccAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTgat           <  1:1938500/98‑1 (MQ=255)
                                                                                          agcCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTgat           >  2:1919630/1‑100 (MQ=255)
                                                                                           gcCCCAGGCCCGGGCCTGGGTTAGAGGCCACCATAAAGTAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTgatc          >  1:2352559/1‑100 (MQ=255)
                                                                                           aacccAGGTCCGGGCCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTgatc          >  1:295562/3‑100 (MQ=255)
                                                                                           aacccAGGTCCGGGCCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTgatc          >  2:641058/3‑100 (MQ=255)
                                                                                           aacccAGGTCCGGGCCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTgatc          >  2:886367/3‑100 (MQ=255)
                                                                                            ccccAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTgatcg         >  2:407925/1‑100 (MQ=255)
                                                                                             cccAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTgatcga        >  2:671212/1‑100 (MQ=255)
                                                                                             cccAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTgatcga        >  1:1535922/1‑100 (MQ=255)
                                                                                                 ggCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGc    <  1:760409/100‑1 (MQ=255)
                                                                                                  gCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCa   <  2:136298/100‑1 (MQ=255)
                                                                                                  gCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCa   <  2:2665274/100‑1 (MQ=255)
                                                                                                   cccGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCat  <  2:1150548/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
TCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCAT  >  NC_000913/3078087‑3078285

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: