Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,593,935 G→C 8.0% A310P (GCG→CCG)  PP_3171 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,593,9350GC 8.0% 226.9 / 8.6 75A310P (GCG→CCG) PP_3171hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (33/36);  new base C (3/3);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

CGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTCCTCTCACACCACCGTGCGTACGGTTCCGTACACGGCGGTTCAGGTTACGCGGCTAAGCCG  >  NC_002947/3593789‑3594084
                                                                                                                                                  |                                                                                                                                                     
cGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGgcgc                                                                                                                                                    <  2:1002304/150‑1 (MQ=255)
cGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGgcgc                                                                                                                                                    >  1:105482/1‑150 (MQ=255)
cGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGgcgc                                                                                                                                                    <  2:319976/150‑1 (MQ=255)
cGGACCACGAGCGCGCCGTGCTCGATCAGGGGATGGCCTGCTCCGTTGCTGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGCGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGgcgc                                                                                                                                                    <  1:674671/150‑1 (MQ=255)
   aCCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGgcgc                                                                                                                                                    >  1:491657/1‑147 (MQ=255)
   aCCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGgcgc                                                                                                                                                    <  2:491657/147‑1 (MQ=255)
     cACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAg                                                                                                                                               >  2:523139/1‑150 (MQ=255)
              gcCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGCCCCGGGGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACAt                                                                                                                                      <  1:87742/150‑1 (MQ=255)
              gcCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACAt                                                                                                                                      >  2:924222/1‑150 (MQ=255)
              gcCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACAt                                                                                                                                      <  2:870827/150‑1 (MQ=255)
                      cGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGc                                                                                                                              <  2:1178679/150‑1 (MQ=255)
                       gATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGcc                                                                                                                             <  1:577834/150‑1 (MQ=255)
                         tCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCcgc                                                                                                                           >  2:298275/1‑150 (MQ=255)
                             gTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAg                                                                                                                       >  2:89658/1‑150 (MQ=255)
                                 tGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGAAGCGGTCACAGGCCTATGCCGCTCAGgcgc                                                                                                                   >  2:1017033/1‑150 (MQ=255)
                                  ggCCTGCTCCGTAGCCGGAGGGCCAGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCgct                                                                                                                  <  2:233512/150‑1 (MQ=255)
                                   gCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGAACGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCgctg                                                                                                                 <  1:990237/150‑1 (MQ=255)
                                    ccTGCTCCGTAGCCGGGGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCgctgg                                                                                                                >  1:409363/1‑150 (MQ=255)
                                      tGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCgctggct                                                                                                              >  2:924620/1‑150 (MQ=255)
                                         tCCGTAGCCGGAGTGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGaa                                                                                                           <  1:291460/150‑1 (MQ=255)
                                         tCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGaa                                                                                                           >  2:1421809/1‑150 (MQ=255)
                                            gTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTg                                                                                                        <  2:50056/150‑1 (MQ=255)
                                             tAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACTCGCGGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTgg                                                                                                       <  2:1058675/150‑1 (MQ=255)
                                               gCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGc                                                                                                     <  1:1024508/150‑1 (MQ=255)
                                               gCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCGGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGc                                                                                                     <  1:612190/150‑1 (MQ=255)
                                                ccGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCt                                                                                                    <  1:131923/150‑1 (MQ=255)
                                                ccGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGCCGCGCCAGCgg                                                                                                                                            >  2:1042558/1‑110 (MQ=255)
                                                ccGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGCCGCGCCAGCgg                                                                                                                                            <  1:1042558/110‑1 (MQ=255)
                                                 cGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTg                                                                                                   >  1:1420394/1‑150 (MQ=255)
                                                   gAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGac                                                                                                 <  2:105482/150‑1 (MQ=255)
                                                    aGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGaca                                                                                                <  2:1063658/150‑1 (MQ=255)
                                                         ccGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATc                                                                                           >  2:693399/1‑150 (MQ=255)
                                                         ccGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATc                                                                                           <  1:1146295/150‑1 (MQ=255)
                                                         ccGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATc                                                                                           >  1:445449/1‑150 (MQ=255)
                                                         ccGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGCCGCGCCAGCgg                                                                                                                                            <  2:151136/101‑1 (MQ=255)
                                                         ccGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGCCGCGCCAGCgg                                                                                                                                            >  1:151136/1‑101 (MQ=255)
                                                           gAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATGGa                                                                                         <  1:1421809/150‑1 (MQ=255)
                                                           gAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGa                                                                                         <  2:1005773/150‑1 (MQ=255)
                                                                 gTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGAGACGGTCGAGTaggg                                                                                   >  2:865138/1‑149 (MQ=255)
                                                                 gTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTagga                                                                                   >  2:1390754/1‑150 (MQ=255)
                                                                   gAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTaggagg                                                                                 <  1:391891/150‑1 (MQ=255)
                                                                    aGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGc                                                                                <  2:302157/150‑1 (MQ=255)
                                                                      ggAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGa                                                                                             <  1:947262/135‑1 (MQ=255)
                                                                      ggAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGa                                                                                             >  2:947262/1‑135 (MQ=255)
                                                                      ggAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCgg                                                                              >  1:575846/1‑150 (MQ=255)
                                                                      ggAAGGGCTGCGCTCGCTGATCGAGCAAACCGGGGCGGAACAGCGGCTGATAGAGTGCCGCATCCATGCCCCGCTGGCACGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCgg                                                                              >  1:574852/1‑150 (MQ=255)
                                                                       gAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGa                                                                             <  1:1279342/150‑1 (MQ=255)
                                                                            gCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTcac                                                                        <  1:924620/150‑1 (MQ=255)
                                                                              tGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGa                                                                                         >  1:1136853/1‑131 (MQ=255)
                                                                              tGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGa                                                                                         <  2:1136853/131‑1 (MQ=255)
                                                                                     gctGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCgctggct                                                                                                              >  2:983563/1‑103 (MQ=255)
                                                                                     gctGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCgctggct                                                                                                              <  1:983563/103‑1 (MQ=255)
                                                                                     gctGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTc                                                                          >  1:483543/1‑139 (MQ=255)
                                                                                     gctGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTc                                                                          <  2:483543/139‑1 (MQ=255)
                                                                                                       ggCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTCCTCTCACACCACCGTGc                                             >  2:563157/1‑150 (MQ=255)
                                                                                                              gAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAg                                                                                        >  1:1215124/1‑100 (MQ=255)
                                                                                                              gAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAg                                                                                        <  2:1215124/100‑1 (MQ=255)
                                                                                                              gAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTcaca                                                                       <  2:1364187/117‑1 (MQ=255)
                                                                                                              gAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTcaca                                                                       >  1:1364187/1‑117 (MQ=255)
                                                                                                               aGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGCCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTCCTCTCACACCACCGTGCGTACGGtt                                     <  1:89658/150‑1 (MQ=255)
                                                                                                                  tgatgaTCGACTGCCACATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATGCGCACCCGCCGTCCTCTCACACCACCGTGCGTACGGTTCCg                                  >  1:727688/1‑150 (MQ=255)
                                                                                                                     tgatCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTCCTCTCACACCACCGTGCGTACGGTTCCGTac                               >  2:597701/1‑150 (MQ=255)
                                                                                                                      gatCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGc                                                                                >  2:957228/1‑100 (MQ=255)
                                                                                                                      gatCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGc                                                                                <  1:957228/100‑1 (MQ=255)
                                                                                                                       atCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACAt                                                                                                                                      <  1:489928/45‑1 (MQ=255)
                                                                                                                       atCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACAt                                                                                                                                      >  2:489928/1‑45 (MQ=255)
                                                                                                                            cTGCCGCCTCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATc                                                                                           >  2:420499/1‑83 (MQ=255)
                                                                                                                            cTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATc                                                                                           <  1:420499/83‑1 (MQ=255)
                                                                                                                                     ccATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTCCTCTCACACCACCGTGCGTACGGTTCCGTACACGGCGGTTCAGGTTa               >  2:1343889/1‑150 (MQ=255)
                                                                                                                                      cATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTCCTCTCACACCACCGTGCGTACGGTTCCGTACACGGCGGTTCAGGTTAc              >  1:1248930/1‑150 (MQ=255)
                                                                                                                                       aTGCCCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGt                                                                                       <  1:1387339/76‑1 (MQ=255)
                                                                                                                                       aTGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGt                                                                                       >  2:1387339/1‑76 (MQ=255)
                                                                                                                                          gaccGCTGCCGCGCCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTc                                                                          >  1:946187/3‑86 (MQ=255)
                                                                                                                                          gaccGCTGCCGCGCCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTc                                                                          <  2:946187/84‑1 (MQ=255)
                                                                                                                                          aCCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTCCTCTCACACCACCGTGCGTACGGTTCCGTACACGGCGGTTCAGGTTACGCgg          <  1:523139/150‑1 (MQ=255)
                                                                                                                                            ccGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGa                                                                                         <  2:323313/69‑1 (MQ=255)
                                                                                                                                            ccGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGa                                                                                         >  1:323313/1‑69 (MQ=255)
                                                                                                                                            ccGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTCCTCTcac                                                       <  2:792021/103‑1 (MQ=255)
                                                                                                                                            ccGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTCCTCTcac                                                       >  1:792021/1‑103 (MQ=255)
                                                                                                                                            ccGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTCCTCTCACACCACCGTGCGTACGGTTCCGTa                                <  2:1099762/126‑1 (MQ=255)
                                                                                                                                            ccGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTCCTCTCACACCACCGTGCGTACGGTTCCGTa                                >  1:1099762/1‑126 (MQ=255)
                                                                                                                                                 ggCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTCCTCTCACACCACCGTGCGTACGGTTCCGTACACGGCGGTTCAGGTTACGCGGCTAAGcc   <  1:597701/150‑1 (MQ=255)
                                                                                                                                                  gcgcGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTCCTCTCACACCACCGTGCGTACGGTTCCGTACACGGCGGTTCAGGTTACGCGGCTAAGCCg  <  1:1343889/150‑1 (MQ=255)
                                                                                                                                                  gcgcGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTCCTCTCACACCACCGTGCGTACGGTTCCGTACACGGCGGTTCAGGTTACGCGGCTAAGCCg  <  1:924222/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGGTCACATGCCTATGCCGCGCAGGCGCTGGCTGAACTGGGCTGACACGATCGAGTAGGAGGCGGATTCACACCCGCCGTCCTCTCACACCACCGTGCGTACGGTTCCGTACACGGCGGTTCAGGTTACGCGGCTAAGCCG  >  NC_002947/3593789‑3594084

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