Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,163,117 T→A 37.8% 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,1170TA37.8% 208.4 / 269.1 222I36I (ATT→ATAbirAbifunctional biotin‑[acetylCoA carboxylase] holoenzyme synthetase/ DNA‑binding transcriptional repressor, bio‑5'‑AMP‑binding
Reads supporting (aligned to +/- strand):  ref base T (74/64);  new base A (49/35);  total (123/99)
Fisher's exact test for biased strand distribution p-value = 5.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

Only 100 of 232 total aligned reads displayed.

ATTGGTGCATCAGGTGAAGTGAGCGCAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGTAGTGTAGCCGTGCTGCCAG  >  CP009273/4162969‑4163265
                                                                                                                                                    |                                                                                                                                                    
  tGGTGCATCAGGTGAAGTGAGCGCAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATtaa                                                                                                                                                    <  2:2938860/149‑1 (MQ=255)
     tGCATCAGGTGAAGTGAGCGCAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATtaataa                                                                                                                                                 <  2:3237031/149‑1 (MQ=255)
      gCATCAGGTGAAGTGAGCGCAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATaaa                                                                                                                                                <  1:3151029/149‑1 (MQ=255)
         tCAGGTGAAGTGAGCGCAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAacac                                                                                                                                             <  1:243440/149‑1 (MQ=255)
          cAGGTGAAGTGAGCGCAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAAcaca                                                                                                                                            <  2:3324295/149‑1 (MQ=255)
            ggTGAAGTGAGCGCAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACAtt                                                                                                                                          <  2:1380183/149‑1 (MQ=255)
                  gTGAGCGCAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGaca                                                                                                                                    <  1:1740282/149‑1 (MQ=255)
                   tGAGCGCAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGacac                                                                                                                                   >  2:3123359/1‑149 (MQ=255)
                     aGCGCAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTg                                                                                                                                 <  1:1785560/149‑1 (MQ=255)
                          aGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGa                                                                                                                            <  1:157001/149‑1 (MQ=255)
                                 aCAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGc                                                                                                                     >  1:1248941/1‑149 (MQ=255)
                                   aaTTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGt                                                                                                                   >  2:182419/1‑149 (MQ=255)
                                      ttCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGTGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGa                                                                                                                <  2:103645/149‑1 (MQ=255)
                                      ttCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGa                                                                                                                <  1:252816/149‑1 (MQ=255)
                                         aTGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGt                                                                                                             >  2:3621755/1‑149 (MQ=255)
                                         aTGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGt                                                                                                             >  2:3753069/1‑149 (MQ=255)
                                          tGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTc                                                                                                            <  2:3977362/149‑1 (MQ=255)
                                          tGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTc                                                                                                            >  2:1638488/1‑149 (MQ=255)
                                          tGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTc                                                                                                            >  1:3654427/1‑149 (MQ=255)
                                           gAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCt                                                                                                           >  2:2299413/1‑149 (MQ=255)
                                           gAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCt                                                                                                           <  2:3694868/149‑1 (MQ=255)
                                            aaGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCtt                                                                                                          >  1:814267/1‑149 (MQ=255)
                                            aaGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCtt                                                                                                          >  1:2871929/1‑149 (MQ=255)
                                              ggATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTa                                                                                                        >  1:1655027/1‑149 (MQ=255)
                                                  aaCACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGt                                                                                                    >  1:4089805/1‑149 (MQ=255)
                                                   acacCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGtt                                                                                                   >  1:234296/1‑149 (MQ=255)
                                                   acacCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGtt                                                                                                   >  1:1460513/1‑149 (MQ=255)
                                                      ccGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCg                                                                                                >  1:548884/1‑149 (MQ=255)
                                                      ccGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCg                                                                                                >  2:2483143/1‑149 (MQ=255)
                                                      ccGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCg                                                                                                >  1:1171940/1‑149 (MQ=255)
                                                        gTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCggg                                                                                              <  1:1162825/149‑1 (MQ=255)
                                                                gAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGAt                                                                                      <  2:1098202/149‑1 (MQ=255)
                                                                 aaaTTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATa                                                                                     <  1:2352483/149‑1 (MQ=255)
                                                                  aaTTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATAc                                                                                    <  2:3692245/149‑1 (MQ=255)
                                                                   attgattgCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACa                                                                                   >  2:182160/1‑149 (MQ=255)
                                                                   attgattgCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACa                                                                                   >  1:2694166/1‑149 (MQ=255)
                                                                      gattgCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAgcc                                                                                >  2:2617279/1‑149 (MQ=255)
                                                                        ttgCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCCGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAgcctg                                                                              >  2:2414228/1‑149 (MQ=255)
                                                                        ttgCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAgcctg                                                                              >  2:2691738/1‑149 (MQ=255)
                                                                          gCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAgcctgcc                                                                            >  1:3417286/1‑149 (MQ=255)
                                                                          gCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAgcctgcc                                                                            <  2:1958318/149‑1 (MQ=255)
                                                                           cccTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAgcctgcct                                                                           >  2:126259/1‑149 (MQ=255)
                                                                             cTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGa                                                                         >  1:1386035/1‑149 (MQ=255)
                                                                             cTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGa                                                                         >  2:1488641/1‑149 (MQ=255)
                                                                             cTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGa                                                                         >  1:2768730/1‑149 (MQ=255)
                                                                              tGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAg                                                                        <  1:2251832/149‑1 (MQ=255)
                                                                               gTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGc                                                                       <  1:1400462/149‑1 (MQ=255)
                                                                                ttAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGCAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACCTTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGcc                                                                      <  1:3271552/149‑1 (MQ=255)
                                                                                 tAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCt                                                                     >  1:532186/1‑149 (MQ=255)
                                                                                   gCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTGGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTAt                                                                   >  2:3650991/1‑149 (MQ=255)
                                                                                   gCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTAt                                                                   <  1:3727467/149‑1 (MQ=255)
                                                                                   gCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTAt                                                                   >  2:2517587/1‑149 (MQ=255)
                                                                                    cGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATc                                                                  <  2:2684739/149‑1 (MQ=255)
                                                                                     gAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATcc                                                                 <  1:1296088/149‑1 (MQ=255)
                                                                                      aaCGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCa                                                                >  2:2245892/1‑149 (MQ=255)
                                                                                      aaCGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCa                                                                <  1:3210034/149‑1 (MQ=255)
                                                                                       aCGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAg                                                               >  1:1406555/1‑149 (MQ=255)
                                                                                       aCGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAg                                                               >  1:1777766/1‑149 (MQ=255)
                                                                                         ggTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGtt                                                             >  2:382527/1‑149 (MQ=255)
                                                                                         ggTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGtt                                                             <  1:3372493/149‑1 (MQ=255)
                                                                                         ggTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGtt                                                             >  2:2590696/1‑149 (MQ=255)
                                                                                           tGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTAc                                                           <  1:163350/149‑1 (MQ=255)
                                                                                           tGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTAc                                                           <  1:2236080/149‑1 (MQ=255)
                                                                                           tGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTAc                                                           <  1:3504097/149‑1 (MQ=255)
                                                                                            gAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACt                                                          >  2:433435/1‑149 (MQ=255)
                                                                                             aaTTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACAGCGTGACTGGGGCGTTGATGTATTTACCGTTCCGGGTAAAGGATACATCCTGCCTGAGCCTATCCAGTTACtt                                                         >  2:1802801/1‑149 (MQ=255)
                                                                                               tttCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTaa                                                       <  1:847021/149‑1 (MQ=255)
                                                                                               tttCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTaa                                                       <  1:847969/149‑1 (MQ=255)
                                                                                                 tCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAAt                                                      <  2:96262/148‑1 (MQ=255)
                                                                                                 tCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAAt                                                      >  1:96262/1‑148 (MQ=255)
                                                                                                  cACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGc                                                    >  1:3910160/1‑149 (MQ=255)
                                                                                                   aCTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCt                                                   >  2:737751/1‑149 (MQ=255)
                                                                                                     tctGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTaa                                                 >  2:2959254/1‑149 (MQ=255)
                                                                                                     tctGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTaa                                                 <  1:641174/149‑1 (MQ=255)
                                                                                                     tctGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGAGGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTaa                                                 >  2:3674283/1‑149 (MQ=255)
                                                                                                      ctGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTaaa                                                <  1:2470946/149‑1 (MQ=255)
                                                                                                          cGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGa                                            <  1:3823334/149‑1 (MQ=255)
                                                                                                             gCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGatat                                         <  2:4089805/149‑1 (MQ=255)
                                                                                                                gTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTgg                                      <  1:45036/149‑1 (MQ=255)
                                                                                                                 ttGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTggg                                     >  1:2058482/1‑149 (MQ=255)
                                                                                                                 ttGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTggg                                     >  1:4091939/1‑149 (MQ=255)
                                                                                                                   gggTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTc                                   <  1:1509958/149‑1 (MQ=255)
                                                                                                                     gTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAg                                 >  2:167862/1‑149 (MQ=255)
                                                                                                                     gTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAg                                 >  2:3507989/1‑149 (MQ=255)
                                                                                                                      tGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGc                                <  1:2331492/149‑1 (MQ=255)
                                                                                                                        aaaCGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTg                              >  2:2854213/1‑149 (MQ=255)
                                                                                                                         aaCGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTgg                             <  2:1757714/149‑1 (MQ=255)
                                                                                                                         aaCGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTgg                             >  2:674185/1‑149 (MQ=255)
                                                                                                                            gCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATg                          >  2:2186470/1‑149 (MQ=255)
                                                                                                                            gCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATg                          >  2:2558373/1‑149 (MQ=255)
                                                                                                                               gggAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCg                       >  1:3684431/1‑149 (MQ=255)
                                                                                                                               gggAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCg                       >  1:3034235/1‑149 (MQ=255)
                                                                                                                                 gAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGt                     <  2:2280275/149‑1 (MQ=255)
                                                                                                                                   aTGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGTAg                   >  2:1328344/1‑149 (MQ=255)
                                                                                                                                   aTGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGTAg                   >  2:2021682/1‑149 (MQ=255)
                                                                                                                                   aTGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGTAg                   >  1:1444777/1‑149 (MQ=255)
                                                                                                                                     gAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGTAgtg                 <  2:1889758/149‑1 (MQ=255)
                                                                                                                                      aGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGTAgtgt                <  1:3252332/149‑1 (MQ=255)
                                                                                                                                                    taataaACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGTAGTGTAGCCGTGCTGCCAg  >  2:3710748/1‑149 (MQ=255)
                                                                                                                                                    taataaACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGTAGTGTAGCCGTGCTGCCAg  >  2:2081289/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ATTGGTGCATCAGGTGAAGTGAGCGCAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGTAGTGTAGCCGTGCTGCCAG  >  CP009273/4162969‑4163265

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: