Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,302,650 C→A 16.2% M57I (ATG→ATT alkB ← oxidative demethylase of N1‑methyladenine or N3‑methylcytosine DNA lesions

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,302,6500CA16.2% 897.9 / 120.3 299M57I (ATG→ATTalkBoxidative demethylase of N1‑methyladenine or N3‑methylcytosine DNA lesions
Reads supporting (aligned to +/- strand):  ref base C (121/129);  new base A (23/26);  total (144/155)
Fisher's exact test for biased strand distribution p-value = 8.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 300 total aligned reads displayed.

CGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGCTCTTGCCACGGTTCAGC  >  CP009273/2302503‑2302799
                                                                                                                                                   |                                                                                                                                                     
cgccgTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACAt                                                                                                                                                     <  1:2408217/150‑1 (MQ=255)
  ccgTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATgg                                                                                                                                                   >  1:3056648/1‑150 (MQ=255)
   cgTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGt                                                                                                                                                  <  2:4316924/150‑1 (MQ=255)
    gTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGta                                                                                                                                                 <  2:4658431/150‑1 (MQ=255)
           cACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATcccc                                                                                                                                          <  1:4255297/150‑1 (MQ=255)
             cGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCg                                                                                                                                        <  1:2617846/150‑1 (MQ=255)
               ttGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCggg                                                                                                                                      >  1:4669225/1‑150 (MQ=255)
               ttGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCggg                                                                                                                                      <  2:4966628/150‑1 (MQ=255)
               ttGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCggg                                                                                                                                      >  2:436965/1‑150 (MQ=255)
               ttGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCcgg                                                                                                                                      >  1:5051061/1‑150 (MQ=255)
                tGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCgggg                                                                                                                                     >  2:4786974/1‑150 (MQ=255)
                   cATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATCTCCCCCGGGGGTg                                                                                                                                  >  2:1785088/1‑150 (MQ=255)
                   cATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTg                                                                                                                                  >  1:2125424/1‑150 (MQ=255)
                   cATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTg                                                                                                                                  >  2:2667930/1‑150 (MQ=255)
                    aTAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGa                                                                                                                                 <  2:5204233/150‑1 (MQ=255)
                    aTAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGa                                                                                                                                 >  1:483821/1‑150 (MQ=255)
                     tAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGAc                                                                                                                                <  1:1359521/150‑1 (MQ=255)
                       aaTTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCa                                                                                                                              >  1:5174456/1‑150 (MQ=255)
                          tATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCt                                                                                                                           <  1:1270838/150‑1 (MQ=255)
                           aTGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTg                                                                                                                          >  1:395528/1‑150 (MQ=255)
                           aTGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTg                                                                                                                          <  1:2112288/150‑1 (MQ=255)
                            tGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTgg                                                                                                                         >  2:1165654/1‑150 (MQ=255)
                             gAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGc                                                                                                                        <  1:3791324/150‑1 (MQ=255)
                              aaaaCTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCg                                                                                                                       >  1:4783519/1‑150 (MQ=255)
                              aaaaCTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCg                                                                                                                       <  1:4328320/150‑1 (MQ=255)
                              aaaaCTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCg                                                                                                                       <  1:1720288/150‑1 (MQ=255)
                                  ctctGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAAc                                                                                                                   <  1:2595225/150‑1 (MQ=255)
                                      gtgGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCg                                                                                                               <  2:925234/150‑1 (MQ=255)
                                          cATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGActg                                                                                                           >  2:1958868/1‑150 (MQ=255)
                                           aTGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGActgg                                                                                                          >  2:2283797/1‑150 (MQ=255)
                                              gCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGActggctg                                                                                                       <  1:623422/150‑1 (MQ=255)
                                                gggCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACtggctggc                                                                                                     >  1:2075758/1‑150 (MQ=255)
                                                gggCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACtggctggc                                                                                                     >  1:4208688/1‑150 (MQ=255)
                                                gggCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACtggctggc                                                                                                     <  2:2380828/150‑1 (MQ=255)
                                                 ggCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCa                                                                                                    >  1:3571933/1‑150 (MQ=255)
                                                 ggCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCa                                                                                                    >  1:948415/1‑150 (MQ=255)
                                                    cACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACg                                                                                                 >  2:4476505/1‑150 (MQ=255)
                                                    cACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACg                                                                                                 >  2:3175655/1‑150 (MQ=255)
                                                    cACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACg                                                                                                 >  1:2725020/1‑150 (MQ=255)
                                                    cACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACg                                                                                                 >  1:2321057/1‑150 (MQ=255)
                                                         tttATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCAtt                                                                                            <  2:1718451/150‑1 (MQ=255)
                                                          ttATTTGTTTGCGGTTCTCTGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATTTCCCCCGGGGGTGACCATCTGGCGTAACGGCGACTGGGTGGCAAGGTGAttt                                                                                           <  2:1642283/150‑2 (MQ=255)
                                                          ttATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTa                                                                                           >  1:4090778/1‑150 (MQ=255)
                                                          ttATTTGTTTGCGGATCAAGGGGCGAAGAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTa                                                                                           <  2:3117984/150‑1 (MQ=255)
                                                           tATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTaa                                                                                          >  1:4651093/1‑150 (MQ=255)
                                                             tttgtttgCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAata                                                                                        <  2:3817147/150‑1 (MQ=255)
                                                              ttgtttgCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAatat                                                                                       <  2:5079082/150‑1 (MQ=255)
                                                              ttgtttgCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCATCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAatat                                                                                       <  2:1396618/150‑1 (MQ=255)
                                                                gtttgCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATcg                                                                                     <  1:4419081/150‑1 (MQ=255)
                                                                  ttgCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATcgcg                                                                                   <  1:2256257/150‑1 (MQ=255)
                                                                    gCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGa                                                                                 >  1:2023810/1‑150 (MQ=255)
                                                                     cGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGAt                                                                                <  1:1588590/150‑1 (MQ=255)
                                                                     cGGATCAATGGGCGAATAGAGAGAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGAt                                                                                <  2:411278/150‑1 (MQ=255)
                                                                      ggATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATc                                                                               <  1:527457/150‑1 (MQ=255)
                                                                      ggATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATc                                                                               <  1:3337716/150‑1 (MQ=255)
                                                                      ggATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACAATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATc                                                                               >  1:1864043/1‑150 (MQ=255)
                                                                       gATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCa                                                                              >  1:3084416/1‑150 (MQ=255)
                                                                          cAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGtt                                                                           <  1:3020418/150‑1 (MQ=255)
                                                                           aaTGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTg                                                                          >  1:2278268/1‑150 (MQ=255)
                                                                               ggCGAATAGAGAGAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTcc                                                                      <  2:2974827/150‑1 (MQ=255)
                                                                                  gAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCa                                                                   >  1:8977/1‑150 (MQ=255)
                                                                                   aaTAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCAGTAATATCGCGGATCAGTTGCTCCGCAg                                                                  >  1:4870698/1‑150 (MQ=255)
                                                                                    aTAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCCCAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTCATATCGCGGATCAGTTGCTCCGCAGc                                                                 >  2:4917158/1‑150 (MQ=255)
                                                                                      agagATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGt                                                               <  1:5223824/150‑1 (MQ=255)
                                                                                       gtgaTAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGtt                                                              <  2:4710122/148‑1 (MQ=255)
                                                                                       gagaTAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGtt                                                              >  1:1908958/1‑150 (MQ=255)
                                                                                       gagaTAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGtt                                                              <  1:1532716/150‑1 (MQ=255)
                                                                                         gaTAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTaa                                                            >  1:1228563/1‑150 (MQ=255)
                                                                                           tatccTTGCCGATGGGTCGTCCAGCCCAGATGCCCACTGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACTGCGACTGGCTGGCATCGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTaaaa                                                          <  2:4625991/147‑1 (MQ=255)
                                                                                            aaCCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATCTCGCGGATCAGTTGCTCCGCAGCGTTAAAAg                                                         >  2:4903000/1‑150 (MQ=255)
                                                                                             aCCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGc                                                        <  1:1009494/150‑1 (MQ=255)
                                                                                               ctagCCGATGGGTCGTCCAGCCCAGATGCCCACAGGTGGTCATCGCCCCCGAAAGGGAATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTCGCTGGCAACGTCCTTCATATCGCGGATCAGGTGCGCCGCAGCGTGACAAGcta                                                      >  2:5240038/4‑148 (MQ=255)
                                                                                                 tGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCaaaa                                                    >  2:3491991/1‑150 (MQ=255)
                                                                                                     gATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAATGCAAAACGcc                                                >  2:4872925/1‑150 (MQ=255)
                                                                                                          gtcgtcCAGCCCAGATGACCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACAATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTaaa                                           >  2:155606/1‑150 (MQ=255)
                                                                                                           tcgtcCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTaaaa                                          <  1:1164249/150‑1 (MQ=255)
                                                                                                             gtcCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAAtt                                        <  2:115031/150‑1 (MQ=255)
                                                                                                             gtcCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAAtt                                        >  1:2811048/1‑150 (MQ=255)
                                                                                                             gtcCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAAtt                                        >  1:3489648/1‑150 (MQ=255)
                                                                                                             gtcCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAAtt                                        <  1:4555023/150‑1 (MQ=255)
                                                                                                                 aGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCg                                    <  1:4704767/150‑1 (MQ=255)
                                                                                                                     cAGTTGCCCACAGTTGGTCATCTCCACCGACATGGTATATCCCCCGTGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACTTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCAGTAAAATTACCGCAcc                                <  2:3032531/150‑1 (MQ=255)
                                                                                                                     cAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCAcc                                <  2:1188499/150‑1 (MQ=255)
                                                                                                                      aGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCa                               >  2:2068068/1‑150 (MQ=255)
                                                                                                                      aGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCa                               >  2:747149/1‑150 (MQ=255)
                                                                                                                        aTGCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGGTAAAAGCAAAACGCCGTACAATTACCGCACCAgc                             >  2:867199/1‑150 (MQ=255)
                                                                                                                            ccACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAgccgcc                         <  2:2125424/150‑1 (MQ=255)
                                                                                                                            ccACAGTTGGTCATCGCCACCGACATGGTATAACCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAgccgcc                         >  1:3759637/1‑150 (MQ=255)
                                                                                                                                aGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTg                     <  1:4786974/150‑1 (MQ=255)
                                                                                                                                aGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTg                     >  1:3913697/1‑150 (MQ=255)
                                                                                                                                 gTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTgg                    <  1:436966/150‑1 (MQ=255)
                                                                                                                                      tCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGCTCtt               <  1:4071842/150‑1 (MQ=255)
                                                                                                                                      tCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGCTCtt               <  2:569713/150‑1 (MQ=255)
                                                                                                                                      tCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGCTCtt               <  2:1487070/150‑1 (MQ=255)
                                                                                                                                      tCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGCTCtt               >  2:4044014/1‑150 (MQ=255)
                                                                                                                                         tCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGCTCTTGcc            <  1:3142482/150‑1 (MQ=255)
                                                                                                                                            ccACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGCTCTTGCCACg         <  1:3745175/150‑1 (MQ=255)
                                                                                                                                            ccACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGCTCTTGCCACg         <  1:2564932/150‑1 (MQ=255)
                                                                                                                                                cGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGCTCTTGCCACGGTTc     <  1:799124/150‑1 (MQ=255)
                                                                                                                                                 gACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGCTCTTGCCACGGTTCa    <  1:5020865/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGCTCTTGCCACGGTTCAGC  >  CP009273/2302503‑2302799

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: