Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,302,650 C→A 6.5% 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,6500CA 6.5% 1053.6 / 36.1 276M57I (ATG→ATTalkBoxidative demethylase of N1‑methyladenine or N3‑methylcytosine DNA lesions
Reads supporting (aligned to +/- strand):  ref base C (122/136);  new base A (10/8);  total (132/144)
Fisher's exact test for biased strand distribution p-value = 6.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 276 total aligned reads displayed.

CCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGCTCTTGCCACGGTTC  >  CP009273/2302502‑2302796
                                                                                                                                                    |                                                                                                                                                  
ccgccgTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACa                                                                                                                                                   <  1:5321406/150‑1 (MQ=255)
   ccgTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATgg                                                                                                                                                <  1:2259538/150‑1 (MQ=255)
     ggaGCCGCACGTTGACATCAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGta                                                                                                                                              <  2:2901169/148‑1 (MQ=255)
     gTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGta                                                                                                                                              <  2:133195/150‑1 (MQ=255)
     gTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGta                                                                                                                                              <  1:1672721/150‑1 (MQ=255)
         ccGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATc                                                                                                                                          >  1:1110863/1‑150 (MQ=255)
          cGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATcc                                                                                                                                         <  1:3354627/150‑1 (MQ=255)
           gCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATccc                                                                                                                                        <  1:1898589/150‑1 (MQ=255)
            cACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATcccc                                                                                                                                       <  1:164458/150‑1 (MQ=255)
                  gACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGAGGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCggggg                                                                                                                                 >  2:4064208/1‑150 (MQ=255)
                   aCATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGt                                                                                                                                <  1:1337684/150‑1 (MQ=255)
                    cATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGAGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTg                                                                                                                               <  2:2540125/150‑1 (MQ=255)
                     aTAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGTTGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGa                                                                                                                              <  1:5405825/150‑1 (MQ=255)
                     aTAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGa                                                                                                                              <  1:837861/150‑1 (MQ=255)
                      tAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGAc                                                                                                                             <  2:741852/150‑1 (MQ=255)
                       aaaTTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGAcc                                                                                                                            >  1:4275732/1‑150 (MQ=255)
                        aaTTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCa                                                                                                                           >  1:2948278/1‑150 (MQ=255)
                           tATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCt                                                                                                                        >  2:1890021/1‑150 (MQ=255)
                            aTGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTg                                                                                                                       >  2:1213600/1‑150 (MQ=255)
                                agaCTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGa                                                                                                                   >  1:3982206/3‑150 (MQ=255)
                                aaaCTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGa                                                                                                                   >  1:2652718/1‑150 (MQ=255)
                                  aCTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGaaa                                                                                                                 <  2:5332544/150‑1 (MQ=255)
                                  aCTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGaaa                                                                                                                 >  1:416680/1‑150 (MQ=255)
                                   ctctGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAAc                                                                                                                <  2:396138/150‑1 (MQ=255)
                                   ctcgGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAAc                                                                                                                <  2:2590837/150‑1 (MQ=255)
                                      tgtgGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGACTAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGc                                                                                                             >  2:4215762/1‑150 (MQ=255)
                                      tgtgGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGc                                                                                                             <  1:1259275/150‑1 (MQ=255)
                                          gCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGAct                                                                                                         <  1:4473433/150‑1 (MQ=255)
                                              ggCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGActggct                                                                                                     <  1:5136787/150‑1 (MQ=255)
                                                cGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGActggctgg                                                                                                   <  1:1745933/150‑1 (MQ=255)
                                                cGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGActggctgg                                                                                                   <  1:1745638/150‑1 (MQ=255)
                                                  ggCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCa                                                                                                 >  2:3519619/1‑150 (MQ=255)
                                                     cACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACg                                                                                              >  2:870969/1‑150 (MQ=255)
                                                       cGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTc                                                                                            >  1:2072274/1‑150 (MQ=255)
                                                         gTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAAc                                                                                                                >  1:2494625/1‑128 (MQ=255)
                                                           ttATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTa                                                                                        <  2:1084440/150‑1 (MQ=255)
                                                            tATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTaa                                                                                       <  2:1571751/150‑1 (MQ=255)
                                                            tATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTaa                                                                                       <  2:2062632/150‑1 (MQ=255)
                                                              tttgtttgCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCAGCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAata                                                                                     >  2:1758832/1‑150 (MQ=255)
                                                               ttgtttgCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAatat                                                                                    >  2:3257404/1‑150 (MQ=255)
                                                                 gtttgCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATcg                                                                                  <  1:959678/150‑1 (MQ=255)
                                                                 gtttgCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATcg                                                                                  <  1:2582401/150‑1 (MQ=255)
                                                                  tttgCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATcgc                                                                                 <  1:336632/150‑1 (MQ=255)
                                                                      cGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGAt                                                                             >  1:2081922/1‑150 (MQ=255)
                                                                        gATCAATGGGCGAATAGAGATAACCTTGTCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCa                                                                           <  1:4699913/150‑1 (MQ=255)
                                                                        gATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCa                                                                           >  2:2277272/1‑150 (MQ=255)
                                                                        gATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCa                                                                           <  1:2307090/150‑1 (MQ=255)
                                                                          tCATTGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGt                                                                         <  2:924854/150‑1 (MQ=255)
                                                                           cAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGtt                                                                        <  2:3643592/150‑1 (MQ=255)
                                                                           cAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGtt                                                                        >  1:3279338/1‑150 (MQ=255)
                                                                            aaTGGGCGAATAGAGATAACCTTTCCTAAGGGTCGACCAGCCCAGATGCCCACAGTTGGTCAACGCCACAGACAAGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTATAATCGCGGATCAGTTg                                                                       >  1:3549349/1‑150 (MQ=255)
                                                                             aTGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGc                                                                      <  1:2785370/150‑1 (MQ=255)
                                                                             aTGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCAGTAATATCGCGGATCAGTTGc                                                                      >  2:4977837/1‑150 (MQ=255)
                                                                               gggCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTc                                                                    <  2:4755617/150‑1 (MQ=255)
                                                                               gggCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTc                                                                    >  2:1145279/1‑150 (MQ=255)
                                                                               gggCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTc                                                                    >  2:2977795/1‑150 (MQ=255)
                                                                                ggCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTcc                                                                   >  2:4430969/1‑150 (MQ=255)
                                                                                ggCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTcc                                                                   <  2:416585/150‑1 (MQ=255)
                                                                                 gCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCg                                                                  <  2:1190018/150‑1 (MQ=255)
                                                                                   gAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCa                                                                >  1:1102/1‑150 (MQ=255)
                                                                                   gAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCa                                                                >  1:3140768/1‑150 (MQ=255)
                                                                                    aaTAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAg                                                               <  1:780141/150‑1 (MQ=255)
                                                                                    aaTAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAg                                                               >  2:3040050/1‑150 (MQ=255)
                                                                                     aTAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGc                                                              >  1:3958561/1‑150 (MQ=255)
                                                                                        gagaTAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGtt                                                           <  2:3854162/150‑1 (MQ=255)
                                                                                        gagaTAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGtt                                                           <  2:4939957/150‑1 (MQ=255)
                                                                                           aTAACCTTGCCGATGGGTCGTCCAGCCCAGATGACCACAGTTGGTCATCGCAAACGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTaaa                                                        >  2:2040108/1‑150 (MQ=255)
                                                                                            tAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTaaaa                                                       <  2:1806231/150‑1 (MQ=255)
                                                                                                 ttGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCaaa                                                  >  2:4140549/1‑150 (MQ=255)
                                                                                                     cGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAATGCAAAACGc                                              >  1:275394/1‑150 (MQ=255)
                                                                                                     cGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGc                                              >  2:2626209/1‑150 (MQ=255)
                                                                                                         gggTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTa                                          <  1:3433593/150‑1 (MQ=255)
                                                                                                           gtcgtcCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTaaa                                        <  1:5366823/150‑1 (MQ=255)
                                                                                                           gtcgtcCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTaaa                                        <  1:3680208/150‑1 (MQ=255)
                                                                                                              gtcCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAAtt                                     <  1:4259319/150‑1 (MQ=255)
                                                                                                               tcCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTa                                    >  1:856827/1‑150 (MQ=255)
                                                                                                                 cAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTCAAATTAcc                                  >  2:2671611/1‑150 (MQ=255)
                                                                                                                 cAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTAcc                                  <  2:4609080/150‑1 (MQ=255)
                                                                                                                  aGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCg                                 <  1:525040/150‑1 (MQ=255)
                                                                                                                  aGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCg                                 >  1:4571208/1‑150 (MQ=255)
                                                                                                                   gCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGc                                <  1:2358160/150‑1 (MQ=255)
                                                                                                                     ccAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCAc                              <  2:3823530/150‑1 (MQ=255)
                                                                                                                       aGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCa                            >  2:2472567/1‑150 (MQ=255)
                                                                                                                       aGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCa                            >  1:3897169/1‑150 (MQ=255)
                                                                                                                          tGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAgcc                         >  2:1703374/1‑150 (MQ=255)
                                                                                                                              cacaGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCa                     <  1:4064208/150‑1 (MQ=255)
                                                                                                                               acaGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAg                    >  1:1514474/1‑150 (MQ=255)
                                                                                                                               acaGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCACCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAg                    >  1:1514147/1‑150 (MQ=255)
                                                                                                                               acaGTAGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAg                    <  2:94817/150‑1 (MQ=255)
                                                                                                                                   ttGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGc                <  2:2357397/150‑1 (MQ=255)
                                                                                                                                      gTCATCGCCACCGACATGGTATATCCACCGGGGGTGACCATCTGGCGAAACGGCGTCTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGACAGTGGctct             <  2:4566654/150‑1 (MQ=255)
                                                                                                                                       tCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGCTCtt            <  2:3706181/150‑1 (MQ=255)
                                                                                                                                       tCATCGCCACCGAAATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGCTCtt            >  2:4367452/1‑150 (MQ=255)
                                                                                                                                         aTCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGCTCTTGc          >  1:4346367/1‑150 (MQ=255)
                                                                                                                                         aTCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGCTCTTGc          >  2:2265697/1‑150 (MQ=255)
                                                                                                                                             ccACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGCTCTTGCCACg      >  2:1516813/1‑150 (MQ=255)
                                                                                                                                                ccGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGCTCTTGCCACGGtt   <  1:4106832/150‑1 (MQ=255)
                                                                                                                                                 cGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGCTCTTGCCACGGTTc  >  1:5261641/1‑150 (MQ=255)
                                                                                                                                                 cGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGCTCTTGCCACGGTTc  <  1:879145/150‑1 (MQ=255)
                                                                                                                                                 cGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGCTCTTGCCACGGTTc  >  1:2551870/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGACATGGTATATCCCCCGGGGGTGACCATCTGGCGAAACGGCGACTGGCTGGCAACGTCATTAATATCGCGGATCAGTTGCTCCGCAGCGTTAAAAGCAAAACGCCGTAAAATTACCGCACCAGCCGCCAGTGGCTCTTGCCACGGTTC  >  CP009273/2302502‑2302796

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: