Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 713,100 T→A 5.6% I29N (ATC→AAC) ‡ lspA → lipoprotein signal peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947713,1000TA 5.6% 1031.7 / 33.4 303I29N (ATC→AAC) ‡lspAlipoprotein signal peptidase
Reads supporting (aligned to +/- strand):  ref base T (134/152);  new base A (8/9);  total (142/161)
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.34e-01

Only 100 of 307 total aligned reads displayed.

TCGGCAGCCATCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCGGGTTGGCAGCGCTGGT  >  NC_002947/712952‑713249
                                                                                                                                                    |                                                                                                                                                      
   gCAGCCATCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGAc                                                                                                                                                    <  2:6436978/150‑1 (MQ=255)
   gCAGCCATCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGAc                                                                                                                                                    <  1:2652531/150‑1 (MQ=255)
     aGCCATCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCa                                                                                                                                                  <  1:3877763/150‑1 (MQ=255)
     aGCCATCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCa                                                                                                                                                  <  2:403128/150‑1 (MQ=255)
        cATCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGGTGGCCTGGTCGATGACCAggac                                                                                                                                              <  2:4089684/150‑3 (MQ=255)
        cATCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCGATGACCAgga                                                                                                                                               <  2:3113811/150‑2 (MQ=255)
          tCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGcca                                                                                                                                            >  1:4983549/1‑150 (MQ=255)
           ccGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGccac                                                                                                                                           <  2:2738294/150‑1 (MQ=255)
             ggAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCCGTGGCTGAGCTTGCTGGTCCTGGTCGATGACCAggaccagc                                                                                                                                          <  1:4183183/150‑7 (MQ=255)
                gCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGc                                                                                                                                      <  1:2873926/150‑1 (MQ=255)
                gCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGc                                                                                                                                      <  1:348430/150‑1 (MQ=255)
                    ccGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTa                                                                                                                                  >  2:4679080/1‑150 (MQ=255)
                     cGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTAc                                                                                                                                 >  1:6490801/1‑150 (MQ=255)
                     cGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTAc                                                                                                                                 >  2:2566473/1‑150 (MQ=255)
                     cGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTAc                                                                                                                                 >  2:3728341/1‑150 (MQ=255)
                     cGGAAATCTGCGGCCGTTGCGTCGACAACCGGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTAc                                                                                                                                 <  1:5625647/150‑1 (MQ=255)
                          aTCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTcaa                                                                                                                            <  2:6248017/150‑1 (MQ=255)
                                 gCCGTTGCGTCGACAACCTGACCGGCTCGGGAGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCAGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGcc                                                                                                                     <  2:5835515/150‑1 (MQ=255)
                                   cGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCt                                                                                                                   >  1:1577919/1‑150 (MQ=255)
                                     ttGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGa                                                                                                                 <  1:5491453/150‑1 (MQ=255)
                                     ttGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGa                                                                                                                 <  1:6716500/150‑1 (MQ=255)
                                       gCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGAcc                                                                                                               >  2:6960092/1‑150 (MQ=255)
                                         gTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCAt                                                                                                             >  1:4044466/1‑150 (MQ=255)
                                          tCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATg                                                                                                            >  1:5329917/1‑150 (MQ=255)
                                           cGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGt                                                                                                           <  2:138076/150‑1 (MQ=255)
                                            gACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTa                                                                                                          <  2:3407221/150‑1 (MQ=255)
                                             aCAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTAc                                                                                                         <  1:2155757/150‑1 (MQ=255)
                                               aaCCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACca                                                                                                       >  2:451539/1‑150 (MQ=255)
                                               aaCCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACca                                                                                                       <  2:6337841/150‑1 (MQ=255)
                                                aCCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACcag                                                                                                      >  1:3572250/1‑150 (MQ=255)
                                                 ccTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACcagc                                                                                                     >  2:893727/1‑150 (MQ=255)
                                                   tGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCAAGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACcagcag                                                                                                   >  1:4587697/1‑150 (MQ=255)
                                                       cGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGGTGGCCTGGTCGATGACCAGGACCAGCAAGCtca                                                                                                                                   >  2:3461808/1‑112 (MQ=25)
                                                       cGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCg                                                                                               >  1:2338974/1‑150 (MQ=255)
                                                         gCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTg                                                                                             <  2:1281391/150‑1 (MQ=255)
                                                         gCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTg                                                                                             <  1:1271522/150‑1 (MQ=255)
                                                         gCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTg                                                                                             <  1:5409726/150‑1 (MQ=255)
                                                         gCTAGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTg                                                                                             <  1:6149072/150‑1 (MQ=255)
                                                          cTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTgg                                                                                            <  2:388949/150‑1 (MQ=255)
                                                             gggTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCa                                                                                         >  1:6381631/1‑150 (MQ=255)
                                                                tGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTc                                                                                      <  1:6317220/150‑1 (MQ=255)
                                                                    gTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTg                                                                                                                  <  2:4628173/118‑1 (MQ=255)
                                                                    gTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGa                                                                                  <  1:5186694/150‑1 (MQ=255)
                                                                    gTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGa                                                                                  >  1:3345695/1‑150 (MQ=255)
                                                                    gTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGa                                                                                  >  1:3287452/1‑150 (MQ=255)
                                                                     tGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTgg                                                                                            <  2:6901113/139‑1 (MQ=255)
                                                                       cgcCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTTCCAGCAGATCGTGGTCATTCCTGacta                                                                               <  2:3580843/150‑1 (MQ=255)
                                                                       cgcCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGacta                                                                               <  2:4156135/150‑1 (MQ=255)
                                                                          cACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGactacct                                                                            >  2:5362629/1‑148 (MQ=255)
                                                                          cACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACtt                                                                            <  2:5053328/150‑1 (MQ=255)
                                                                          cACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACtt                                                                            >  2:1668044/1‑150 (MQ=255)
                                                                           aCTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGactactac                                                                           >  2:523216/1‑148 (MQ=255)
                                                                                 ccTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTgg                                                                     >  2:7189185/1‑150 (MQ=255)
                                                                                   tAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGAc                                                                   >  1:2106440/1‑150 (MQ=255)
                                                                                       ccTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTg                                                               <  1:4445399/150‑1 (MQ=255)
                                                                                        cagCAGCGGGGCGCTTCGGGCGCCTAGCATGGCTGTGGCAGAGCTTGCTGGTCCTGCTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATCTACCACCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTgg                                                              <  2:6772538/148‑1 (MQ=255)
                                                                                         tGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGc                                                             <  1:6213846/150‑1 (MQ=255)
                                                                                         tGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGc                                                             >  1:7242649/1‑150 (MQ=255)
                                                                                          gCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCt                                                            <  1:3849983/150‑1 (MQ=255)
                                                                                          gCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCt                                                            <  1:4352684/150‑1 (MQ=255)
                                                                                            aGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTa                                                          <  2:6693591/150‑1 (MQ=255)
                                                                                             gCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTac                                                         <  1:7395081/150‑1 (MQ=255)
                                                                                               ggggCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTacaa                                                       <  2:1201084/150‑1 (MQ=255)
                                                                                               ggggCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTacaa                                                       <  1:5800031/150‑1 (MQ=255)
                                                                                                gggCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTacaac                                                      >  1:4669513/1‑150 (MQ=255)
                                                                                                         gggCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTTcc                                             <  2:2016794/150‑1 (MQ=255)
                                                                                                          ggCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCcg                                            >  1:1896943/1‑150 (MQ=255)
                                                                                                              ccTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGtt                                        >  2:5493288/1‑150 (MQ=255)
                                                                                                               cTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTc                                       >  1:705780/1‑150 (MQ=255)
                                                                                                               cTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTc                                       >  2:4769601/1‑150 (MQ=255)
                                                                                                                ttGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCa                                      <  2:6061734/150‑1 (MQ=255)
                                                                                                                ttGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCa                                      <  2:2934871/150‑1 (MQ=255)
                                                                                                                     tggctgtggctgAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTc                                 <  2:1167351/150‑1 (MQ=255)
                                                                                                                     tggctgtggctgAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTc                                 >  2:1944999/1‑150 (MQ=255)
                                                                                                                         tgtggctgAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTgg                             <  2:4396402/150‑1 (MQ=255)
                                                                                                                         tgtggctgAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTgg                             <  1:4893608/150‑1 (MQ=255)
                                                                                                                         tgtggctgAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTgg                             <  2:7127705/150‑1 (MQ=255)
                                                                                                                           tggctgAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCt                           <  1:1175271/150‑1 (MQ=255)
                                                                                                                             gctgAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGa                         <  2:1689714/150‑1 (MQ=255)
                                                                                                                                    ttGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCg                  >  2:1849394/1‑150 (MQ=255)
                                                                                                                                     ttggtggcctggtcGATGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCCGTGCCGCGTTCAGCTTCCTGGCGGACGGAGCgg                 >  2:1091717/8‑150 (MQ=255)
                                                                                                                                     tGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCgg                 <  2:82037/150‑1 (MQ=255)
                                                                                                                                     tGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCgg                 <  2:5154670/150‑1 (MQ=255)
                                                                                                                                      gCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCggg                >  1:5958334/1‑150 (MQ=255)
                                                                                                                                        tggtcctggtcATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTaca                                                        <  1:3031310/108‑1 (MQ=255)
                                                                                                                                        tggtcctggtcATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTaca                                                        >  2:3031310/1‑108 (MQ=255)
                                                                                                                                        tggtcctggtcATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCGGGtt              <  2:3813709/150‑1 (MQ=255)
                                                                                                                                         tggcctggtcGATGACCAGGCCACCAAGCTGTACTTCGACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCGGGTTg             >  1:1377074/4‑150 (MQ=255)
                                                                                                                                          gtcctggtcATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCt                               <  2:2374728/131‑1 (MQ=255)
                                                                                                                                          gtcctggtcATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCGGGTTgg            >  2:1225454/1‑150 (MQ=255)
                                                                                                                                           gcctggtcGATGACCAGGACCAGCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCGGGTTgg            >  2:1528916/2‑150 (MQ=255)
                                                                                                                                            cctggtcATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCGGGTTGGCa          >  1:2401579/1‑150 (MQ=255)
                                                                                                                                             ctggtcATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCGGGTTGGCAg         >  1:4718958/1‑150 (MQ=255)
                                                                                                                                              ttgtcATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCGGGTTGGCAgc        >  2:2289592/3‑150 (MQ=255)
                                                                                                                                              tggtcATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCGGGTTGGCAgc        >  2:733375/1‑150 (MQ=255)
                                                                                                                                                 tcATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCGGGTTGGCAGCGCt     >  1:4857819/1‑150 (MQ=255)
                                                                                                                                                  cATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTGCAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCGGGTTGGCAGCGCTg    <  2:2821150/150‑1 (MQ=255)
                                                                                                                                                   aTCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCGGGGTGGCAGCGCTgg   >  1:4802633/1‑150 (MQ=255)
                                                                                                                                                    tCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCGGGTTGGCAGCGCTGgt  <  2:4983549/150‑1 (MQ=255)
                                                                                                                                                    tCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCGGGTTGGCAGCGCTGgt  >  1:3094686/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
TCGGCAGCCATCCGGAGCACCCGGAAATCTGCGGCCGTTGCGTCGACAACCTGACCGGCTCGGGTGAGGTGCGCCACTATGCCTAACCCTGCAGCGGGGCGCTTCGGGCGCCTTGCATGGCTGTGGCTGAGCTTGCTGGTCCTGGTCATCGACCAGGCCACCAAGCTGTACTTCAACAACGCCCTGACCATGTACCAGCAGATCGTGGTCATTCCTGACTACTTCAGCTGGACCCTGGCTTACAACACCGGTGCCGCGTTCAGCTTCCTGGCTGACGGAGCGGGTTGGCAGCGCTGGT  >  NC_002947/712952‑713249

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: