Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,646,327 G→T 9.1% D159Y (GAT→TAT)  sseA → 3‑mercaptopyruvate sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,646,3270GT 9.1% 854.9 / 39.7 237D159Y (GAT→TAT) sseA3‑mercaptopyruvate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base G (112/103);  new base T (11/10);  total (124/113)
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.97e-01

Only 100 of 246 total aligned reads displayed.

GCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTG  >  CP009273/2646178‑2646476
                                                                                                                                                     |                                                                                                                                                     
gcACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCg                                                                                                                                                       <  1:1545722/150‑1 (MQ=255)
 cACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGa                                                                                                                                                      >  1:3064903/1‑150 (MQ=255)
  aCCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGAt                                                                                                                                                     <  2:1202374/150‑1 (MQ=255)
   ccTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATg                                                                                                                                                    <  1:1679178/150‑1 (MQ=255)
   ccTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATg                                                                                                                                                    <  1:1652081/150‑1 (MQ=255)
         gtgtAGAGAAAGTGTCGATTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCTATGtattat                                                                                                                                              <  2:4157212/150‑1 (MQ=255)
            tAGAGAAAGTGTCGATTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTgg                                                                                                                                           >  1:2144545/1‑150 (MQ=255)
            tAGAGAAAGTGTCGATTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTgg                                                                                                                                           >  2:1263047/1‑150 (MQ=255)
             agagAAAGTGTCGATTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGc                                                                                                                                          >  2:3827338/1‑150 (MQ=255)
               agaAAGTGTCGATTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCaa                                                                                                                                        >  2:3103685/1‑150 (MQ=255)
               agaAAGTGTCGATTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCaa                                                                                                                                        <  2:747963/150‑1 (MQ=255)
                gaAAGTGTCGAGTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAg                                                                                                                                       <  2:3340495/150‑1 (MQ=255)
                 aaaGTGTCGATTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGc                                                                                                                                      >  1:3436087/1‑150 (MQ=255)
                 aaaGTGTCGATTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGc                                                                                                                                      >  1:1382060/1‑150 (MQ=255)
                  aaGTGTCGATTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGcc                                                                                                                                     <  1:3965750/150‑1 (MQ=255)
                   aGTGTCGATTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCa                                                                                                                                    <  2:1807621/150‑1 (MQ=255)
                    gtgtCGATTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCAt                                                                                                                                   <  2:362230/150‑1 (MQ=255)
                     tgtCGATTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCTATGTATTATTGGCAAGCCATg                                                                                                                                  >  2:3450466/1‑150 (MQ=255)
                          aTTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGATAGTATAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAAt                                                                                                                             >  2:997381/1‑150 (MQ=255)
                             cTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACg                                                                                                                          >  2:666848/1‑150 (MQ=255)
                              tGGGGGGTTGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTATAGCTGCCGGAAGGAGAGTTTAAAGACGCGATTAATCCTGAAGCCGTGGTGAAAGTAACCAATGTATTATTTGCAAGCCATGAAAATActg                                                                                                                         >  1:3190281/1‑148 (MQ=255)
                              tGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCTATGTATTATTGGCAAGCCATGAAAATACgg                                                                                                                         <  2:2933690/150‑1 (MQ=255)
                              tGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCTATGTATTATTGGCAAGCCATGAAAATACgg                                                                                                                         >  1:907410/1‑150 (MQ=255)
                              tGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCTATGTATTATTGGCAAGCCATGAAAATACgg                                                                                                                         <  2:651765/150‑1 (MQ=255)
                              tGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACgg                                                                                                                         >  2:2186938/1‑150 (MQ=255)
                              tGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACgg                                                                                                                         >  1:1792456/1‑150 (MQ=255)
                              tGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACgg                                                                                                                         >  2:2186907/1‑150 (MQ=255)
                               ggggggTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGgc                                                                                                                        <  2:1912584/150‑1 (MQ=255)
                                 ggggTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGgcgc                                                                                                                      <  1:1599232/150‑1 (MQ=255)
                                 ggggTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGgcgc                                                                                                                      <  1:3860233/150‑1 (MQ=255)
                                     tGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATAGGGCGCAAat                                                                                                                  >  2:1120125/1‑150 (MQ=255)
                                     tGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAat                                                                                                                  <  2:3275488/150‑1 (MQ=255)
                                     tGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAat                                                                                                                  <  2:2975300/150‑1 (MQ=255)
                                      ggACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAatt                                                                                                                 >  1:300282/1‑150 (MQ=255)
                                          ttGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTg                                                                                                             >  1:3767552/1‑150 (MQ=255)
                                            gCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGAt                                                                                                           <  2:1446407/150‑1 (MQ=255)
                                             cAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATg                                                                                                          >  2:1309304/1‑150 (MQ=255)
                                                    gCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCTATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATgcccgccc                                                                                                   <  2:4400530/150‑1 (MQ=255)
                                                      aGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCTATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCgg                                                                                                 <  1:1359954/150‑1 (MQ=255)
                                                          cgATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGc                                                                                             >  1:2375364/1‑150 (MQ=255)
                                                            atgatCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCAc                                                                                           <  1:1438396/150‑1 (MQ=255)
                                                               atCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGtt                                                                                        >  2:1472923/1‑150 (MQ=255)
                                                               atCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGtt                                                                                        >  1:122338/1‑150 (MQ=255)
                                                                tCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGttt                                                                                       <  1:986792/150‑1 (MQ=255)
                                                                 ctgctgTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGtttt                                                                                      <  2:4518985/150‑1 (MQ=255)
                                                                   gctgTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTaa                                                                                    >  1:46377/1‑150 (MQ=255)
                                                                       ttAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCa                                                                                >  1:2621658/1‑150 (MQ=255)
                                                                        tAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAg                                                                               >  2:3943712/1‑150 (MQ=255)
                                                                            agaagGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGt                                                                           <  1:3155523/150‑1 (MQ=255)
                                                                            agaagGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGt                                                                           >  1:2001408/1‑150 (MQ=255)
                                                                             gaagGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGtt                                                                          <  2:2321355/150‑1 (MQ=255)
                                                                              aagGTGCAGTAGAGCTGCCGGAATGAGAGTTTAACGCCGCTTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTtg                                                                         <  2:2446553/150‑1 (MQ=255)
                                                                              aagGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTtg                                                                         <  1:1472923/150‑1 (MQ=255)
                                                                              aagGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTtg                                                                         <  2:3005558/150‑1 (MQ=255)
                                                                               agGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTtga                                                                        >  2:2846388/1‑150 (MQ=255)
                                                                                ggTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTtgat                                                                       <  1:3561298/150‑1 (MQ=255)
                                                                                   gCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGaa                                                                    <  1:1350701/150‑1 (MQ=255)
                                                                                    cAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAAc                                                                   >  2:2268526/1‑150 (MQ=255)
                                                                                     aGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAAcc                                                                  >  1:1235531/1‑150 (MQ=255)
                                                                                      gTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCt                                                                 <  1:2413649/150‑1 (MQ=255)
                                                                                        agagCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCg                                                               >  2:632324/1‑150 (MQ=255)
                                                                                         gagCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGc                                                              >  2:1809600/1‑150 (MQ=255)
                                                                                             tGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAg                                                          >  1:515302/1‑150 (MQ=255)
                                                                                              gCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAgg                                                         >  1:864238/1‑150 (MQ=255)
                                                                                                cGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGtt                                                       <  1:3319997/150‑1 (MQ=255)
                                                                                                 ggAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGttt                                                      <  2:259732/150‑1 (MQ=255)
                                                                                                    aGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACg                                                   >  2:2760046/1‑150 (MQ=255)
                                                                                                     ggAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGt                                                  >  2:2960486/1‑150 (MQ=255)
                                                                                                     ggAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGt                                                  >  2:1482288/1‑150 (MQ=255)
                                                                                                      gagaGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTc                                                 <  2:4077462/150‑1 (MQ=255)
                                                                                                         agTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTcgcg                                              >  2:95632/1‑150 (MQ=255)
                                                                                                          gTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCgg                                             <  2:3897553/150‑1 (MQ=255)
                                                                                                            ttAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGAc                                           <  2:1882722/150‑1 (MQ=255)
                                                                                                             tAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACa                                          <  1:1285710/150‑1 (MQ=255)
                                                                                                              aaCGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCTATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACat                                         >  2:197478/1‑150 (MQ=255)
                                                                                                              aaCGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCTATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACat                                         >  2:197452/1‑150 (MQ=255)
                                                                                                                     cgTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCg                                  <  2:3387566/150‑1 (MQ=255)
                                                                                                                      gTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCgg                                 <  2:3651776/150‑1 (MQ=255)
                                                                                                                      gTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCgg                                 >  1:2325583/1‑150 (MQ=255)
                                                                                                                            tCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACt                           >  2:3672449/1‑150 (MQ=255)
                                                                                                                            tCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACt                           >  1:422805/1‑150 (MQ=255)
                                                                                                                            tCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACt                           <  1:2311786/150‑1 (MQ=255)
                                                                                                                             ccTCAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTg                          >  2:3844910/1‑150 (MQ=255)
                                                                                                                               tGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGaa                        >  1:2275335/1‑150 (MQ=255)
                                                                                                                                    ccGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTc                   >  1:1056931/1‑150 (MQ=255)
                                                                                                                                    ccGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTc                   >  2:1148489/1‑150 (MQ=255)
                                                                                                                                     cGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTcc                  <  1:4515789/150‑1 (MQ=255)
                                                                                                                                        ggtgAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTg               >  2:2504973/1‑150 (MQ=255)
                                                                                                                                        ggtgAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTg               >  1:2645912/1‑150 (MQ=255)
                                                                                                                                         gtgAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTgg              <  2:4268493/150‑1 (MQ=255)
                                                                                                                                           gAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGAc            <  2:3432075/150‑1 (MQ=255)
                                                                                                                                              aGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGa         >  2:1570087/1‑150 (MQ=255)
                                                                                                                                              aGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGa         >  2:3404625/1‑150 (MQ=255)
                                                                                                                                                tAACCGATGTATTATTGGCAAGCCAAGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACAGAATGATCCGTGGACGGAAc       >  1:2494726/1‑150 (MQ=255)
                                                                                                                                                 aaCCTATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACt      <  2:1414357/150‑1 (MQ=255)
                                                                                                                                                 aaCCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTTCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGATTACGTCGCGGACATATACACGGAGCACTGAATGTTCCGTGGACGGAACt      >  1:2888443/1‑150 (MQ=255)
                                                                                                                                                 aaCCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACt      <  1:367024/150‑1 (MQ=255)
                                                                                                                                                     gATGTATTATTGGCAATCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTg  <  1:4408401/150‑1 (MQ=255)
                                                                                                                                                     gATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTg  >  1:4007663/1‑150 (MQ=255)
                                                                                                                                                     gATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTGACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTg  <  2:3666376/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAATCCTGAAGCCGTGGTGAAAGTAACCGATGTATTATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTG  >  CP009273/2646178‑2646476

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: