Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,863,398 Δ1 bp 5.6% coding (1042/1068 nt) PP_1666 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,863,3980T. 5.6% 227.9 / 7.1 72coding (1042/1068 nt)PP_1666hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (29/39);  new base . (2/2);  total (31/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 = 5.20e-01

GCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTCCTGGTCGGCATCATGCTGGC  >  NC_002947/1863249‑1863546
                                                                                                                                                     |                                                                                                                                                     
gCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGt                                                                                                                                                       <  1:907126/150‑1 (MQ=255)
  aaTCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTc                                                                                                                                                     >  2:700007/1‑150 (MQ=255)
          aCAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCtt                                                                                                                                            <  1:770954/150‑1 (MQ=255)
            aGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGc                                                                                                                                          >  1:890391/1‑150 (MQ=255)
              cTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGt                                                                                                                                        >  1:3489/1‑150 (MQ=255)
                    gCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGc                                                                                                                                  <  2:561492/150‑1 (MQ=255)
                     cccAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCt                                                                                                                                 >  2:804630/1‑150 (MQ=255)
                     cccAGCTTGGCCTCGCCAAGCTGCAGGAAGTGCCCGTCGAGCAGGCACCGGCGGTGCTGGCCACGCCAGTACCGGCAGTCGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCt                                                                                                                                 <  2:265282/150‑1 (MQ=255)
                      ccAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTg                                                                                                                                <  1:402089/150‑1 (MQ=255)
                        aGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGt                                                                                                                              <  2:890391/150‑1 (MQ=255)
                          cTTGGCCAGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGc                                                                                                                                              >  1:1087847/1‑132 (MQ=255)
                          cTTGGCCAGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGc                                                                                                                                              <  2:1087847/132‑1 (MQ=255)
                             ggcctggccAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGa                                                                                                                         <  2:969016/150‑1 (MQ=255)
                             ggcctggccAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGa                                                                                                                         <  2:434365/150‑1 (MQ=255)
                              gcctggccAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAg                                                                                                                        >  2:700356/1‑150 (MQ=255)
                                 tggccAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAgcgc                                                                                                                     <  2:237048/150‑1 (MQ=255)
                                      aaGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGcc                                                                                                                <  2:726127/150‑1 (MQ=255)
                                          tGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACGCAAGCCGTCGAACGGCTTgggt                                                                                                                                        >  2:82422/1‑119 (MQ=255)
                                          tGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACGCAAGCCGTCGAACGGCTTGGGt                                                                                                                                        <  1:82422/122‑1 (MQ=255)
                                          tGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGtga                                                                                                            <  2:994768/150‑1 (MQ=255)
                                           gCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCAGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGtgac                                                                                                           >  2:1006767/1‑150 (MQ=255)
                                              ggAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGttt                                                                                                                                      >  1:566397/1‑120 (MQ=255)
                                              ggAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGttt                                                                                                                                      <  2:566397/120‑1 (MQ=255)
                                              ggAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGtgactga                                                                                                        <  1:700356/150‑1 (MQ=255)
                                               gAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGtgactgac                                                                                                       >  2:35768/1‑150 (MQ=255)
                                                aaGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACa                                                                                                      <  2:246760/150‑1 (MQ=255)
                                                       ggTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCg                                                                                               >  1:1261982/1‑150 (MQ=255)
                                                          cGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTg                                                                                            <  1:700007/150‑1 (MQ=255)
                                                           gAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTc                                                                                                                                                     >  1:1002995/1‑93 (MQ=255)
                                                           gAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTc                                                                                                                                                     <  2:1002995/93‑1 (MQ=255)
                                                           gAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTgg                                                                                           >  1:539285/1‑150 (MQ=255)
                                                             gccggCACCGCCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCCCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCCCTTGCCGTGACTGACATGCGTCCCTGGAt                                                                                         >  2:661229/4‑150 (MQ=255)
                                                             gCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGAt                                                                                         >  1:983333/1‑150 (MQ=255)
                                                                 gCACCGGCGGTGCCGGCCATGCCAGTACCGTCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTg                                                                                                                                >  1:1534024/1‑107 (MQ=255)
                                                                 gCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTg                                                                                                                                <  2:1534024/107‑1 (MQ=255)
                                                                     cggcggTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATctggctgg                                                                                 <  1:499195/150‑1 (MQ=255)
                                                                               ggCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCAct                                                                       >  1:907134/1‑150 (MQ=255)
                                                                                  cATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCActgct                                                                    <  2:539285/150‑1 (MQ=255)
                                                                                   aTGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCActgctg                                                                   <  2:1261982/150‑1 (MQ=255)
                                                                                      ccAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATc                                                                >  1:498694/1‑150 (MQ=255)
                                                                                       cAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTg                                                                                                                  <  1:1094205/99‑1 (MQ=255)
                                                                                       cAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTg                                                                                                                  >  2:1094205/1‑99 (MQ=255)
                                                                                         gTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCt                                                             <  1:223682/150‑1 (MQ=255)
                                                                                          tACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCtg                                                            <  2:614328/150‑1 (MQ=255)
                                                                                            ccGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCtgtg                                                          <  2:3489/150‑1 (MQ=255)
                                                                                             cGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGtgc                                                         >  1:420509/1‑150 (MQ=255)
                                                                                                 aGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACACGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTgctgct                                                     <  1:8041/150‑1 (MQ=255)
                                                                                                        gccgcAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGc                                              >  2:1120800/1‑150 (MQ=255)
                                                                                                        gccgcAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGc                                              >  1:602529/1‑150 (MQ=255)
                                                                                                          cgcAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTaa                                                                                                                            >  1:437902/1‑70 (MQ=255)
                                                                                                          cgcAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTaa                                                                                                                            <  2:437902/70‑1 (MQ=255)
                                                                                                          cgcAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCt                                            <  2:498694/150‑1 (MQ=255)
                                                                                                           gcAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGAGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTg                                           >  2:279659/1‑150 (MQ=255)
                                                                                                                 gcgcAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAat                                     <  1:378384/150‑1 (MQ=255)
                                                                                                                 gcgcAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAat                                     >  1:762221/1‑150 (MQ=255)
                                                                                                                   gcAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAatat                                   <  1:1120800/150‑1 (MQ=255)
                                                                                                                        ccGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTa                              >  1:1309968/1‑150 (MQ=255)
                                                                                                                        ccGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAAGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTa                              >  1:245637/1‑150 (MQ=255)
                                                                                                                         cgcgCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTAc                             <  2:550206/150‑1 (MQ=255)
                                                                                                                            gcCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGcc                          <  1:1006767/150‑1 (MQ=255)
                                                                                                                            gcCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGcc                          <  2:383679/150‑1 (MQ=255)
                                                                                                                            gcCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGcc                          >  1:975192/1‑150 (MQ=255)
                                                                                                                                gcgcAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGAGCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTTCGCCGGTc                      <  2:563817/150‑1 (MQ=255)
                                                                                                                                 ccatctcagcACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTCCTTCCTGCACAATATTCTTACGCCGTTcc                     <  1:43259/144‑1 (MQ=255)
                                                                                                                                  gcAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGAACTGCTGTTCGCTGTCCTGCTATATACCCTTGACAATATTCTTACGCCGTTCCt                    <  1:35768/150‑1 (MQ=255)
                                                                                                                                   cagccagcACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTCCTg                   <  2:762221/150‑1 (MQ=255)
                                                                                                                                        agcACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTATTTGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTCCtagatcg              >  2:43259/1‑144 (MQ=255)
                                                                                                                                         gcACCCAAGCCGTTCGACGGCTTGCGTCTTCTCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTCCTGGTCGGc             <  1:661229/150‑1 (MQ=255)
                                                                                                                                         aaacgcAAGCCGTCGAACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTCCTGGTCGGc             >  2:581475/6‑150 (MQ=255)
                                                                                                                                           aCCCAAGCCGTTCGACGGCTTGTGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAata                                    >  2:1038698/1‑125 (MQ=255)
                                                                                                                                           aCCCAAGCCGTTCGACGGCTTGTGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAata                                    <  1:1038698/125‑1 (MQ=255)
                                                                                                                                            cgcAAGCCGTCGAACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTCCTGGTCGGcatc          <  1:581475/148‑1 (MQ=255)
                                                                                                                                                    gTTCGACGTCGTTCGTTTTCGCCGATCAGCCGCCATTCCCGCGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTCCTGGTCGGCATCATGCTGGc  <  2:676274/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTACCGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGAGCGCTTGCCGTGACTGACATGCGTCGCTGGATCTGGCTGGGCGCTGCACTGCTGATCGCTGTGCTGCTTTATAGCCTGCACAATATTCTTACGCCGTTCCTGGTCGGCATCATGCTGGC  >  NC_002947/1863249‑1863546

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