Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,485,365 T→A 6.0% L89Q (CTG→CAG)  puuA‑I → glutamate‑putrescine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,485,3650TA 6.0% 776.0 / 28.5 232L89Q (CTG→CAG) puuA‑Iglutamate‑putrescine ligase
Reads supporting (aligned to +/- strand):  ref base T (126/90);  new base A (2/12);  total (129/103)
Fisher's exact test for biased strand distribution p-value = 1.57e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 233 total aligned reads displayed.

CGCTACCCCATGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGCCAAGCGTGGCATTC  >  NC_002947/2485217‑2485514
                                                                                                                                                    |                                                                                                                                                     
cgcTACCCCATGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTg                                                                                                                                                      <  1:3429166/150‑1 (MQ=255)
cgcTACCCCATGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTg                                                                                                                                                      <  2:36941/150‑1 (MQ=255)
   tACCCCATGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTg                                                                                                                                                   >  2:1237531/1‑150 (MQ=255)
   tACCCCATGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTg                                                                                                                                                   >  2:408564/1‑150 (MQ=255)
     ccccATGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGcc                                                                                                                                                 >  1:3634548/1‑150 (MQ=255)
       ccATGGACCTGCTGGAAAAAGTCGCTGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCGGGTGCCgg                                                                                                                                               >  2:672826/1‑150 (MQ=255)
        cATGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGTCGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCggg                                                                                                                                              >  2:4061271/1‑150 (MQ=255)
        cATGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCggg                                                                                                                                              >  2:3227976/1‑150 (MQ=255)
        cATGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCggg                                                                                                                                              >  2:2289105/1‑150 (MQ=255)
          tGGACCTGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTACCCGGCACCAGGCGCCgggca                                                                                                                                            <  1:3312137/150‑3 (MQ=255)
          tGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCACCAGGCGCCgggca                                                                                                                                            <  1:1006080/150‑3 (MQ=255)
          tGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCACCAGGCGCCgggca                                                                                                                                            <  2:315819/150‑3 (MQ=255)
          tGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTACCCGGCACCAGGCGCCgggca                                                                                                                                            <  1:1754835/150‑3 (MQ=255)
                 gCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGaa                                                                                                                                     >  2:3068477/1‑150 (MQ=255)
                 gCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGaa                                                                                                                                     <  1:2282560/150‑1 (MQ=255)
                 gCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGaa                                                                                                                                     <  2:1110855/150‑1 (MQ=255)
                 gCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGaa                                                                                                                                     <  2:4521685/150‑1 (MQ=255)
                 gCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGaa                                                                                                                                     <  2:3286313/150‑1 (MQ=255)
                   tGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGc                                                                                                                                   >  2:796396/1‑150 (MQ=255)
                    ggAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGcc                                                                                                                                  <  1:527706/150‑1 (MQ=255)
                    ggAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGcc                                                                                                                                  >  2:3586207/1‑150 (MQ=255)
                     gAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCg                                                                                                                                 >  2:4123636/1‑150 (MQ=255)
                     gAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCg                                                                                                                                 >  1:3795501/1‑150 (MQ=255)
                     gAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCg                                                                                                                                 >  1:463315/1‑150 (MQ=255)
                        aaaGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTg                                                                                                                              <  2:4344279/150‑1 (MQ=255)
                             cGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTg                                                                                                                         <  1:189040/150‑1 (MQ=255)
                              gCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTgg                                                                                                                        >  1:1741466/1‑150 (MQ=255)
                              gCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTACCCGGCACCAGGCGcc                                                                                                                                                 >  2:2034835/1‑125 (MQ=255)
                              gCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTACCCGGCACCAGGCGcc                                                                                                                                                 <  1:2034835/125‑1 (MQ=255)
                                   cGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGaaaa                                                                                                                   <  1:1825602/150‑1 (MQ=255)
                                   cGGTGCGCCTCGGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCCGGGCGGCCGTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGaaaa                                                                                                                   <  1:1213948/150‑1 (MQ=255)
                                     gTGCGCCCCTGAAGCGGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAgg                                                                                                                 <  2:2890941/150‑1 (MQ=255)
                                       gcgcCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGaa                                                                                                               <  2:1323634/150‑1 (MQ=255)
                                       gcgcCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGaa                                                                                                               <  2:1436546/150‑1 (MQ=255)
                                       gcgcCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGaa                                                                                                               <  2:4285331/150‑1 (MQ=255)
                                         gcCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGt                                                                                                             <  1:4167795/150‑1 (MQ=255)
                                         gcCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGt                                                                                                             >  1:2892840/1‑150 (MQ=255)
                                          ccTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTg                                                                                                            >  1:2265051/1‑150 (MQ=255)
                                           ctctGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGa                                                                                                           >  1:3529297/1‑150 (MQ=255)
                                           ctctGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGa                                                                                                           >  2:2776797/1‑150 (MQ=255)
                                            tctGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGAt                                                                                                          >  1:1217316/1‑150 (MQ=255)
                                                aaGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGcc                                                                                                      >  1:4379491/1‑150 (MQ=255)
                                                 aGCTGCCGCAGAACTGCGTGCGGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTCTTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGccc                                                                                                     <  2:3736567/150‑1 (MQ=255)
                                                    tGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGa                                                                                                  >  2:2353175/1‑150 (MQ=255)
                                                         cAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTg                                                                                             <  2:758271/150‑1 (MQ=255)
                                                                gCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGt                                                                                      >  1:2592211/1‑150 (MQ=255)
                                                                 cGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTc                                                                                     <  1:218019/150‑1 (MQ=255)
                                                                 cGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGGGACTAGTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTc                                                                                     <  2:4250343/150‑1 (MQ=255)
                                                                   tgctgcTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCg                                                                                   >  1:1452654/1‑150 (MQ=255)
                                                                   tgctgcTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCg                                                                                   >  1:1272629/1‑150 (MQ=255)
                                                                        ctgGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGc                                                                              <  2:1152212/150‑1 (MQ=255)
                                                                         tgGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCa                                                                             <  1:3512162/150‑1 (MQ=255)
                                                                           ggTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACt                                                                           >  2:4181455/1‑150 (MQ=255)
                                                                           ggTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACt                                                                           >  1:3042128/1‑150 (MQ=255)
                                                                           ggTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACt                                                                           <  2:1261747/150‑1 (MQ=255)
                                                                           ggTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACt                                                                           >  2:3324680/1‑150 (MQ=255)
                                                                           ggTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCAGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACt                                                                           >  2:3803589/1‑150 (MQ=255)
                                                                            gTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGGCGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTg                                                                          >  2:1726103/1‑150 (MQ=255)
                                                                                ccAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCaaa                                                                      >  2:3742966/1‑150 (MQ=255)
                                                                                ccAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCaaa                                                                      <  1:3974717/150‑1 (MQ=255)
                                                                                ccAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCaaa                                                                      <  1:4158276/150‑1 (MQ=255)
                                                                                ccAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCaaa                                                                      <  2:525978/150‑1 (MQ=255)
                                                                                ccAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCaaa                                                                      >  1:2950365/1‑150 (MQ=255)
                                                                                ccAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCaaa                                                                      >  1:680664/1‑150 (MQ=255)
                                                                                ccAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCaaa                                                                      >  1:825127/1‑150 (MQ=255)
                                                                                     gcggcCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGa                                                                 >  2:234302/1‑150 (MQ=255)
                                                                                           tttACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAAt                                                           >  2:4512084/1‑150 (MQ=255)
                                                                                           tttACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAAt                                                           >  2:351455/1‑150 (MQ=255)
                                                                                              aCCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCgg                                                                                                                                               <  1:4476059/63‑1 (MQ=255)
                                                                                              aCCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCg                                                        >  2:4071980/1‑150 (MQ=255)
                                                                                                ccGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAg                                                      >  2:219829/1‑150 (MQ=255)
                                                                                                ccGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAg                                                      >  1:3373014/1‑150 (MQ=255)
                                                                                                ccGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAg                                                      >  2:1565245/1‑150 (MQ=255)
                                                                                                      ggTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGt                                                >  1:2612861/1‑150 (MQ=255)
                                                                                                       gTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTg                                               <  1:3227976/150‑1 (MQ=255)
                                                                                                            tATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTg                                          <  2:2577501/150‑1 (MQ=255)
                                                                                                            tATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTg                                          >  1:1235074/1‑150 (MQ=255)
                                                                                                               tGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTg                                       <  1:2353175/150‑1 (MQ=255)
                                                                                                                gCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTTCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGTAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTgg                                      <  1:646101/150‑1 (MQ=255)
                                                                                                                 cTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGgc                                     >  2:717717/1‑150 (MQ=255)
                                                                                                                 cTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGgc                                     >  1:1178538/1‑150 (MQ=255)
                                                                                                                 cTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGgc                                     >  1:1640163/1‑150 (MQ=255)
                                                                                                                  ttCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGgcg                                    >  1:2810385/1‑150 (MQ=255)
                                                                                                                   tCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGgcgc                                   <  2:409874/150‑1 (MQ=255)
                                                                                                                        gACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGtt                              >  2:2698029/1‑150 (MQ=255)
                                                                                                                           ggCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGttctt                           >  2:89411/1‑150 (MQ=255)
                                                                                                                           ggCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGttctt                           <  1:4417183/150‑1 (MQ=255)
                                                                                                                            gCGACCCGGAGGCTGCCCGGCGCCGGGGGCCGGGTACGCTGCAGCCGGTGCGCTGGGACAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCGGGCGCGTGttcttc                          <  2:1957692/150‑1 (MQ=255)
                                                                                                                                  cGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTc                    >  1:1612943/1‑150 (MQ=255)
                                                                                                                                   ggATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCt                   <  1:4253433/150‑1 (MQ=255)
                                                                                                                                   ggATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCt                   <  1:4061280/150‑1 (MQ=255)
                                                                                                                                   ggATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCt                   >  1:2047137/1‑150 (MQ=255)
                                                                                                                                      tgctgcCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGc                <  1:4123636/150‑1 (MQ=255)
                                                                                                                                        gtacccGGCACCAGGCGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGCCa              >  2:4479311/4‑150 (MQ=255)
                                                                                                                                         tgcCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGCCaa             >  1:2063378/1‑150 (MQ=255)
                                                                                                                                          gcCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGCCAAg            <  1:3100344/150‑1 (MQ=255)
                                                                                                                                           cccGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGCCAAGc           >  2:3402971/1‑150 (MQ=255)
                                                                                                                                             cGGCGCCTGGTGCCGGGTACGCTGACGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGCCAAGCGt         <  2:2466200/150‑1 (MQ=255)
                                                                                                                                             cGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGCCAAGCGt         <  2:1402583/150‑1 (MQ=255)
                                                                                                                                                    tGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGCCAAGCGTGGCATTc  >  1:700030/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGCTACCCCATGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGCCAAGCGTGGCATTC  >  NC_002947/2485217‑2485514

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: