Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,556,425 G→C 6.9% E26Q (GAG→CAG)  PP_3139 → group 1 family glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,556,4250GC 6.9% 567.9 / 24.5 173E26Q (GAG→CAG) PP_3139group 1 family glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base G (86/75);  new base C (6/6);  total (92/81)
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.33e-01

Only 100 of 174 total aligned reads displayed.

GCCCGTAGTCGACAGCCTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGGTTGTTGGCACCGACGTGGCAGGTA  >  NC_002947/3556276‑3556571
                                                                                                                                                     |                                                                                                                                                  
gCCCGTAGTCGACAGCCTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCAGCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGACAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTgg                                                                                                                                                    >  2:1202113/1‑150 (MQ=255)
   cGTAGTCGACAGCCTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAgc                                                                                                                                                 <  2:3525141/150‑1 (MQ=255)
       gTCGACAGCCTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCa                                                                                                                                             >  1:1860767/1‑150 (MQ=255)
          gACAGCCTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAgg                                                                                                                                          >  1:4614111/1‑150 (MQ=255)
              gCCTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGtt                                                                                                                                      <  2:3316047/150‑1 (MQ=255)
              gCCTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGtt                                                                                                                                      <  1:5309344/150‑1 (MQ=255)
              gCCTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGtt                                                                                                                                      <  1:4012258/150‑1 (MQ=255)
              gCCTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTCCAGCGCCagtatttc                                                                                                                                      <  1:2254724/150‑8 (MQ=255)
               ccTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTg                                                                                                                                     >  2:3711179/1‑150 (MQ=255)
               ccTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTg                                                                                                                                     >  2:3554861/1‑150 (MQ=255)
                 tgaagtgaagGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAg                                                                                                                                   >  2:666300/1‑150 (MQ=255)
                 tgaagtgaagGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAg                                                                                                                                   >  1:5179457/1‑150 (MQ=255)
                 tgaagtgaagGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAg                                                                                                                                   <  1:3711179/150‑1 (MQ=255)
                   aagtgaagGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCGGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGt                                                                                                                                 <  1:4482113/150‑1 (MQ=255)
                       gaagGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCaa                                                                                                                             <  1:381159/150‑1 (MQ=255)
                           gTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGta                                                                                                                         >  2:696888/1‑150 (MQ=255)
                           gTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGta                                                                                                                         >  2:3894795/1‑150 (MQ=255)
                              gCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCAGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAg                                                                                                                      <  2:2105335/150‑1 (MQ=255)
                                 aGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGccc                                                                                                                   >  1:484087/1‑150 (MQ=255)
                                    gAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGc                                                                                                                >  2:660906/1‑150 (MQ=255)
                                       ggTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCgtg                                                                                                             <  2:2912083/150‑1 (MQ=255)
                                       ggTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTCCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCgtg                                                                                                             <  2:3799316/150‑1 (MQ=255)
                                         ttGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCgtggg                                                                                                           <  2:4614111/150‑1 (MQ=255)
                                          tGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCgtgggt                                                                                                          >  2:22249/1‑150 (MQ=255)
                                          tGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATAGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTCGCCCTGGAGCGCCAAGCGGTTGAGGTGCAACGTATAGCCCTGCgagggt                                                                                                          >  1:876498/1‑150 (MQ=255)
                                           ggTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCgtgggtg                                                                                                         >  2:2467010/1‑150 (MQ=255)
                                           ggTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCgtgggtg                                                                                                         >  2:5194617/1‑150 (MQ=255)
                                            gTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCgtgggtgg                                                                                                        <  2:757907/150‑1 (MQ=255)
                                                  aaGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGcc                                                                                                  >  2:3097069/1‑150 (MQ=255)
                                                     cGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCACTCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGaa                                                                                               <  2:2698914/150‑1 (MQ=255)
                                                     cGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGaa                                                                                               <  1:3124603/150‑1 (MQ=255)
                                                      gggTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAg                                                                                              <  2:1860767/150‑1 (MQ=255)
                                                        gTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTCCAGCGCCa                                                                                                                                             <  2:867416/101‑1 (MQ=255)
                                                        gTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTCCAGCGCCa                                                                                                                                             >  1:867456/1‑101 (MQ=255)
                                                         tttGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGc                                                                                           >  2:4636287/1‑150 (MQ=255)
                                                         tttGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGc                                                                                           >  2:4203127/1‑150 (MQ=255)
                                                         tttGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGc                                                                                           >  2:1456468/1‑150 (MQ=255)
                                                              aTCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGAt                                                                                      >  2:4971963/1‑150 (MQ=255)
                                                              aTCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGAt                                                                                      >  1:4246131/1‑150 (MQ=255)
                                                               tCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATg                                                                                     >  1:727143/1‑150 (MQ=255)
                                                               tCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATg                                                                                     >  2:2049886/1‑150 (MQ=255)
                                                                 tGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGcc                                                                                   >  2:4320151/1‑150 (MQ=255)
                                                                   gATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGa                                                                                 <  1:110970/150‑1 (MQ=255)
                                                                     ttCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAgg                                                                               >  2:2348145/1‑150 (MQ=255)
                                                                        gTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGGAGGATGCCGAGGACg                                                                            >  1:911294/1‑150 (MQ=255)
                                                                           tGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCaa                                                                         >  1:1666324/1‑150 (MQ=255)
                                                                           tGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCaa                                                                         <  1:2851226/150‑1 (MQ=255)
                                                                             cgcATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGc                                                                       <  1:391851/150‑1 (MQ=255)
                                                                                 ttGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAAc                                                                   <  2:1763522/150‑1 (MQ=255)
                                                                                   gCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACgg                                                                 <  1:4113184/150‑1 (MQ=255)
                                                                                    cTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACggg                                                                <  1:3357745/150‑1 (MQ=255)
                                                                                     ttACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACgggg                                                               >  2:4517360/1‑150 (MQ=255)
                                                                                     ttACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACgggg                                                               <  1:2481569/150‑1 (MQ=255)
                                                                                       agttCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGTCGCAAGCGAACGGGGCa                                                             >  2:3287849/3‑150 (MQ=255)
                                                                                       aCTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCa                                                             >  1:4203797/1‑150 (MQ=255)
                                                                                         ttCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAg                                                           >  1:1046272/1‑150 (MQ=255)
                                                                                         ttCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAg                                                           >  2:4371532/1‑150 (MQ=255)
                                                                                         ttCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAg                                                           >  1:1338517/1‑150 (MQ=255)
                                                                                          tcatcaATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGa                                                          <  1:1940958/150‑1 (MQ=255)
                                                                                             tcaatcaGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGAccc                                                       <  1:1057090/150‑1 (MQ=255)
                                                                                             tcaatcaGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGAccc                                                       >  1:4376434/1‑150 (MQ=255)
                                                                                                 tcaGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTa                                                   <  2:136637/150‑1 (MQ=255)
                                                                                                    gTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGt                                                >  1:5015341/1‑150 (MQ=255)
                                                                                                    gTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGt                                                >  1:1306610/1‑150 (MQ=255)
                                                                                                     tACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTa                                               >  2:3177339/1‑150 (MQ=255)
                                                                                                       cccGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTAtt                                             <  2:1053526/150‑1 (MQ=255)
                                                                                                       cccGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTAtt                                             <  2:3064069/150‑1 (MQ=255)
                                                                                                       cccGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTAtt                                             <  2:4574690/150‑1 (MQ=255)
                                                                                                        ccGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCgtgggtgg                                                                                                        <  2:869378/90‑1 (MQ=255)
                                                                                                          gAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCa                                          >  1:4939749/1‑150 (MQ=255)
                                                                                                              gTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGt                                      <  1:2348145/150‑1 (MQ=255)
                                                                                                              gTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGt                                      <  1:623534/150‑1 (MQ=255)
                                                                                                              gTCAGCCACAGGTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGt                                      <  2:3231387/150‑1 (MQ=255)
                                                                                                                cAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTgg                                    <  2:5160583/150‑1 (MQ=255)
                                                                                                                cAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTgg                                    <  1:493605/150‑1 (MQ=255)
                                                                                                                 aggcattgTTTCATCCGTCGCGAAATTCTGGCGCTGGAGCGCCTAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGAGGCCGAAGTGCAGGTTGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTggg                                   <  2:952051/143‑1 (MQ=255)
                                                                                                                  gCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGc                                  >  2:842974/1‑150 (MQ=255)
                                                                                                                   ccACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCt                                 >  2:2034868/1‑150 (MQ=255)
                                                                                                                    cacaGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGATGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTc                                >  2:5079612/1‑150 (MQ=255)
                                                                                                                      caGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCaa                              <  2:4939749/150‑1 (MQ=255)
                                                                                                                        gTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAgg                            <  2:2483408/150‑1 (MQ=255)
                                                                                                                         tttCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAggg                           <  2:1758467/150‑1 (MQ=255)
                                                                                                                          ttCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAgggg                          <  1:1588594/150‑1 (MQ=255)
                                                                                                                              tCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGgttgt                      <  2:3755232/150‑1 (MQ=255)
                                                                                                                               ccGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGgttgtt                     <  2:186948/150‑1 (MQ=255)
                                                                                                                                 gTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGGTTGTTgg                   <  2:1735701/150‑1 (MQ=255)
                                                                                                                                 gTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGGTTGTTgg                   >  1:414352/1‑150 (MQ=255)
                                                                                                                                    gcgAAATACTGGCGCTGGAGCGCCAAGGGGATGAGCAGCAACGAATAGCCCTGCGTGGTTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGAAAGGTGGGCTCAGGGGGTTGTTGGCAc                >  2:1868526/1‑150 (MQ=255)
                                                                                                                                       aaaTACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGGTTGTTGGCACCGa             >  2:1310415/1‑150 (MQ=255)
                                                                                                                                         cccttGGCGCTCCAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGGTTGTTGGCACCGACg           >  2:1344037/5‑150 (MQ=255)
                                                                                                                                          tACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGGTTGTTGGCACCGACGt          <  1:1132403/150‑1 (MQ=255)
                                                                                                                                           cttGGCGCTCCAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTAtt                                             >  1:3098253/3‑114 (MQ=255)
                                                                                                                                            ttGGCGCTCCAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGgttgtt                     >  2:3571825/2‑137 (MQ=255)
                                                                                                                                             tGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACgggg                                                               >  1:3120060/1‑94 (MQ=255)
                                                                                                                                             tGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACgggg                                                               <  2:3120060/94‑1 (MQ=255)
                                                                                                                                               gcgcTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGGTTGTTGGCACCGACGTGGCAg     >  1:2016976/1‑150 (MQ=255)
                                                                                                                                               gcgcTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGGTTGTTGGCACCGACGTGGCAg     >  2:5227472/1‑150 (MQ=255)
                                                                                                                                               gcgcTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGGTTGTGGGCACCGACGTGGCAg     >  1:2016436/1‑150 (MQ=255)
                                                                                                                                                 gcTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGGTTGTTGGCACCGACGTGGCAGGt   <  1:1032856/150‑1 (MQ=255)
                                                                                                                                                  cTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGGTTGTTGGCACCGACGTGGCAGGTa  <  1:2775708/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GCCCGTAGTCGACAGCCTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGGTTGTTGGCACCGACGTGGCAGGTA  >  NC_002947/3556276‑3556571

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