Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,302,650 C→A 18.9% 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,6500CA18.9% 940.3 / 155.1 323M57I (ATG→ATTalkBoxidative demethylase of N1‑methyladenine or N3‑methylcytosine DNA lesions
Reads supporting (aligned to +/- strand):  ref base C (143/118);  new base A (27/35);  total (170/153)
Fisher's exact test for biased strand distribution p-value = 1.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 326 total aligned reads displayed.

GCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGCTCTTGCCACGGTTCA  >  CP009273/2302501‑2302797
                                                                                                                                                     |                                                                                                                                                   
gccgccgTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCAACGCCACCGAc                                                                                                                                                     <  1:1446606/150‑1 (MQ=255)
gccgccgTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCAACGCCACCGAc                                                                                                                                                     <  1:1912291/150‑1 (MQ=255)
  cgccgTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACAt                                                                                                                                                   >  1:576711/1‑150 (MQ=255)
   gccgTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATg                                                                                                                                                  >  1:3609402/1‑150 (MQ=255)
   gccgTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATg                                                                                                                                                  <  2:4321393/150‑1 (MQ=255)
        aGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGtata                                                                                                                                             >  1:805599/1‑150 (MQ=255)
        aGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGtata                                                                                                                                             >  1:3236500/1‑150 (MQ=255)
        aGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGtata                                                                                                                                             >  1:4247992/1‑150 (MQ=255)
           cGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATATcc                                                                                                                                          <  1:1248223/150‑1 (MQ=255)
             cACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATcccc                                                                                                                                        <  1:3219941/150‑1 (MQ=255)
               cGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCg                                                                                                                                      <  1:3886354/150‑1 (MQ=255)
                gTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCgg                                                                                                                                     <  1:976408/150‑1 (MQ=255)
                 ttGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCggg                                                                                                                                    <  1:1146778/150‑1 (MQ=255)
                   gACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCggggg                                                                                                                                  <  2:4427539/150‑1 (MQ=255)
                    aCATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGt                                                                                                                                 <  1:2439300/150‑1 (MQ=255)
                     cATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTg                                                                                                                                >  2:3444792/1‑150 (MQ=255)
                      aTAAATTATGAAAACTCTGTGTCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGa                                                                                                                               >  1:1589859/1‑150 (MQ=255)
                       tAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGAc                                                                                                                              >  2:2498552/1‑150 (MQ=255)
                        aaaTTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGAcc                                                                                                                             >  2:4657805/1‑150 (MQ=255)
                          aTTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCAt                                                                                                                           >  2:391372/1‑150 (MQ=255)
                          aTTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCAt                                                                                                                           >  2:3066829/1‑150 (MQ=255)
                             aTGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTg                                                                                                                        <  2:334954/150‑1 (MQ=255)
                              tGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTgg                                                                                                                       <  1:207023/150‑1 (MQ=255)
                                   aCTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGaaa                                                                                                                  <  1:2822500/150‑1 (MQ=255)
                                     tctGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACg                                                                                                                <  1:329655/150‑1 (MQ=255)
                                       tgtgGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGc                                                                                                              >  2:4252925/1‑150 (MQ=255)
                                       tgtgGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGc                                                                                                              <  1:4701733/150‑1 (MQ=255)
                                       tgtgGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGc                                                                                                              >  2:4095692/1‑150 (MQ=255)
                                          ggCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGAc                                                                                                           >  2:523747/1‑150 (MQ=255)
                                            cATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGActg                                                                                                         >  2:488342/1‑150 (MQ=255)
                                            cATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGActg                                                                                                         >  2:3363200/1‑150 (MQ=255)
                                            cATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGActg                                                                                                         >  1:2484213/1‑150 (MQ=255)
                                              tGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGActggc                                                                                                       <  1:27418/150‑1 (MQ=255)
                                                 cGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGActggctgg                                                                                                    <  1:3850782/150‑1 (MQ=255)
                                                  gggCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACtggctggc                                                                                                   >  2:2997225/1‑150 (MQ=255)
                                                  gggCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACtggctggc                                                                                                   >  2:563665/1‑150 (MQ=255)
                                                  gggCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACtggctggc                                                                                                   >  1:311633/1‑150 (MQ=255)
                                                   ggCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCa                                                                                                  <  2:1591638/150‑1 (MQ=255)
                                                    gCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTCTATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCaa                                                                                                 >  2:4792336/1‑150 (MQ=255)
                                                     ccACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAAc                                                                                                >  1:4727947/1‑150 (MQ=255)
                                                        cGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTc                                                                                             >  1:4442553/1‑150 (MQ=255)
                                                         ggTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCa                                                                                            <  1:1195854/150‑1 (MQ=255)
                                                         ggTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCa                                                                                            >  2:4539283/1‑150 (MQ=255)
                                                         ggTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCa                                                                                            >  1:4783752/1‑150 (MQ=255)
                                                         ggTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCa                                                                                            <  1:4095683/150‑1 (MQ=255)
                                                         ggTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCa                                                                                            >  1:2797118/1‑150 (MQ=255)
                                                           tttATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCAtt                                                                                          <  1:4191265/150‑1 (MQ=255)
                                                           tttATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCAtt                                                                                          >  1:233959/1‑150 (MQ=255)
                                                            ttATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTa                                                                                         >  2:2402474/1‑150 (MQ=255)
                                                            ttATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTa                                                                                         <  1:2048740/150‑1 (MQ=255)
                                                              aTTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAat                                                                                       <  1:4824701/150‑1 (MQ=255)
                                                              aTTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAat                                                                                       >  2:3234054/1‑150 (MQ=255)
                                                               tttgtttgCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAata                                                                                      <  2:1454915/150‑1 (MQ=255)
                                                                   tttgCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATcgc                                                                                  >  1:897407/1‑150 (MQ=255)
                                                                       cGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGAt                                                                              >  2:2795703/1‑150 (MQ=255)
                                                                        ggATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACATCATTAATATCGCGGATc                                                                             >  2:1961982/1‑150 (MQ=255)
                                                                         gATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCa                                                                            >  2:3695145/1‑150 (MQ=255)
                                                                         gATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGCTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCa                                                                            <  2:2889792/150‑1 (MQ=255)
                                                                          aTCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAg                                                                           <  1:1747783/150‑1 (MQ=255)
                                                                            cAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGtt                                                                         <  2:1351761/150‑1 (MQ=255)
                                                                            cAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGtt                                                                         >  1:2203019/1‑150 (MQ=255)
                                                                            cAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGtt                                                                         >  2:2759169/1‑150 (MQ=255)
                                                                            cAATGGGCGAATAGAGATAACCGTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGtt                                                                         <  2:3831945/150‑1 (MQ=255)
                                                                               tGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCt                                                                      >  1:3455340/1‑150 (MQ=255)
                                                                                  gCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCg                                                                   <  2:4256730/150‑1 (MQ=255)
                                                                                    gAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCa                                                                 <  2:1230392/150‑1 (MQ=255)
                                                                                     aaTAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAg                                                                >  2:3989653/1‑150 (MQ=255)
                                                                                       tAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCACCGTCATTAATATCGCGGATCAGTTGCTCCGCAGCg                                                              >  1:1856266/1‑150 (MQ=255)
                                                                                         gagaTAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGtt                                                            >  2:1110832/1‑150 (MQ=255)
                                                                                         gagaTAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGtt                                                            >  1:3444058/1‑150 (MQ=255)
                                                                                                  ttGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCaaa                                                   <  1:1001258/150‑1 (MQ=255)
                                                                                                      cGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGc                                               >  2:1599463/1‑150 (MQ=255)
                                                                                                       gATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGcc                                              <  1:2656793/150‑1 (MQ=255)
                                                                                                       gATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGcc                                              >  1:3394333/1‑150 (MQ=255)
                                                                                                       gATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGcc                                              >  1:4156862/1‑150 (MQ=255)
                                                                                                        aTGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCg                                             <  2:3147718/150‑1 (MQ=255)
                                                                                                           ggTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTaa                                          <  1:3991163/150‑1 (MQ=255)
                                                                                                           ggTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTaa                                          >  2:3891036/1‑150 (MQ=255)
                                                                                                              cgtcCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAAt                                       >  1:3471989/1‑150 (MQ=255)
                                                                                                               gtcCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAAtt                                      >  1:2675434/1‑150 (MQ=255)
                                                                                                               gtcCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAAtt                                      <  2:3807298/150‑1 (MQ=255)
                                                                                                                tcCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTa                                     >  2:2724093/1‑150 (MQ=255)
                                                                                                                tcCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTa                                     >  2:290394/1‑150 (MQ=255)
                                                                                                                  cAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTAcc                                   >  1:1650934/1‑150 (MQ=255)
                                                                                                                  cAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTAcc                                   >  2:3546103/1‑150 (MQ=255)
                                                                                                                   aGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCg                                  >  1:3675139/1‑150 (MQ=255)
                                                                                                                      ccAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCAc                               >  1:4365892/1‑150 (MQ=255)
                                                                                                                        aGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCa                             <  2:2150272/150‑1 (MQ=255)
                                                                                                                                  aGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAACACGCCGTAAAATTACCGCACCAGCCGCCAGTg                   <  2:2615054/150‑1 (MQ=255)
                                                                                                                                  aGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTg                   <  1:853216/150‑1 (MQ=255)
                                                                                                                                  aGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTg                   <  1:2707881/150‑1 (MQ=255)
                                                                                                                                   gTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTgg                  <  1:2913235/150‑1 (MQ=255)
                                                                                                                                   gTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTgg                  <  2:1882800/150‑1 (MQ=255)
                                                                                                                                      ggTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGctc               >  2:1759829/1‑150 (MQ=255)
                                                                                                                                         cATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTAATTAATATCGCGGATCAGTTGCTCCGCATCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGCTCTTg            <  2:1172100/150‑1 (MQ=255)
                                                                                                                                         cATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGCTCTTg            >  2:1818854/1‑150 (MQ=255)
                                                                                                                                          aTCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGCTCTTGc           >  2:1256131/1‑150 (MQ=255)
                                                                                                                                           tCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGCTCTTGcc          <  1:4541158/150‑1 (MQ=255)
                                                                                                                                               cACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGCTCTTGCCACgg      <  1:713239/150‑1 (MQ=255)
                                                                                                                                                aCCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGCTCTTGCCACGGt     >  1:888681/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGCTCTTGCCACGGTTCA  >  CP009273/2302501‑2302797

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: