Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 295,250 T→C 5.2% L218P (CTC→CCC)  ssuB → aliphatic sulfonates ABC transporter ATP‑binding protein SsuB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947295,2500TC 5.2% 609.6 / 13.5 190L218P (CTC→CCC) ssuBaliphatic sulfonates ABC transporter ATP‑binding protein SsuB
Reads supporting (aligned to +/- strand):  ref base T (99/80);  new base C (7/3);  total (107/83)
Fisher's exact test for biased strand distribution p-value = 5.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.26e-01

Only 100 of 192 total aligned reads displayed.

GATGCGTTGACCCGTATCGAGATGCAGCAACTGATCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCGGGCACTGCGCCCGAGCCGGATCCTGTAGCCCCTCTGCCCACGCAACTGCGCT  >  NC_002947/295102‑295399
                                                                                                                                                    |                                                                                                                                                     
gATGCGTTGACCCGTATCGAGATGCAGCAACTGATCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTc                                                                                                                                                      <  2:107288/150‑1 (MQ=255)
 aTGCGTTGACCCGTATCGAGATGCAGCAACTGATCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCg                                                                                                                                                     <  1:5373221/150‑1 (MQ=255)
  tGCGTTGACCCGTATCGAGATGCAGCAACTGATCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGa                                                                                                                                                    <  1:3001909/150‑1 (MQ=255)
    cGTTGACCCGTATCGAGATGCAGCAACTGATCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGAcc                                                                                                                                                  <  2:4967921/150‑1 (MQ=255)
     gTTGACCCGTATCGAGATGCAGCAACTGATCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCt                                                                                                                                                 >  1:3112455/1‑150 (MQ=255)
         aCCCGTATCGAGATGCAGCAACTGATCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACt                                                                                                                                             <  1:56172/150‑1 (MQ=255)
         aCCCGTATCGAGATGCAGCAACTGATCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACAGTGAGGTCGAGCCCGACCTcgcc                                                                                                                                             <  2:719242/150‑4 (MQ=255)
           ccttATCGAGATGCAGCAACAGATCGAACGCCTGTCGCGTAAGCACGGCATCACTTTACTGCTCGTCAACCACGACCTCAGCGAGTCCGTTGCCGTGGCTGACAGGGAGATCCTGAACGAGGACGGCTAGGTCGGGCTCGACCTCACTGt                                                                                                                                           >  2:3080088/4‑150 (MQ=255)
           ccGTATCGAGATGCAGCAACTGATCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGt                                                                                                                                           >  1:3493452/1‑150 (MQ=255)
                 cGAGATGCAGCAACTGATCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCt                                                                                                                                     >  2:1570618/1‑150 (MQ=255)
                 cGAGATGCAGCAACTGATCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCt                                                                                                                                     >  2:1348765/1‑150 (MQ=255)
                  gagaTGCAGCAACTGATCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTg                                                                                                                                    >  1:1055260/1‑150 (MQ=255)
                  gagaTGCAGCAACTGATCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTg                                                                                                                                    <  1:1868640/150‑1 (MQ=255)
                  gagaTGCAGCAACTGATCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTg                                                                                                                                    <  2:1612481/150‑1 (MQ=255)
                  gagaTGCAGCAACTGATCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTAGACCTCACTGTTGACCTg                                                                                                                                    <  1:3151009/150‑1 (MQ=255)
                           caACTGATCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGccc                                                                                                                           >  1:3016852/1‑150 (MQ=255)
                                gATCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCTGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGgc                                                                                                                      <  2:1109262/150‑1 (MQ=255)
                                gATCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGgc                                                                                                                      >  2:4187998/1‑150 (MQ=255)
                                gATCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGgc                                                                                                                      >  1:4071207/1‑150 (MQ=255)
                                 aTCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGTGAGGTCGAGCCCGACCTcgcc                                                                                                                                             <  2:780584/126‑4 (MQ=255)
                                  tCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGACTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGgcgc                                                                                                                    <  2:3385926/150‑1 (MQ=255)
                                     aaCGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTg                                                                                                                 >  1:232305/1‑150 (MQ=255)
                                          cTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCg                                                                                                            >  1:977649/1‑150 (MQ=255)
                                          cTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCg                                                                                                            <  2:2619586/150‑1 (MQ=255)
                                          cTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCg                                                                                                            >  1:3548480/1‑150 (MQ=255)
                                             tgGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCAc                                                                                                         <  2:3805762/150‑1 (MQ=255)
                                                 gTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTc                                                                                                     >  2:935412/1‑150 (MQ=255)
                                                 gTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTc                                                                                                     >  1:5270071/1‑150 (MQ=255)
                                                  tCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCt                                                                                                    >  1:4860473/1‑150 (MQ=255)
                                                  tCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCt                                                                                                    >  1:2561582/1‑150 (MQ=255)
                                                  tCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCGGGCGCGTGGTTCGCACCggggc                                                                                                   <  1:3948228/150‑5 (MQ=255)
                                                    aGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTgg                                                                                                  <  1:2916322/150‑1 (MQ=255)
                                                    aGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTgg                                                                                                  <  1:5425348/150‑1 (MQ=255)
                                                     gCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGc                                                                                                 >  1:67349/1‑150 (MQ=255)
                                                      cACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGcc                                                                                                >  2:3176141/1‑150 (MQ=255)
                                                       aCGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCcg                                                                                               >  2:529893/1‑150 (MQ=255)
                                                       aCGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCcg                                                                                               >  1:5127054/1‑150 (MQ=255)
                                                         ggCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCcgcg                                                                                             >  1:750702/1‑150 (MQ=255)
                                                         ggCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCcgcg                                                                                             >  2:5340309/1‑150 (MQ=255)
                                                            ttCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTg                                                                                          >  1:1783678/1‑150 (MQ=255)
                                                            ttCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTg                                                                                          <  2:525876/150‑1 (MQ=255)
                                                            gtCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGCCCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTg                                                                                          <  2:712573/149‑1 (MQ=255)
                                                             tCACTGTACGGCTGGTCACCCACGACGTCAGCGCGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTgg                                                                                         <  1:5209922/150‑1 (MQ=255)
                                                               aCTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGaa                                                                                       >  2:118051/1‑150 (MQ=255)
                                                                cTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGaaa                                                                                      <  2:3514992/150‑1 (MQ=255)
                                                                cTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGaaa                                                                                      <  1:4128259/150‑1 (MQ=255)
                                                                 tGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAg                                                                                     <  1:3640365/150‑1 (MQ=255)
                                                                   tACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGTCGCGCTGGAAAGCg                                                                                   >  1:448178/1‑150 (MQ=255)
                                                                     ctgctgGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGaa                                                                                 <  1:3802472/150‑1 (MQ=255)
                                                                      tgctgGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAg                                                                                >  1:218180/1‑150 (MQ=255)
                                                                         tgGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGAc                                                                                                                                                   <  1:3373043/80‑1 (MQ=255)
                                                                             cACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCaa                                                                         >  2:5198298/1‑150 (MQ=255)
                                                                              aCCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAAc                                                                        >  1:694040/1‑150 (MQ=255)
                                                                                ccACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCg                                                                      >  2:4649101/1‑150 (MQ=255)
                                                                                ccACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCg                                                                      >  2:5446226/1‑150 (MQ=255)
                                                                                 cACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCgt                                                                     >  2:3018946/1‑150 (MQ=255)
                                                                                  acgacgTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCgtg                                                                    >  1:1080027/1‑150 (MQ=255)
                                                                                  acgacgTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCgtg                                                                    <  1:1472314/150‑1 (MQ=255)
                                                                                  acgacgTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCgtg                                                                    >  1:2497052/1‑150 (MQ=255)
                                                                                       gTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCATGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTg                                                               <  2:2485769/150‑1 (MQ=255)
                                                                                       gTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTg                                                               >  2:2163044/1‑150 (MQ=255)
                                                                                       gTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTg                                                               <  1:118051/150‑1 (MQ=255)
                                                                                        tCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGt                                                              >  1:5236553/1‑150 (MQ=255)
                                                                                        tCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGt                                                              >  2:3471251/1‑150 (MQ=255)
                                                                                            cGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCAc                                                          <  2:1897827/150‑1 (MQ=255)
                                                                                                gCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCg                                                      <  2:1705480/150‑1 (MQ=255)
                                                                                                  ggTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTTAACCGTGTTCTGTc                                                             >  2:860494/1‑141 (MQ=255)
                                                                                                  ggTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCggg                                                    <  2:2152287/150‑1 (MQ=255)
                                                                                                   gTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGGTCTGTCCACCCCGGGc                                                   <  2:4837497/150‑1 (MQ=255)
                                                                                                    ttGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCGGGCa                                                  <  1:2713589/150‑1 (MQ=255)
                                                                                                           ggCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGTCCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCGGGCACTGCGcc                                           >  2:3685974/1‑150 (MQ=255)
                                                                                                               gACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGACTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCGGGCACTGCGCCCGAg                                       >  2:3130174/1‑150 (MQ=255)
                                                                                                                aCCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCGGGCACTGCGCCCGAGc                                      <  2:2893620/150‑1 (MQ=255)
                                                                                                                 ccGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCGGGCACTGCGCCCGAGcc                                     <  2:3997323/150‑1 (MQ=255)
                                                                                                                  cGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCGGGCACTGCGCCCGAGCCg                                    >  2:5231704/1‑150 (MQ=255)
                                                                                                                      tGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCGGGCACTGCGCCCGAGCCGGATc                                <  1:361930/150‑1 (MQ=255)
                                                                                                                       gATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCGGGCACTGCGCCCGAGCCGGATcc                               <  2:4230375/150‑1 (MQ=255)
                                                                                                                          ccTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCGGGCACTGCGCCCGAGCCGGATCCTGt                            >  2:1499032/1‑150 (MQ=255)
                                                                                                                           cTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCGGGCACTGCGCCCGAGCCGGATCCTGTa                           <  2:67585/150‑1 (MQ=255)
                                                                                                                           cTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCGGGCACTGCGCCCGAGCCGGATCCTGTa                           <  1:3430410/150‑1 (MQ=255)
                                                                                                                           cTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCGGGCACTGCGCCCGAGCCGGATCCTGTa                           <  2:2794520/150‑1 (MQ=255)
                                                                                                                            tGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCGGGCACTGCGCCCGAGCCGGATCCTGTAg                          >  1:616628/1‑150 (MQ=255)
                                                                                                                            tGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCGGGCACTGCGCCCGAGCCGGATCCTGTAg                          >  2:4922311/1‑150 (MQ=255)
                                                                                                                             gATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCGGGCACTGCGCCCGAGCCGGATCCTGTAGc                         >  2:602942/1‑150 (MQ=255)
                                                                                                                                 gAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCGGGCACTGCGCCCGAGCCGGATCCTGTAGCCCct                     >  1:806462/1‑150 (MQ=255)
                                                                                                                                  aGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCCCCCCGGGCACTGCGCCCGAGCCGGATCCTGTAGCCCctc                    >  2:4180837/1‑150 (MQ=255)
                                                                                                                                   ggACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCGGGCACTGCGCCCGAGCCGGATCCTGTAGCCCctct                   <  1:724597/150‑1 (MQ=255)
                                                                                                                                    gACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCGGGCACTGCGCCCGAGCCGGATCCTGTAGCCCCTCTg                  >  2:4367341/1‑150 (MQ=255)
                                                                                                                                     aCGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCGGGCACTGCGCCCGAGCCGGATCCTGTAGCCCCTCTGc                 <  2:5062794/150‑1 (MQ=255)
                                                                                                                                       ggCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCGGGCACTGCGCCCGAGCCGGATCCTGTAGCCCCTCTGccc               >  1:867624/1‑150 (MQ=255)
                                                                                                                                       agtgAGGTCGAGCCCGACCTCGCTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCGGGCACTGCGCCCGAGCCGGATCCTGTAGCCCCTCTGccc               >  2:140507/4‑150 (MQ=255)
                                                                                                                                       agtgAGGTCGAGCCCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCGGGCACTGCGCCCGAGCCGGATCCTGTAGCCCCTCTGccc               >  1:4920451/4‑150 (MQ=255)
                                                                                                                                        gtgAGGTCGAGCCCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCGGGCACTGCGCCCGAGCCGGATCCTGTAGCCCCTCTGCCCa              >  2:4873319/3‑150 (MQ=255)
                                                                                                                                         cGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCGGGCACTGCGCCCGAGCCGGATCCTGTAGCCCCTCTGCCCAc             <  2:747904/150‑1 (MQ=255)
                                                                                                                                            ggTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGaa                                                                                 <  2:3284901/79‑1 (MQ=255)
                                                                                                                                            ggTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGaa                                                                                 >  1:3284901/1‑79 (MQ=255)
                                                                                                                                             gTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCGGGCACTGCGCCCGAGCCGGATCCTGTAGCCCCTCTGCCCACGCaa         >  2:1648228/1‑150 (MQ=255)
                                                                                                                                                   cTCGACCTCACTGTTGACCTGGCACGGCCCCCGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCGGGCACTGCGCCCGAGCCGGATCCTGTAGCCCCTCTGCCCACGCAACTgcgc   >  2:2692486/1‑150 (MQ=255)
                                                                                                                                                    tCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCGGGCACTGCGCCCGAGCCGGATCCTGTAGCCCCTCTGCCCACGCAACTGCGCt  >  2:4170615/1‑150 (MQ=255)
                                                                                                                                                    tCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCGGGCACTGCGCCCGAGCCGGATCCTGTAGCCCCTCTGCCCACGCAACTGCGCt  >  2:1827102/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GATGCGTTGACCCGTATCGAGATGCAGCAACTGATCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCGGGCACTGCGCCCGAGCCGGATCCTGTAGCCCCTCTGCCCACGCAACTGCGCT  >  NC_002947/295102‑295399

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