Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,302,650 C→A 19.8% 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,6500CA19.8% 884.0 / 169.0 310M57I (ATG→ATTalkBoxidative demethylase of N1‑methyladenine or N3‑methylcytosine DNA lesions
Reads supporting (aligned to +/- strand):  ref base C (125/123);  new base A (35/27);  total (160/150)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 311 total aligned reads displayed.

CGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGCTCTTGCCACGGTTCAG  >  CP009273/2302503‑2302798
                                                                                                                                                   |                                                                                                                                                    
cgccgTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACa                                                                                                                                                     >  1:1021911/1‑149 (MQ=255)
 gccgTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACAt                                                                                                                                                    <  1:1619063/149‑1 (MQ=255)
  ccgTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATg                                                                                                                                                   <  2:3392879/149‑1 (MQ=255)
    gTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGt                                                                                                                                                 <  1:1475402/149‑1 (MQ=255)
      aGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGtat                                                                                                                                               >  2:4281161/1‑149 (MQ=255)
      aGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGtat                                                                                                                                               <  1:3432267/149‑1 (MQ=255)
               ttGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCgg                                                                                                                                      >  2:3815154/1‑149 (MQ=255)
               ttGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCgg                                                                                                                                      >  2:2060447/1‑149 (MQ=255)
               ttGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCgg                                                                                                                                      >  2:634266/1‑149 (MQ=255)
               ttGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCgg                                                                                                                                      >  2:691496/1‑149 (MQ=255)
                tGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCggg                                                                                                                                     <  1:3068138/149‑1 (MQ=255)
                   cATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGt                                                                                                                                  >  1:2623599/1‑149 (MQ=255)
                   cATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGt                                                                                                                                  <  1:229594/149‑1 (MQ=255)
                    aTAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTg                                                                                                                                 >  1:3566062/1‑149 (MQ=255)
                     tAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGa                                                                                                                                <  1:3859145/149‑1 (MQ=255)
                       aaTTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGAcc                                                                                                                              <  1:4824694/149‑1 (MQ=255)
                        aTTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCa                                                                                                                             >  2:2736646/1‑149 (MQ=255)
                        aTTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCa                                                                                                                             >  2:2208305/1‑149 (MQ=255)
                        aTTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTCTATCCCCCGGGGGTGACCa                                                                                                                             >  2:3588597/1‑149 (MQ=255)
                          tATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATc                                                                                                                           >  1:1916803/1‑149 (MQ=255)
                            tGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTg                                                                                                                         >  1:3863273/1‑149 (MQ=255)
                            tGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTg                                                                                                                         >  2:609195/1‑149 (MQ=255)
                            tGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTg                                                                                                                         >  2:2235288/1‑149 (MQ=255)
                              aaaaCTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGc                                                                                                                       <  1:2748659/149‑1 (MQ=255)
                              aaaaCTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGc                                                                                                                       <  1:177313/149‑1 (MQ=255)
                                aaCTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGa                                                                                                                     >  2:4960554/1‑149 (MQ=255)
                                aaCTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGa                                                                                                                     >  1:1739819/1‑149 (MQ=255)
                                 aCTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGaa                                                                                                                    >  2:3959058/1‑149 (MQ=255)
                                  ctctGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGaaa                                                                                                                   >  2:5108884/1‑149 (MQ=255)
                                   tctGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAAc                                                                                                                  <  1:503438/149‑1 (MQ=255)
                                        ggCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGa                                                                                                             >  2:2637648/1‑149 (MQ=255)
                                         gCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGGGACCATCTGGCGAAACGGCGAc                                                                                                            <  1:3342899/149‑1 (MQ=255)
                                           aTGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGActg                                                                                                          >  1:4654666/1‑149 (MQ=255)
                                                gggCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGActggctgg                                                                                                     >  2:5053273/1‑149 (MQ=255)
                                                   ccACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCaa                                                                                                  >  1:2568909/1‑149 (MQ=255)
                                                    cACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAAc                                                                                                 >  1:587085/1‑149 (MQ=255)
                                                    cACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAAc                                                                                                 >  2:4860181/1‑149 (MQ=255)
                                                        gTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCa                                                                                             >  2:1261892/1‑149 (MQ=255)
                                                              ttgtttgCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAata                                                                                       <  1:657587/149‑1 (MQ=255)
                                                               tgtttgCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAatat                                                                                      <  1:3772683/149‑1 (MQ=255)
                                                               tgtttgCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTGGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAatat                                                                                      <  2:194335/149‑1 (MQ=255)
                                                                gtttgCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATc                                                                                     <  2:2884581/149‑1 (MQ=255)
                                                                gtttgCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATc                                                                                     <  2:4114497/149‑1 (MQ=255)
                                                                gtttgCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATc                                                                                     <  2:4816412/149‑1 (MQ=255)
                                                                    gCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCgg                                                                                 <  2:2337638/149‑1 (MQ=255)
                                                                    gCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCgg                                                                                 >  2:4797686/1‑149 (MQ=255)
                                                                     cGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGa                                                                                <  2:753391/149‑1 (MQ=255)
                                                                      ggATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGAt                                                                               <  2:523704/149‑1 (MQ=255)
                                                                      ggATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGAt                                                                               <  2:4298118/149‑1 (MQ=255)
                                                                       gATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATc                                                                              <  2:1730145/149‑1 (MQ=255)
                                                                         tCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAg                                                                            <  1:4015835/149‑1 (MQ=255)
                                                                          cAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGt                                                                           >  2:2914238/1‑149 (MQ=255)
                                                                          cAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGt                                                                           >  2:1555947/1‑149 (MQ=255)
                                                                          cAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGt                                                                           <  2:1285113/149‑1 (MQ=255)
                                                                           aaTGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGtt                                                                          <  2:3029670/149‑1 (MQ=255)
                                                                            aTGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTg                                                                         >  1:2685464/1‑149 (MQ=255)
                                                                            aTGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTg                                                                         >  1:1226871/1‑149 (MQ=255)
                                                                              gggCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCt                                                                       >  1:3618609/1‑149 (MQ=255)
                                                                               ggCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTc                                                                      >  1:2026302/1‑149 (MQ=255)
                                                                               ggCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTc                                                                      >  1:4246834/1‑149 (MQ=255)
                                                                                   aaTAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCa                                                                  >  2:3474444/1‑149 (MQ=255)
                                                                                    aTAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAg                                                                 <  1:2245771/149‑1 (MQ=255)
                                                                                       gagaTAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGt                                                              >  2:3020273/1‑149 (MQ=255)
                                                                                       gagaTAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGt                                                              <  2:198980/149‑1 (MQ=255)
                                                                                       gagaTAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGt                                                              >  2:2197070/1‑149 (MQ=255)
                                                                                       gagaTAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGt                                                              <  1:337862/149‑1 (MQ=255)
                                                                                       gagaTAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGt                                                              >  1:3096208/1‑149 (MQ=255)
                                                                                        agaTAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGtt                                                             <  2:4523350/149‑1 (MQ=255)
                                                                                             aCCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAg                                                        >  1:4387959/1‑149 (MQ=255)
                                                                                             aCCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAg                                                        <  2:1475078/149‑1 (MQ=255)
                                                                                               cTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCa                                                      <  1:2502703/149‑1 (MQ=255)
                                                                                               cTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCaa                                                     <  1:877172/149‑1 (MQ=255)
                                                                                               cTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTCTATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCAGTCATATCGCGGATCAGTTGCTCCGCCGCGTTAAAAGCa                                                      >  2:1812566/1‑149 (MQ=255)
                                                                                                    cGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACg                                                 >  2:439555/1‑149 (MQ=255)
                                                                                                    cGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACg                                                 >  2:1032128/1‑149 (MQ=255)
                                                                                                     gATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGc                                                >  1:1347526/1‑149 (MQ=255)
                                                                                                     gATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGc                                                <  1:1454579/149‑1 (MQ=255)
                                                                                                     gATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGc                                                >  1:3250543/1‑149 (MQ=255)
                                                                                                     gATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGc                                                >  2:825727/1‑149 (MQ=255)
                                                                                                             gtcCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAAt                                        >  2:3173924/1‑149 (MQ=255)
                                                                                                              tcCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAAtt                                       <  1:2205854/149‑1 (MQ=255)
                                                                                                                cAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTAc                                     <  2:4623315/149‑1 (MQ=255)
                                                                                                                 aGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTAcc                                    >  1:4104057/1‑149 (MQ=255)
                                                                                                                     cAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCAc                                >  2:1906689/1‑149 (MQ=255)
                                                                                                                       gATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCa                              >  2:2560675/1‑149 (MQ=255)
                                                                                                                       gATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCa                              >  2:1189873/1‑149 (MQ=255)
                                                                                                                         tGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCACAACGCCGTAAAATTACCGCACCAgc                            >  2:3755234/1‑149 (MQ=255)
                                                                                                                          gCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAgcc                           <  2:3296690/149‑1 (MQ=255)
                                                                                                                                aGTTGGTCATCGCCACCGACAGGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGt                     <  1:3279702/149‑1 (MQ=255)
                                                                                                                                 gTTGGTCATCGCCTCCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGAGACTGGCTGGCAAGTTCATTATTATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTg                    <  1:4895076/149‑1 (MQ=255)
                                                                                                                                   tGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGc                  >  2:1069324/1‑149 (MQ=255)
                                                                                                                                      tCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGctct               <  1:1218519/149‑1 (MQ=255)
                                                                                                                                      tCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGctct               <  1:634266/149‑1 (MQ=255)
                                                                                                                                       cATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGCTCtt              >  1:978184/1‑149 (MQ=255)
                                                                                                                                       cATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGCTCtt              >  2:122864/1‑149 (MQ=255)
                                                                                                                                       cATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGCTCtt              <  2:244163/149‑1 (MQ=255)
                                                                                                                                         tCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGCTCTTGc            >  2:401663/1‑149 (MQ=255)
                                                                                                                                               ccGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGCTCTTGCCACGGt      >  2:65484/1‑149 (MQ=255)
                                                                                                                                                  aCATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGCTCTTGCCACGGTTCa   >  1:2956446/1‑149 (MQ=255)
                                                                                                                                                   cATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGCTCTTGCCACGGTTCAg  >  1:1402424/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGCTCTTGCCACGGTTCAG  >  CP009273/2302503‑2302798

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: