Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,513,023 T→A 6.1% Q48L (CAG→CTG)  PP_2201 ← alpha/beta family hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,513,0230TA 6.1% 852.8 / 29.8 255Q48L (CAG→CTG) PP_2201alpha/beta family hydrolase
Reads supporting (aligned to +/- strand):  ref base T (115/124);  new base A (5/11);  total (120/135)
Fisher's exact test for biased strand distribution p-value = 2.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 259 total aligned reads displayed.

GCGCGATGACCCGGTAGCCGGCCTTGCTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAGGGCGGTTGCGCGCAGCATGGGGGCA  >  NC_002947/2512874‑2513172
                                                                                                                                                     |                                                                                                                                                     
gcgcGATGACCCGGTAGCCGGCCTTGCTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCt                                                                                                                                                       >  1:1912680/1‑150 (MQ=255)
gcgcGATGACCCGGTAGCCGGCCTTGCTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGAATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCt                                                                                                                                                       >  1:5467253/1‑150 (MQ=255)
       gACCCGGGAGCCGGCCTTGCTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTg                                                                                                                                                <  2:1536180/150‑1 (MQ=255)
        aCCCGGTAGCCGGCCTTGCTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTgg                                                                                                                                               <  1:837316/150‑1 (MQ=255)
        aCCCGGTAGCCGGCCTTGCTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTgg                                                                                                                                               >  1:3611262/1‑150 (MQ=255)
         cccGGTAGCCGGCCTTGCTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCAGGCCCTgca                                                                                                                                              <  2:5379431/150‑3 (MQ=255)
          ccGGTAGCCGGCCTTGCTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGa                                                                                                                                             <  2:2923660/150‑1 (MQ=255)
          ccGGTAGCCGGCCTTGCTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCAGGCCCTgcag                                                                                                                                             <  2:2116811/150‑4 (MQ=255)
           cGGTAGCCGGCCTTGCTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGAc                                                                                                                                            <  2:610074/150‑1 (MQ=255)
           cGGTAGCCGGCCTTGCTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCAGGCCCTgcaga                                                                                                                                            <  1:112548/150‑5 (MQ=255)
              tAGCCGGCCTTGCTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTg                                                                                                                                         <  1:1281226/150‑1 (MQ=255)
                gCCGGCCTTGCTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCCGCAGGGCCTGGCCCTGGGACTTGaa                                                                                                                                       <  2:2558897/150‑1 (MQ=255)
                   ggCCTTGCTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTc                                                                                                                                    <  1:3905736/150‑1 (MQ=255)
                     ccTTGCTCAGGGCGTCGTTGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGa                                                                                                                                  <  1:2931139/150‑1 (MQ=255)
                     ccTTGCTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGa                                                                                                                                  >  2:4735634/1‑150 (MQ=255)
                           tCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCt                                                                                                                            >  2:4266136/1‑150 (MQ=255)
                             aGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTg                                                                                                                          <  2:3997706/150‑1 (MQ=255)
                              gggCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGa                                                                                                                         <  2:5497180/150‑1 (MQ=255)
                              gggCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGa                                                                                                                         <  1:1019245/150‑1 (MQ=255)
                                  gTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAgggg                                                                                                                     <  2:15171/150‑1 (MQ=255)
                                  gTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCATGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCGTGAgggg                                                                                                                     >  1:4554253/1‑150 (MQ=255)
                                   tCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAggggg                                                                                                                    >  1:4077295/1‑150 (MQ=255)
                                    cGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGt                                                                                                                   <  2:1087745/150‑1 (MQ=255)
                                     gATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTg                                                                                                                  >  1:121863/1‑150 (MQ=255)
                                     gATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTg                                                                                                                  >  1:5048402/1‑150 (MQ=255)
                                     gATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTg                                                                                                                  >  1:1988831/1‑150 (MQ=255)
                                        ggtggtTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGTGGGTGCgg                                                                                                               <  1:1695514/150‑1 (MQ=255)
                                         gtggtTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCggg                                                                                                              >  2:4302659/1‑150 (MQ=255)
                                            gtTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAg                                                                                                           <  2:2291108/150‑1 (MQ=255)
                                             tttCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGc                                                                                                          >  1:3472507/1‑150 (MQ=255)
                                              ttCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTGGAGGGGGTGCGGGTAGCt                                                                                                         >  2:221549/1‑150 (MQ=255)
                                                  cAGGCGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTGGAAGTCGAAATGCTTGAGGGGGTGCGGGTCGCTGAAg                                                                                                     >  2:5329557/1‑150 (MQ=255)
                                                       ggcggcGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTc                                                                                                <  1:2377912/150‑1 (MQ=255)
                                                       ggcggcGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTc                                                                                                <  2:1301682/150‑1 (MQ=255)
                                                       ggcggcGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTc                                                                                                <  1:2630963/150‑1 (MQ=255)
                                                       ggcggcGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTc                                                                                                <  2:199673/150‑1 (MQ=255)
                                                        gcggcgCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCg                                                                                               >  2:1523168/1‑150 (MQ=255)
                                                              cAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTgg                                                                                         <  2:1605696/150‑1 (MQ=255)
                                                              cAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTgg                                                                                         <  2:3376018/150‑1 (MQ=255)
                                                              cAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTgg                                                                                         <  2:1017700/150‑1 (MQ=255)
                                                              cAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTgg                                                                                         <  2:1370517/150‑1 (MQ=255)
                                                              cAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTgg                                                                                         <  2:919019/150‑1 (MQ=255)
                                                                gAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCAGGCCCTgca                                                                                                                                              >  1:2898334/1‑93 (MQ=255)
                                                                 aaGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTggggg                                                                                      <  2:1912680/150‑1 (MQ=255)
                                                                   gTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCGCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGcc                                                                                    >  2:93511/1‑150 (MQ=255)
                                                                   gTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGcc                                                                                    >  1:372782/1‑150 (MQ=255)
                                                                      cTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCgta                                                                                 >  1:1633740/1‑150 (MQ=255)
                                                                      attGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGTCACGTCAATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCTAGCTGGGGGCCgta                                                                                 <  2:3736822/149‑1 (MQ=255)
                                                                        tGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCgtagg                                                                               <  1:745220/150‑1 (MQ=255)
                                                                           cGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCgtaggtag                                                                            >  1:523253/1‑150 (MQ=255)
                                                                             tGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGc                                                                          >  1:3863509/1‑150 (MQ=255)
                                                                               cATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGa                                                                        >  2:2967908/1‑150 (MQ=255)
                                                                                  cAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTc                                                                     <  1:2168545/150‑1 (MQ=255)
                                                                                  cAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTc                                                                     >  2:1779392/1‑150 (MQ=255)
                                                                                  cAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTc                                                                     <  1:498285/150‑1 (MQ=255)
                                                                                  cAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTc                                                                     <  1:4532076/150‑1 (MQ=255)
                                                                                   aGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCg                                                                    >  1:1060707/1‑150 (MQ=255)
                                                                                     caccacGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGc                                                                  >  2:4981397/1‑150 (MQ=255)
                                                                                      accacGTTGGGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCa                                                                 <  2:2894408/150‑1 (MQ=255)
                                                                                       ccacGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAg                                                                >  1:1591221/1‑150 (MQ=255)
                                                                                        nacGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGc                                                               >  1:630299/2‑150 (MQ=255)
                                                                                        cacGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGc                                                               >  2:2420471/1‑150 (MQ=255)
                                                                                        cacGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGACAGc                                                               >  2:330452/1‑150 (MQ=255)
                                                                                            ttGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATg                                                           <  1:4325023/150‑1 (MQ=255)
                                                                                              gtgGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGc                                                         >  2:1726066/1‑150 (MQ=255)
                                                                                                ggCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCg                                                       >  1:1447369/1‑150 (MQ=255)
                                                                                                     ttGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAActg                                                  >  1:3109034/1‑150 (MQ=255)
                                                                                                       ggcctggcctTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAActgct                                                <  1:5184662/150‑1 (MQ=255)
                                                                                                        gcctggcctTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAActgctg                                               >  2:5272789/1‑150 (MQ=255)
                                                                                                        gcctggcctTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAActgctg                                               <  2:1447369/150‑1 (MQ=255)
                                                                                                         cctggcctTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTgg                                              >  1:4693533/1‑150 (MQ=255)
                                                                                                          ctggcctTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGc                                             >  2:115671/1‑150 (MQ=255)
                                                                                                           tggcctTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCa                                            >  1:2827561/1‑150 (MQ=255)
                                                                                                           tggcctTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCa                                            >  1:2827603/1‑150 (MQ=255)
                                                                                                             gcctTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAag                                          <  2:273309/150‑1 (MQ=255)
                                                                                                             gcctTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAag                                          <  1:4594370/150‑1 (MQ=255)
                                                                                                               ctTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAaggc                                        >  1:5360238/1‑150 (MQ=255)
                                                                                                               ctTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAaggc                                        <  2:4543752/150‑1 (MQ=255)
                                                                                                                 tGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAaggcag                                      <  2:471015/150‑1 (MQ=255)
                                                                                                                     cgggcACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCa                                  >  1:3164292/1‑150 (MQ=255)
                                                                                                                     cgggcACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCa                                  <  1:2667802/150‑1 (MQ=255)
                                                                                                                            gTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAg                           >  1:4255121/1‑150 (MQ=255)
                                                                                                                               cATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAGGGc                        >  1:5617271/1‑150 (MQ=255)
                                                                                                                                aTGTAGCCCATCTGCAGGGCCTGGCCCTGGGCCTTGAAGTCGATATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGAGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAGGGCg                       >  2:5135037/1‑150 (MQ=255)
                                                                                                                                       ccATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGa                                                                                                       >  2:687844/1‑63 (MQ=255)
                                                                                                                                       ccATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAGGGCGGTTgcgc                >  1:69737/1‑150 (MQ=255)
                                                                                                                                       ccATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAGGGCGGTTgcgc                <  2:3407064/150‑1 (MQ=255)
                                                                                                                                       ccATCTGCAGGGCCTGGCCCTGGGACTGGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCCAGGGCGGTTgcgc                >  2:1210864/1‑150 (MQ=255)
                                                                                                                                       acATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTCGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAGGGCGGTTgcgc                >  2:5456802/2‑150 (MQ=255)
                                                                                                                                        cATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAGGGCGGTTgcgcg               >  1:4873831/1‑150 (MQ=255)
                                                                                                                                        cATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAGGGCGGTTgcgcg               >  2:3249050/1‑150 (MQ=255)
                                                                                                                                         aTCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAGGGCGGTTgcgcgc              >  1:2781896/1‑150 (MQ=255)
                                                                                                                                           tcccAGGGCCAGGCCCTGGGACTGGAAGTCGAAATGCTGGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCACGGCAGGCCAGGGCAAGGGCGGGTGCGCgcag            >  2:462937/4‑150 (MQ=255)
                                                                                                                                            tGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAGGGCGGTTGCGCgcagc           >  2:726154/1‑150 (MQ=255)
                                                                                                                                            tGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAGGGCGGTTGCGCgcagc           >  1:25403/1‑150 (MQ=255)
                                                                                                                                            cccAGGGCCAGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAGGGCGGTTGCGCgca             >  2:5009971/3‑148 (MQ=255)
                                                                                                                                              cAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGTCAGGGCAAGGGCGGTTGCGCGCAGCAt         >  1:4641260/1‑150 (MQ=255)
                                                                                                                                                gggCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAGGGCGGTTGCGCGCAGCATgg       <  2:4663535/150‑1 (MQ=255)
                                                                                                                                                gggCCTGGCCCGGGGACTTGAAGTCGAACTGCTTGAGGGGGTGCGGGTAGCTGAAGCCTGCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAGGGCGGTTGCGCGCAGCATgg       <  2:3855010/150‑1 (MQ=255)
                                                                                                                                                   ccTGGCCCTGGGACTTGAAGTCGAAATGCGTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAGGGCGGTTGCGCGCAGCATggggg    >  2:127183/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCGCGATGACCCGGTAGCCGGCCTTGCTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAGGGCGGTTGCGCGCAGCATGGGGGCA  >  NC_002947/2512874‑2513172

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: