Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,163,117 T→A 66.7% I36I (ATT→ATA birA → bifunctional biotin‑[acetylCoA carboxylase] holoenzyme synthetase/ DNA‑binding transcriptional repressor, bio‑5'‑AMP‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,163,1170TA66.7% 293.9 / 221.5 222I36I (ATT→ATAbirAbifunctional biotin‑[acetylCoA carboxylase] holoenzyme synthetase/ DNA‑binding transcriptional repressor, bio‑5'‑AMP‑binding
Reads supporting (aligned to +/- strand):  ref base T (34/40);  new base A (68/80);  total (102/120)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 232 total aligned reads displayed.

TGGTGCATCAGGTGAAGTGAGCGCAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGTAGTGTAGCCGTGCTG  >  CP009273/4162971‑4163261
                                                                                                                                                  |                                                                                                                                                
tGGTGCATCAGGTGAAGTGAGCGCAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATta                                                                                                                                                 >  2:689242/1‑148 (MQ=255)
tGGTGCATCAGGTGAAGTGAGCGCAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATta                                                                                                                                                 >  2:1751184/1‑148 (MQ=255)
  gTGCATCATGTGAAGTGAGCGCAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATaaat                                                                                                                                               <  2:555659/148‑1 (MQ=255)
  gTGCATCAGGTGAAGTGAGCGCAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATaaat                                                                                                                                               <  2:2346412/148‑1 (MQ=255)
        cAGGTGAAGTGAGCGCAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAacac                                                                                                                                         <  2:1938203/148‑1 (MQ=255)
         aGGTGAAGTGAGCGCAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAAcaca                                                                                                                                        >  1:2969967/1‑148 (MQ=255)
          ggTGAAGTGAGCGCAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAACTTGATTGCCCAGATAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAAcacag                                                                                                                                       >  2:2096993/1‑147 (MQ=255)
          ggTGAAGTGAGCGCAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACAt                                                                                                                                       <  1:1101100/148‑1 (MQ=255)
          ggTGAAGTGAGCGCAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACAt                                                                                                                                       <  1:2782117/148‑1 (MQ=255)
          ggTGAAGTGAGCGCAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACAt                                                                                                                                       >  2:1296901/1‑148 (MQ=255)
            tGAAGTGAGCGCAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTc                                                                                                                                     <  1:483707/148‑1 (MQ=255)
                gTGAGCGCAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGac                                                                                                                                 <  2:1182397/148‑1 (MQ=255)
                gTGAGCGCAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGac                                                                                                                                 >  1:305518/1‑148 (MQ=255)
                gTGAGCGCAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGac                                                                                                                                 >  2:2385083/1‑148 (MQ=255)
                 tGAGCGCAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGaca                                                                                                                                >  1:3305423/1‑148 (MQ=255)
                 tGAGCGCAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGaca                                                                                                                                <  2:2006727/148‑1 (MQ=255)
                 tGAGCGCAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGaca                                                                                                                                <  1:205415/148‑1 (MQ=255)
                  gAGCGCAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGacac                                                                                                                               >  2:806467/1‑148 (MQ=255)
                     cgcAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGc                                                                                                                            <  2:2407907/148‑1 (MQ=255)
                      gcAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCg                                                                                                                           <  2:3012617/148‑1 (MQ=255)
                       cAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGt                                                                                                                          >  2:755379/1‑148 (MQ=255)
                         gTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGa                                                                                                                        <  1:1904339/148‑1 (MQ=255)
                          tGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGAc                                                                                                                       >  1:110875/1‑148 (MQ=255)
                            gagaCAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTg                                                                                                                     >  1:831912/1‑148 (MQ=255)
                            gagaCAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTg                                                                                                                     >  2:873395/1‑148 (MQ=255)
                            gagaCAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCACTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTg                                                                                                                     >  2:1239604/1‑148 (MQ=255)
                             agaCAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTgg                                                                                                                    >  1:1885567/1‑148 (MQ=255)
                             agaCAATTTCATGAAGGATAACACAGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTgg                                                                                                                    >  2:1637649/1‑148 (MQ=255)
                              gaCAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTggg                                                                                                                   <  2:2368564/148‑1 (MQ=255)
                              gaCAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTggg                                                                                                                   >  1:1694242/1‑148 (MQ=255)
                               aCAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTgggg                                                                                                                  >  2:3109663/1‑148 (MQ=255)
                                cAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGc                                                                                                                 <  1:1309007/148‑1 (MQ=255)
                                 aaTTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCg                                                                                                                <  1:806845/148‑1 (MQ=255)
                                 aaTTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCg                                                                                                                <  1:1729764/148‑1 (MQ=255)
                                  aTTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGt                                                                                                               <  2:2793961/148‑1 (MQ=255)
                                  aTTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGt                                                                                                               <  2:190775/148‑1 (MQ=255)
                                   tttCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGtt                                                                                                              >  1:3317350/1‑148 (MQ=255)
                                    ttCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTg                                                                                                             <  2:2798097/148‑1 (MQ=255)
                                    ttCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTg                                                                                                             <  1:2812090/148‑1 (MQ=255)
                                        tGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGt                                                                                                         >  2:670149/1‑148 (MQ=255)
                                          aaGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCt                                                                                                       <  1:2096994/148‑1 (MQ=255)
                                           aGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCtt                                                                                                      <  1:3076035/148‑1 (MQ=255)
                                               tAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTAcc                                                                                                  >  2:2937019/1‑148 (MQ=255)
                                                 acacCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGt                                                                                                <  2:2438064/148‑1 (MQ=255)
                                                  ccccGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGtt                                                                                               <  2:2174763/146‑1 (MQ=255)
                                                  cacCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGtt                                                                                               >  2:729131/1‑148 (MQ=255)
                                                       tGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCggg                                                                                          <  2:3390682/148‑1 (MQ=255)
                                                         ccACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTa                                                                                        >  2:1608699/1‑148 (MQ=255)
                                                          cACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTaa                                                                                       >  2:2512844/1‑148 (MQ=255)
                                                          cACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTaa                                                                                       >  2:2508502/1‑148 (MQ=255)
                                                           aCTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTaaa                                                                                      <  1:266999/148‑1 (MQ=255)
                                                            cTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAg                                                                                     >  1:1014463/1‑148 (MQ=255)
                                                              gAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGa                                                                                   >  2:484211/1‑148 (MQ=255)
                                                               aaaTTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGAt                                                                                  <  2:2101190/148‑1 (MQ=255)
                                                                   tgattgCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAg                                                                              <  1:2561009/148‑1 (MQ=255)
                                                                   tgattgCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAg                                                                              >  1:1392226/1‑148 (MQ=255)
                                                                     attgCCCCGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAgcc                                                                            >  2:1374019/1‑148 (MQ=255)
                                                                        gCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAgcctgc                                                                         <  2:1443164/148‑1 (MQ=255)
                                                                             gTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAg                                                                    <  2:831912/148‑1 (MQ=255)
                                                                              ttAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGc                                                                   <  2:2905697/148‑1 (MQ=255)
                                                                               tAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGcc                                                                  >  1:113558/1‑148 (MQ=255)
                                                                               tAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGcc                                                                  >  2:1416139/1‑148 (MQ=255)
                                                                                  cGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTAt                                                               <  2:523803/148‑1 (MQ=255)
                                                                                    aaCGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATcc                                                             >  1:562945/1‑148 (MQ=255)
                                                                                      cGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAg                                                           <  1:965852/148‑1 (MQ=255)
                                                                                      cGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAg                                                           <  2:2688164/148‑1 (MQ=255)
                                                                                         tGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTa                                                        >  2:3295719/1‑148 (MQ=255)
                                                                                            aTTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACtt                                                     >  1:3190785/1‑148 (MQ=255)
                                                                                               tCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAAt                                                  >  2:1736072/1‑148 (MQ=255)
                                                                                                 aCTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGc                                                <  1:755379/148‑1 (MQ=255)
                                                                                                   tctGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTa                                              <  1:2001988/148‑1 (MQ=255)
                                                                                                       gCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACa                                          <  1:3528632/148‑1 (MQ=255)
                                                                                                       gCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACa                                          <  1:1768396/148‑1 (MQ=255)
                                                                                                       gCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACa                                          <  2:752330/148‑1 (MQ=255)
                                                                                                           gCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACAGCGTGACTGGGGCGTTGATTTCTTTACCGTTCCGGGTAAATGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGata                                      >  2:1822482/1‑148 (MQ=255)
                                                                                                             aGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATAtt                                    <  2:1014463/148‑1 (MQ=255)
                                                                                                                 gggTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGt                                <  1:1696823/148‑1 (MQ=255)
                                                                                                                            tGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATgg                     <  1:1298294/148‑1 (MQ=255)
                                                                                                                            tGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATgg                     <  1:1637649/148‑1 (MQ=255)
                                                                                                                             gggAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGc                    <  2:3305530/148‑1 (MQ=255)
                                                                                                                              ggAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCg                   >  2:2345458/1‑148 (MQ=255)
                                                                                                                              ggAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCg                   <  1:3402957/148‑1 (MQ=255)
                                                                                                                              ggAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCg                   <  2:1080555/148‑1 (MQ=255)
                                                                                                                              ggAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCg                   >  1:1971031/1‑148 (MQ=255)
                                                                                                                                aaTGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGt                 <  1:2303489/148‑1 (MQ=255)
                                                                                                                                aaTGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGt                 <  1:2508401/148‑1 (MQ=255)
                                                                                                                                aaTGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGt                 <  1:2512844/148‑1 (MQ=255)
                                                                                                                                 aTGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGTa                <  1:1750362/148‑1 (MQ=255)
                                                                                                                                 aTGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGTa                >  1:2865076/1‑148 (MQ=255)
                                                                                                                                  tGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGTAg               >  1:38394/1‑148 (MQ=255)
                                                                                                                                   gAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTACTGCTAAACAGATATTGGGTCAGCTGGATGGCGGTAgt              <  2:3119636/148‑1 (MQ=255)
                                                                                                                                      ccGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGTAGTGTa           >  2:3496236/1‑148 (MQ=255)
                                                                                                                                      ccGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGTAGTGTa           >  1:1838885/1‑148 (MQ=255)
                                                                                                                                         ggcggcTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGTAGTGTAGcc        >  1:3018998/1‑148 (MQ=255)
                                                                                                                                          gcggcTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGATCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGTAGTGTAGCCg       <  2:1768174/148‑1 (MQ=255)
                                                                                                                                          gcggcTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGTAGTGTAGCCg       >  1:690679/1‑148 (MQ=255)
                                                                                                                                            ggcTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGTAGTGTAGCCGtg     <  1:3436120/148‑1 (MQ=255)
                                                                                                                                             gcTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGTAGTGTAGCCGtgc    <  2:3140747/148‑1 (MQ=255)
                                                                                                                                              cTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGTAGTGTAGCCGtgct   >  2:2804713/1‑148 (MQ=255)
                                                                                                                                               tATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGTAGTGTAGCCGtgctg  >  2:2623794/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
TGGTGCATCAGGTGAAGTGAGCGCAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGTAGTGTAGCCGTGCTG  >  CP009273/4162971‑4163261

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: