Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,078,191 A→G 7.3% G223G (GGT→GGC cmtA ← mannitol‑specific PTS enzyme IICB component CmtA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,078,1910AG 7.3% 234.5 / 6.7 82G223G (GGT→GGCcmtAmannitol‑specific PTS enzyme IICB component CmtA
Reads supporting (aligned to +/- strand):  ref base A (38/38);  new base G (3/3);  total (41/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

AATGGTCAGCGGCTTCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATTGATGACTGACAACAATGGCAG  >  NC_000913/3078047‑3078337
                                                                                                                                                |                                                                                                                                                  
aatGGTCAGCGGCTTCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTgg                                                                                                                                               >  1:1090896/1‑150 (MQ=255)
aatGGTCAGCGGCTTCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTgg                                                                                                                                               <  2:1902354/150‑1 (MQ=255)
     tCAGCGGCTTCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGTCCGGGCCTggggctg                                                                                                                                          <  1:931880/150‑5 (MQ=255)
     tCAGCGGCTTCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGTCCGGGCCTggggctg                                                                                                                                          <  1:261217/150‑5 (MQ=255)
        gCGGCTTCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAgg                                                                                                                                       >  2:618338/1‑150 (MQ=255)
          ggCTTCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGcc                                                                                                                                     >  1:998929/1‑150 (MQ=255)
          ggCTTCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGcc                                                                                                                                     <  1:801576/150‑1 (MQ=255)
            cTTCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGccac                                                                                                                                   >  1:1925963/1‑150 (MQ=255)
             ttCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGccacc                                                                                                                                  <  2:1090896/150‑1 (MQ=255)
             ttCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGccacc                                                                                                                                  <  2:1721078/150‑1 (MQ=255)
                 tcaGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATaa                                                                                                                              <  2:1081986/150‑1 (MQ=255)
                  caGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATaaa                                                                                                                             >  2:1099903/1‑150 (MQ=255)
                  caGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATaaa                                                                                                                             >  2:1782922/1‑150 (MQ=255)
                    gCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTag                                                                                                                                          >  1:937505/1‑135 (MQ=255)
                    gCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTag                                                                                                                                          <  2:937505/135‑1 (MQ=255)
                     cacaTACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAaagaa                                                                                                                          <  1:1824719/150‑1 (MQ=255)
                         tACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATg                                                                                                                      >  2:1532828/1‑150 (MQ=255)
                           cGGGAAATACAGTTCGTGGATACCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGa                                                                                                                    >  2:811993/1‑150 (MQ=255)
                             gtaaaTAAAGTTAGTTGATCCAACCGAGTAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAAAAAGGTAAACGCCATCAGCAGCCCCAGGCCAGGACCTGGGTTAGAGGCCATAATAAAGAAGATGGAtt                                                                                                                  <  2:515943/148‑1 (MQ=255)
                                    cAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAgg                                                                                                                                       <  2:314769/122‑1 (MQ=255)
                                    cAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAgg                                                                                                                                       >  1:314769/1‑122 (MQ=255)
                                    cAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGtt                                                                                                           >  2:1489985/1‑150 (MQ=255)
                                      gTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTaa                                                                                                         <  1:1290900/150‑1 (MQ=255)
                                         cGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACTAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCg                                                                                                      >  2:1833193/1‑150 (MQ=255)
                                           tGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGaa                                                                                                    >  2:207765/1‑150 (MQ=255)
                                            ggATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAg                                                                                                   <  2:1187882/150‑1 (MQ=255)
                                                   aCCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTg                                                                                            <  1:436736/150‑1 (MQ=255)
                                                    ccGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGc                                                                                           >  2:1470423/1‑150 (MQ=255)
                                                      gAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCAt                                                                                         <  1:811993/150‑1 (MQ=255)
                                                         gAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTcc                                                                                      <  1:57178/150‑1 (MQ=255)
                                                           aGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCa                                                                                    >  2:1112119/1‑150 (MQ=255)
                                                               aataatCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCgg                                                                                <  1:933407/150‑1 (MQ=255)
                                                                 taatCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGat                                                                              >  2:1849279/1‑150 (MQ=255)
                                                                  aatCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGata                                                                             <  1:1489985/150‑1 (MQ=255)
                                                                  aatCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGata                                                                             >  1:1184497/1‑150 (MQ=255)
                                                                    tcatcGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGataat                                                                           <  2:485734/150‑1 (MQ=255)
                                                                     catcGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGataata                                                                          >  1:1425185/1‑150 (MQ=255)
                                                                      atcGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAg                                                                         <  1:1475985/150‑1 (MQ=255)
                                                                      atcGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAg                                                                         >  1:652334/1‑150 (MQ=255)
                                                                          ccccGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGc                                                                     >  1:356699/1‑150 (MQ=255)
                                                                          ccccGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGc                                                                     >  1:98004/1‑150 (MQ=255)
                                                                             cGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCt                                                                  <  1:146245/150‑1 (MQ=255)
                                                                               ggCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTga                                                                >  2:1736146/1‑150 (MQ=255)
                                                                                   cAGAAAGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTgatcga                                                            >  1:481095/1‑150 (MQ=255)
                                                                                    aGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTgatcgat                                                           <  2:1148861/150‑1 (MQ=255)
                                                                                       aCGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGc                                                        <  1:1354886/150‑1 (MQ=255)
                                                                                       aCGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGc                                                        >  2:1354886/1‑150 (MQ=255)
                                                                                         gTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCat                                                      >  1:1918901/1‑150 (MQ=255)
                                                                                              cTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTg                                                 >  2:245462/1‑150 (MQ=255)
                                                                                                cATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAg                                               <  1:401795/150‑1 (MQ=255)
                                                                                                  tCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGa                                             >  1:253929/1‑150 (MQ=255)
                                                                                                    ccTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGaaa                                           >  1:1905775/1‑150 (MQ=255)
                                                                                                      tttACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAg                                         <  2:761374/150‑1 (MQ=255)
                                                                                                        tACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCa                                       <  2:1389207/150‑1 (MQ=255)
                                                                                                         aCCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCAc                                      >  2:1328044/1‑150 (MQ=255)
                                                                                                              aCAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCg                                 >  1:1355987/1‑150 (MQ=255)
                                                                                                              aCAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGGTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCg                                 >  1:1349121/1‑150 (MQ=255)
                                                                                                                 aGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCg                              >  1:1914996/1‑150 (MQ=255)
                                                                                                                   gTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCTGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGt                            <  2:221444/150‑1 (MQ=255)
                                                                                                                   gTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGt                            <  2:1425185/150‑1 (MQ=255)
                                                                                                                        cGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCAtt                       <  1:1736146/150‑1 (MQ=255)
                                                                                                                        cGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCAtt                       <  2:1184497/150‑1 (MQ=255)
                                                                                                                        cGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCAAAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCAtt                       >  2:1363754/1‑150 (MQ=255)
                                                                                                                             gcagcagcCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATtgatga                  <  1:503819/150‑1 (MQ=255)
                                                                                                                              cctctaacccAGGTCCGGGCCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATTGAtgac                 >  2:1498230/8‑150 (MQ=255)
                                                                                                                              cagcagcCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATTGAtgac                 >  1:1003925/1‑150 (MQ=255)
                                                                                                                              cagcagcCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATTGAtgac                 <  1:207765/150‑1 (MQ=255)
                                                                                                                               ctctaacccAGGTCCGGGCCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTGACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGCTCGATAGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATTGAtgact                >  2:1799357/7‑150 (MQ=255)
                                                                                                                                 ctaacccAGGTCCGGGCCTGGGTTAGAGGCCACAATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATTGAtgactga              >  2:1430646/5‑150 (MQ=255)
                                                                                                                                  agcCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATTGAtgactgac             <  1:40797/150‑1 (MQ=255)
                                                                                                                                  agcCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATTGAtgactgac             <  1:184198/150‑1 (MQ=255)
                                                                                                                                  agcCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATTGAtgactgac             <  2:1918901/150‑1 (MQ=255)
                                                                                                                                   gcCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGCCGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATTGATGACTGaca            >  1:1520625/1‑150 (MQ=255)
                                                                                                                                   gcCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATTGATGACTGaca            <  2:1349121/150‑1 (MQ=255)
                                                                                                                                   gcCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATTGATGACTGaca            <  2:1355987/150‑1 (MQ=255)
                                                                                                                                     cccAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATTGATGACTGacaac          >  1:1511882/1‑150 (MQ=255)
                                                                                                                                       cAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATTGATGACTGAcaacaa        >  2:742185/1‑150 (MQ=255)
                                                                                                                                        aGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATTGATGACTGACAACAAt       <  1:1849279/150‑1 (MQ=255)
                                                                                                                                        aGGCCCGGACCTGCGTTAGAGGCAACGATAAAGAAGAAGGATTTACCGTAAACCGAAGCCTGTTGCATTCCCAGCGCATAAAAGACGCCCTCATCCATCCCATAATTCAGGAAAAGCCCTTTCGCCGCTTCATTGATGAATGACAACAAt       >  2:969474/1‑150 (MQ=255)
                                                                                                                                         ggCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATTGATGACTGACAACAATg      <  1:618338/150‑1 (MQ=255)
                                                                                                                                            ccGGGCCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATTGATGACTGACAACAATGGCa   <  1:1498230/150‑1 (MQ=255)
                                                                                                                                             cGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATTGATGACTGACAACAATGGCAg  <  2:98004/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
AATGGTCAGCGGCTTCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATTGATGACTGACAACAATGGCAG  >  NC_000913/3078047‑3078337

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