Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,227,613 C→A 38.6% N324K (AAC→AAA umuC → DNA polymerase V, subunit C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,227,6130CA38.6% 325.7 / 270.8 246N324K (AAC→AAAumuCDNA polymerase V, subunit C
Reads supporting (aligned to +/- strand):  ref base C (74/77);  new base A (54/41);  total (128/118)
Fisher's exact test for biased strand distribution p-value = 2.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 248 total aligned reads displayed.

CGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTTCAGTCAGGGAGTCGCGCAGCTCAATTTATTCGATGACAACGCACCGCGCCCCGGGAGTGAGCAATTGAT  >  CP009273/1227465‑1227762
                                                                                                                                                    |                                                                                                                                                     
cggAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACg                                                                                                                                                      <  1:2840285/150‑1 (MQ=255)
cggAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACg                                                                                                                                                      >  1:909084/1‑150 (MQ=255)
cggAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACg                                                                                                                                                      <  1:162129/150‑1 (MQ=255)
 ggAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACgc                                                                                                                                                     <  1:3361274/150‑1 (MQ=255)
 ggAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACgc                                                                                                                                                     >  2:2399298/1‑150 (MQ=255)
 ggAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAAAgc                                                                                                                                                     >  1:3610887/1‑150 (MQ=255)
    aaaaCTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAAAgctgc                                                                                                                                                  <  1:2137890/150‑1 (MQ=255)
    aaaaCTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCAAGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAAAgctgc                                                                                                                                                  <  1:1556532/150‑1 (MQ=255)
     aaaCTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACgctgct                                                                                                                                                 >  1:1994618/1‑150 (MQ=255)
     aaaCTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAAAgctgct                                                                                                                                                 <  1:1781340/150‑1 (MQ=255)
      aaCTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTa                                                                                                                                                <  1:307438/150‑1 (MQ=255)
      aaCTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTa                                                                                                                                                <  1:591805/150‑1 (MQ=255)
        cTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTAcg                                                                                                                                              <  1:3496865/150‑1 (MQ=255)
             cAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAAAGCTGCTACGCGAtc                                                                                                                                         <  2:368606/150‑1 (MQ=255)
              aGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTCCGCGAtct                                                                                                                                        >  2:2270211/1‑150 (MQ=255)
              aGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTACGCGAtct                                                                                                                                        >  1:3760221/1‑150 (MQ=255)
                cGAGCATCACTATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAAAGCTGCTACGCGATctct                                                                                                                                      <  2:1897614/150‑1 (MQ=255)
                 gAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTg                                                                                                                                     >  2:2598029/1‑150 (MQ=255)
                    cATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAAAGCTGCTACGCGATCTCTGGAt                                                                                                                                  >  2:3347388/1‑150 (MQ=255)
                    cATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAAAGCTGCTACGCGATCTCTGGAt                                                                                                                                  <  2:3557357/150‑1 (MQ=255)
                     aTCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAAAGCTGCTACGCGATCTCTGGATg                                                                                                                                 >  1:2360697/1‑150 (MQ=255)
                       cAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAAAGCTGCTACGCGATCTCTGGATGcc                                                                                                                               >  1:2183897/1‑150 (MQ=255)
                          tatTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATc                                                                                                                            >  1:1112495/1‑150 (MQ=255)
                           atTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCt                                                                                                                           >  1:2786892/1‑150 (MQ=255)
                           atTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAAAGCTGCTACGCGATCTCTGGATGCCATCt                                                                                                                           >  1:2908809/1‑150 (MQ=255)
                            ttGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTTCGGCAATAGCGCGTCTGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTg                                                                                                                          <  1:791214/150‑1 (MQ=255)
                             tGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTgg                                                                                                                         <  1:2436210/150‑1 (MQ=255)
                              gTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAAAGCTGCTACGCGATCTCTGGATGCCATCTGGc                                                                                                                        <  2:1317441/150‑1 (MQ=255)
                              gTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAAAGCTGCTACGCGATCTCTGGATGCCATCTGGc                                                                                                                        >  1:3534253/1‑150 (MQ=255)
                                cGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCaa                                                                                                                      >  2:4097553/1‑150 (MQ=255)
                                cGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAAAGCTGCTACGCGATCTCTGGATGCCATCTGGCaa                                                                                                                      >  2:2879701/1‑150 (MQ=255)
                                   tttATCTCCACGTTTATTAAGACGTCACCCTTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGATACGCGATCTATGGATGCCATCTGGCTAGCg                                                                                                                   <  2:2659186/150‑1 (MQ=255)
                                   tttATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCg                                                                                                                   <  1:2544803/150‑1 (MQ=255)
                                   tttATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAAATGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCg                                                                                                                   <  2:212656/150‑1 (MQ=255)
                                    ttATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCgg                                                                                                                  <  2:1976900/150‑1 (MQ=255)
                                       tctcCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGcc                                                                                                               <  1:3977269/150‑1 (MQ=255)
                                         tcCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCAt                                                                                                             <  1:1766734/150‑1 (MQ=255)
                                         tcCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAAAGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCAt                                                                                                             <  1:934412/150‑1 (MQ=255)
                                           cACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCg                                                                                                           <  2:3108903/150‑1 (MQ=255)
                                              gTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAAAGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTa                                                                                                        >  2:846267/1‑150 (MQ=255)
                                              gTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAAAGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTa                                                                                                        <  1:2951017/150‑1 (MQ=255)
                                               tttATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTAc                                                                                                       >  1:3844637/1‑150 (MQ=255)
                                                     aaGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAAAGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCaaaaa                                                                                                 >  2:822855/1‑150 (MQ=255)
                                                      aGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAAAGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAg                                                                                                <  1:2661576/150‑1 (MQ=255)
                                                      aGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAAAGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAg                                                                                                <  1:3672793/150‑1 (MQ=255)
                                                       gACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGc                                                                                               <  2:670232/150‑1 (MQ=255)
                                                       gACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAAAGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGc                                                                                               <  1:717630/150‑1 (MQ=255)
                                                          gTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCggg                                                                                            >  1:1746423/1‑150 (MQ=255)
                                                           tCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAAAGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGc                                                                                           <  2:4116552/150‑1 (MQ=255)
                                                                  ttGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGc                                                                                    >  2:1520445/1‑150 (MQ=255)
                                                                     cgcTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTgg                                                                                 >  2:368237/1‑150 (MQ=255)
                                                                       cTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTgggg                                                                               <  2:3181292/150‑1 (MQ=255)
                                                                         cAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGa                                                                             >  1:3903785/1‑150 (MQ=255)
                                                                         cAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGa                                                                             >  1:84632/1‑150 (MQ=255)
                                                                          aaTGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGAt                                                                            >  2:1555570/1‑150 (MQ=255)
                                                                           aTGAACCTTAGTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGAtt                                                                           <  2:3020702/150‑1 (MQ=255)
                                                                              aaCCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATttct                                                                        <  2:1081416/150‑1 (MQ=255)
                                                                              aaCCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATttct                                                                        <  2:3405725/150‑1 (MQ=255)
                                                                                ccTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATttcttc                                                                      <  1:1206070/150‑1 (MQ=255)
                                                                                  ttattaCGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTtcag                                                                    <  2:2810838/150‑1 (MQ=255)
                                                                                       aCGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAAAGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTTCAGTCAgg                                                               <  2:4079130/150‑1 (MQ=255)
                                                                                        cGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTTCAGTCAggg                                                              >  1:499843/1‑150 (MQ=255)
                                                                                          gCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAAAGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTTCAGTCAGGGAg                                                            <  1:4122241/150‑1 (MQ=255)
                                                                                           cAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAAAGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTTCAGTCAGGGAGt                                                           >  2:2935566/1‑150 (MQ=255)
                                                                                            aaTAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTTCAGTCAGGGAGTc                                                          <  2:134661/150‑1 (MQ=255)
                                                                                            aaTAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAAAGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTTCAGTCAGGGAGTc                                                          <  1:57452/150‑1 (MQ=255)
                                                                                                   cgTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAAAGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTTCAGTCAGGGAGTCGCGCAGc                                                   >  1:1116169/1‑150 (MQ=255)
                                                                                                      cGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAAAGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTTCAGTCAGGGAGTCGCGCAGCTCa                                                >  2:2346357/1‑150 (MQ=255)
                                                                                                       ggTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAAAGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTTCAGTCAGGGAGTcgcg                                                       <  2:2945537/142‑1 (MQ=255)
                                                                                                       ggTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAAAGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTTCAGTCAGGGAGTCGCGCAGCTCaa                                               >  2:2118747/1‑150 (MQ=255)
                                                                                                         tAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTTCAGTCAGGGAGTCGCGCAGCTCAAtt                                             >  2:1789606/1‑150 (MQ=255)
                                                                                                          aaaaCTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAAAGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTTCAGTCAGGGAGTCGCGCAGCTCAAttt                                            >  1:3668724/1‑150 (MQ=255)
                                                                                                             aCTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTTCAGTCAGGGAGTCGCGCAGCTCAATTTAtt                                         >  1:3323056/1‑150 (MQ=255)
                                                                                                              ctgctgACGCCCACTCAGGACAGCAGGGATATCATTAAAGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTTCAGTCAGGGAGTCGCGCAGCTCAATTTATTc                                        <  1:1344119/150‑1 (MQ=255)
                                                                                                                 ctgACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCTGGCGTGATGCTGGGGGATTTCTTCAGTCAGGGAGTCGCGCAGCTCAATTTATTCGAt                                     <  2:965507/150‑1 (MQ=255)
                                                                                                                 ctgACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTTCAGTCAGGGAGTCGCGCAGCTCAATTTATTCGAt                                     >  2:1692082/1‑150 (MQ=255)
                                                                                                                   gACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTTCAGTCAGGGAGTCGCGCAGCTCAATTTATTCGATGa                                   <  1:2159316/150‑1 (MQ=255)
                                                                                                                   gACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTTCAGTCAGGGAGTCGCGCAGCTCAATTTATTCGATGa                                   <  1:3467488/150‑1 (MQ=255)
                                                                                                                     cGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTTCAGTCAGGGAGTCGCGCAGCTCAATTTATTCGATGACa                                 <  2:2517902/150‑1 (MQ=255)
                                                                                                                      gCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTTCAGTCAGGGAGTCGCGCAGCTCAATTTATTCGATGACaa                                >  2:1004117/1‑150 (MQ=255)
                                                                                                                       cccACTCAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTTCAGTCAGGGAGTCGCGCAGCTCAATTTATTCGATGACAAc                               <  1:2345654/150‑1 (MQ=255)
                                                                                                                       cccACTCAGGACAGCAGGGATATCATTAAAGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTTCAGTCAGGGAGTCGCGCAGCTCAATTTATTCGATGACAAc                               <  2:1194014/150‑1 (MQ=255)
                                                                                                                          aCTCAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTTCAGTCAGGGAGTCGCGCAGCTCAATTTATTCGATGACAACGCa                            <  1:2208177/150‑1 (MQ=255)
                                                                                                                            tCAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTTCAGTCAGGGAGTCGCGCAGCTCAATTTATTCGATGACAACGCAcc                          >  2:4012765/1‑150 (MQ=255)
                                                                                                                             cAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTTCAGTCAGGGAGTCGCGCAGCTCAATTTATTCGATGACAACGCACcg                         <  2:2973026/150‑1 (MQ=255)
                                                                                                                             cAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTTCAGTCAGGGAGTCGCGCAGCTCAATTTATTCGATGACAACGCACcg                         >  2:3824178/1‑150 (MQ=255)
                                                                                                                              aGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTTCAGTCAGGGAGTCGCGCAGCTCAATTTATTCGATGACAACGCACcgc                        >  1:2423584/1‑150 (MQ=255)
                                                                                                                               ggACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTTCAGTCAGGGAGTCGCGCAGCTCAATTTATTCGATGACAACGCACcgcg                       <  1:3933228/150‑1 (MQ=255)
                                                                                                                               ggACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTTCAGTCAGGGAGTCGCGCAGCTCAATTTATTCGATGACAACGCACcgcg                       >  1:3819075/1‑150 (MQ=255)
                                                                                                                                  cagcagGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTTCAGTCAGGGAGTCGCGCAGCTCAATTTATTCGATGACAACGCACCGCGccc                    <  1:2600808/150‑1 (MQ=255)
                                                                                                                                   agcagGGATATCATTAAAGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTTCAGTCAGGGAGTCGCGCAGCTCAATTTATTCGATGACAACGCACCGCGcccc                   <  1:3517732/150‑1 (MQ=255)
                                                                                                                                    gcagGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTTCAGTCAGGGAGTCGCGCAGCTCAATTTATTCGATGACAACGCACCGCGCCCCg                  <  1:2453777/150‑1 (MQ=255)
                                                                                                                                     cagGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTTCAGTCAGGGAGTCGCGCAGCTCAATTTATTCGATGACAACGCACCGCGCCCCgg                 >  1:3594921/1‑150 (MQ=255)
                                                                                                                                      agGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTTCAGTCAGGGAGTCGCGCAGCTCAATTTATTCGATGACAACGCACCGCGCCCCggg                <  1:2737598/150‑1 (MQ=255)
                                                                                                                                      agGGATATCATTAAAGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTTCAGTCAGGGAGTCGCGCAGCTCAATTTATTCGATGACAACGCACCGCGCCCCggg                >  1:4081827/1‑150 (MQ=255)
                                                                                                                                        ggATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTTCAGTCAGGGAGTCGCGCAGCTCAATTTATTCGATGACAACGCACCGCGCCCCGGGAg              <  2:349163/150‑1 (MQ=255)
                                                                                                                                         gATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTTCAGTCAGGGAGTCGCGCAGCTCAATTTATTCGATGACAACGCACCGCGCCCCGGGAGt             >  2:3741339/1‑150 (MQ=255)
                                                                                                                                          atatCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTTCAGTCAGGGAGTCGCGCAGCTCAATTTATTCGATGACAACGCACCGCGCCCCGGGAGTg            <  1:2924527/150‑1 (MQ=255)
                                                                                                                                              cATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTTCAGTCAGGGAGTCGCGCAGCTCAATTTATTCGATGACAACGCACCGCGCCCCGGGAGTGAGCa        <  1:1883530/150‑1 (MQ=255)
                                                                                                                                                    cGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTTCAGTCAGGGAGTCGCGCAGCTCAATTTATTCGATGACAACGCACCGCGCCCCGGGAGTGAGCAATtgat  <  2:2462844/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGTAAAACTGCTGACGCCCACTCAGGACAGCAGGGATATCATTAACGCTGCTACGCGATCTCTGGATGCCATCTGGCAAGCGGGCCATCGTTACCAAAAAGCGGGCGTGATGCTGGGGGATTTCTTCAGTCAGGGAGTCGCGCAGCTCAATTTATTCGATGACAACGCACCGCGCCCCGGGAGTGAGCAATTGAT  >  CP009273/1227465‑1227762

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: