Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,078,191 A→G 16.1% G223G (GGT→GGC 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,1910AG16.1% 383.2 / 57.9 166G223G (GGT→GGCcmtAputative mannitol‑specific PTS IIB and IIC components
Reads supporting (aligned to +/- strand):  ref base A (69/70);  new base G (19/8);  total (88/78)
Fisher's exact test for biased strand distribution p-value = 5.86e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 169 total aligned reads displayed.

GATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTG  >  NC_000913/3078092‑3078290
                                                                                                   |                                                                                                   
 aTCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGAc                                                                                                    >  2:1238833/1‑100 (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)
      aCCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTggg                                                                                               >  2:851187/1‑100 (MQ=255)
       ccGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTggg                                                                                               <  1:3544540/99‑1 (MQ=255)
        ctaGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGtt                                                                                             <  1:2664913/98‑1 (MQ=255)
        cGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGTCCGGGCCTggggc                                                                                             <  1:1171229/100‑3 (MQ=255)
        cGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGTCCGGGCCTggggc                                                                                             <  2:128258/100‑3 (MQ=255)
        cGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGtt                                                                                             <  2:2697497/100‑1 (MQ=255)
         gAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGTCCGGGCCTggggct                                                                                            <  1:1957464/100‑4 (MQ=255)
          aGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTag                                                                                           >  2:2664926/1‑100 (MQ=255)
            gAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTagag                                                                                         <  1:260072/100‑1 (MQ=255)
             atctGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAgg                                                                                        <  2:24440/97‑1 (MQ=255)
             aaGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAgg                                                                                        <  2:1714022/100‑1 (MQ=255)
              aGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGc                                                                                       <  2:1909190/100‑1 (MQ=255)
              aGTGAATAATCAGCGCCCCGGGCGCAAAACGTTTACTCATCCCTTTACCAAGCAAGGTAAACGCCAGCAGCAGCACCAGGCCCGGACCTGGGTTAGAgg                                                                                        <  1:1446632/99‑1 (MQ=255)
               gTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGcc                                                                                      <  2:2065647/100‑1 (MQ=255)
               gTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGcc                                                                                      >  1:1925046/1‑100 (MQ=255)
               gTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGcc                                                                                      >  2:1665121/1‑100 (MQ=255)
                 gAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGccac                                                                                    <  2:2885880/100‑1 (MQ=255)
                 gAATAATCATCGCACCGGGCGCAGAACGTTTACTCATCACTTTACCAAACAAGGTAAACGTCAGCAGCAGCCCCAGCCCCGGACCTGGGTTAGAgg                                                                                        >  1:24439/1‑96 (MQ=255)
                  natcaccatcGCCCCGGGAGCAGAACGTTTACTCAAACCTTAACCAAACACGGCAAACGCCAGCAGCAGCCCCAGGCCCGCACCTGGGTTAGAGGccacc                                                                                   >  1:1660612/7‑100 (MQ=255)
                  aataatCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGccacc                                                                                   <  2:1473591/100‑1 (MQ=255)
                   ataatCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGccacca                                                                                  >  1:1001169/1‑100 (MQ=255)
                     aatCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATa                                                                                <  2:3417776/100‑1 (MQ=255)
                      atcatcGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATaa                                                                               <  2:1692477/100‑1 (MQ=255)
                        catcGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAaag                                                                             <  2:533769/100‑1 (MQ=255)
                        catcGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAaag                                                                             >  2:3561838/1‑100 (MQ=255)
                         atcGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAaaga                                                                            >  2:1595439/1‑100 (MQ=255)
                            gCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAagaaga                                                                         >  1:3153524/1‑100 (MQ=255)
                            gCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAagaaga                                                                         >  1:3445540/1‑100 (MQ=255)
                                 gggCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGAt                                                                    <  2:2815197/100‑1 (MQ=255)
                                    cgcAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTa                                                                 <  2:24519/100‑1 (MQ=255)
                                      cAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTAcc                                                               <  1:1590820/100‑1 (MQ=255)
                                       aGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCg                                                              >  2:2644267/1‑100 (MQ=255)
                                       aGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCg                                                              <  2:211568/100‑1 (MQ=255)
                                           cGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTaa                                                          >  2:792738/1‑100 (MQ=255)
                                            gTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAAc                                                         >  1:812714/1‑100 (MQ=255)
                                            gTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAAc                                                         <  1:1899576/100‑1 (MQ=255)
                                             tttACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGGCCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACTGTTAAcc                                                        >  2:3302048/1‑100 (MQ=255)
                                             tttACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAAcc                                                        >  2:570923/1‑100 (MQ=255)
                                             tttACTCATCCCTTTACCAAACAAGGTAAACGACAGCAGCAGACCCAGGCTCGGAGCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAAcc                                                        <  1:1434061/100‑1 (MQ=255)
                                               tACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGa                                                      <  2:2171614/100‑1 (MQ=255)
                                                   cATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGcc                                                  >  2:2774958/1‑100 (MQ=255)
                                                   cATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGcc                                                  >  2:3419669/1‑100 (MQ=255)
                                                    aTCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCt                                                 >  2:3596186/1‑100 (MQ=255)
                                                     tCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTg                                                >  2:2200324/1‑100 (MQ=255)
                                                      cccTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTg                                                >  1:3255624/1‑99 (MQ=255)
                                                       ccTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGtt                                              <  1:522964/100‑1 (MQ=255)
                                                       ccTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGCCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTgat                                              <  2:2578232/100‑3 (MQ=255)
                                                        cTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGACGATGGATTTACCGTTAACCGAAGCCTGTTg                                             <  1:715815/100‑1 (MQ=255)
                                                        cTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTg                                             >  2:1405478/1‑100 (MQ=255)
                                                         tttACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTATAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTg                                             <  1:1003852/99‑1 (MQ=255)
                                                          ttACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGATTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCa                                           >  1:2603143/1‑100 (MQ=255)
                                                           tACCAAACAAGGTAATCGCCAGCAGCAGCCCTAGGTCTGGACCTGGGTTAGAGGCCACTATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCAt                                          <  1:3086484/100‑1 (MQ=255)
                                                            aCCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCAtt                                         >  1:3752651/1‑100 (MQ=255)
                                                             ccAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGccagatcggaag                                        >  2:1922271/1‑90 (MQ=255)
                                                             ccAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTc                                        >  2:760417/1‑100 (MQ=255)
                                                               aaaCAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTccc                                      >  2:1742709/1‑100 (MQ=255)
                                                               aaaCAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTccc                                      >  2:3656789/1‑100 (MQ=255)
                                                                  cAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGc                                   >  2:3337679/1‑100 (MQ=255)
                                                                     ggTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGa                                >  1:3367995/1‑100 (MQ=255)
                                                                      gTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGat                               <  2:3035541/100‑1 (MQ=255)
                                                                      gTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGan                               <  2:3142548/100‑2 (MQ=255)
                                                                          aCGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGataata                           <  1:3647518/100‑1 (MQ=255)
                                                                           cGCAAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGCCCACCATAAAGAAGAGGGATTTACCATTAACCGAAGCCTGTTGCATTCCCAGCGGataata                           <  1:629993/99‑1 (MQ=255)
                                                                            gCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGa                         <  1:851180/100‑1 (MQ=255)
                                                                             ccAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGAc                        >  2:3817711/1‑100 (MQ=255)
                                                                                gcagcagcCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGc                      >  1:3054037/1‑99 (MQ=255)
                                                                                 cctctaacccAGGTCCGGGCCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGccc                    >  2:1734340/8‑100 (MQ=255)
                                                                                 cctctaacccAGGTCCGGGCCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGccc                    >  2:1097653/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                   >  1:390839/7‑100 (MQ=255)
                                                                                  agcagcCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGccc                    >  1:347004/1‑99 (MQ=255)
                                                                                   tctaacccAGGTCCGGGCCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTg                  >  2:3019881/6‑100 (MQ=255)
                                                                                   tctaacccAGGTCCGGGCCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTg                  >  2:209116/6‑100 (MQ=255)
                                                                                    ctaacccAGGTCCGGGCCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTga                 >  2:2948922/5‑100 (MQ=255)
                                                                                    ctaacccAGGTCCGGGCCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCACTga                 >  2:591162/5‑100 (MQ=255)
                                                                                     agcCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTgat                >  2:3695865/1‑100 (MQ=255)
                                                                                      gcCCCAGGCCCGGGCCTGGGTTAGAGGCCACCATAAAGTAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTgatc               >  1:2352559/1‑100 (MQ=255)
                                                                                      gcCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTgatc               >  2:903184/1‑100 (MQ=255)
                                                                                      aacccAGGTCCGGGCCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTgatc               >  1:295562/3‑100 (MQ=255)
                                                                                      aacccAGGTCCGGGCCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTgatc               >  2:641058/3‑100 (MQ=255)
                                                                                      aacccAGGTCCGGGCCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTgatc               >  1:2135437/3‑100 (MQ=255)
                                                                                      aacccAGGTCCGGGCCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTga                 >  1:295610/3‑98 (MQ=255)
                                                                                       ccccAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTgatcg              >  2:407925/1‑100 (MQ=255)
                                                                                        cccAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTgatcga             >  1:1535922/1‑100 (MQ=255)
                                                                                        cccAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTgatcga             >  2:671212/1‑100 (MQ=255)
                                                                                          cAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTgatcgatc           >  2:677655/1‑100 (MQ=255)
                                                                                            ggCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGc         <  1:760409/100‑1 (MQ=255)
                                                                                            ggCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGc         >  1:1375060/1‑100 (MQ=255)
                                                                                            ggCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGc         <  2:811088/100‑1 (MQ=255)
                                                                                             gCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCCCa        <  2:2233850/100‑1 (MQ=255)
                                                                                              cccGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCat       >  2:3535192/1‑100 (MQ=255)
                                                                                               ccGGGCCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCat       <  1:1097646/99‑1 (MQ=255)
                                                                                                   aCCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGCAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTg  <  1:2994454/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTG  >  NC_000913/3078092‑3078290

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: