Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,163,117 T→A 100% 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,1170TA96.5% 736.0 / 11.7 229I36I (ATT→ATAbirAbifunctional biotin‑[acetylCoA carboxylase] holoenzyme synthetase/ DNA‑binding transcriptional repressor, bio‑5'‑AMP‑binding
Reads supporting (aligned to +/- strand):  ref base T (4/4);  new base A (117/104);  total (121/108)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 238 total aligned reads displayed.

GGTGCATCAGGTGAAGTGAGCGCAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGTAGTGTAGCCGTGCTGC  >  CP009273/4162972‑4163262
                                                                                                                                                 |                                                                                                                                                 
    cATCAGGTGAAGTGAGCGCAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTAtaaataaa                                                                                                                                             >  1:27167/1‑148 (MQ=255)
    cATCAGGTGAAGTGAGCGCAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTAtaaataaa                                                                                                                                             <  1:1846991/148‑1 (MQ=255)
    cATCAGGTGAAGTGAGCGCAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATaaa                                                                                                                                             >  1:895371/1‑148 (MQ=255)
      tCAGGTGAAGTGAGCGCAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAaca                                                                                                                                           >  1:855841/1‑148 (MQ=255)
        aGGTGAAGTGAGCGCAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAAcaca                                                                                                                                         <  2:2169013/148‑1 (MQ=255)
         ggTGAAGTGAGCGCAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACAt                                                                                                                                        <  1:2329504/148‑1 (MQ=255)
          gTGAAGTGAGCGCAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACAtt                                                                                                                                       >  1:442475/1‑148 (MQ=255)
          gTGAAGTGAGCGCAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACAtt                                                                                                                                       >  1:2402881/1‑148 (MQ=255)
                   gcgcAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTg                                                                                                                              >  1:1294002/1‑148 (MQ=255)
                   gcgcAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTg                                                                                                                              >  2:3381312/1‑148 (MQ=255)
                     gcAGTGGAGACAATTTCATGAAGGATAACACCGTGTCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCg                                                                                                                            <  2:1928851/148‑1 (MQ=255)
                     gcAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCg                                                                                                                            >  1:2907861/1‑148 (MQ=255)
                         tGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGAc                                                                                                                        >  1:1941716/1‑148 (MQ=255)
                          ggAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACt                                                                                                                       >  1:870721/1‑148 (MQ=255)
                                  tttCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGtt                                                                                                               >  2:1534408/1‑148 (MQ=255)
                                  tttCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGtt                                                                                                               >  1:200369/1‑148 (MQ=255)
                                   ttCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTg                                                                                                              >  1:2641197/1‑148 (MQ=255)
                                      aTGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATg                                                                                                           <  2:3132627/148‑1 (MQ=255)
                                      aTGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCAGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATg                                                                                                           <  1:942536/148‑1 (MQ=255)
                                        gAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTc                                                                                                         <  2:1378243/148‑1 (MQ=255)
                                         aaGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCt                                                                                                        >  1:1572937/1‑148 (MQ=255)
                                           ggATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCttt                                                                                                      >  1:2853088/1‑148 (MQ=255)
                                           ggATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCttt                                                                                                      <  2:2204087/148‑1 (MQ=255)
                                               aaCACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCg                                                                                                  <  2:3467940/148‑1 (MQ=255)
                                                acacCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGt                                                                                                 >  2:634505/1‑148 (MQ=255)
                                                acacCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGt                                                                                                 >  2:3073494/1‑148 (MQ=255)
                                                 cacCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGtt                                                                                                >  2:1689666/1‑148 (MQ=255)
                                                 cacCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGtt                                                                                                >  2:2069976/1‑148 (MQ=255)
                                                   ccGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTcc                                                                                              >  1:2701358/1‑148 (MQ=255)
                                                      tGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCggg                                                                                           <  1:2812881/148‑1 (MQ=255)
                                                      tGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCggg                                                                                           <  1:2132999/148‑1 (MQ=255)
                                                      tGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTGGATGTCTTTACCGTTCCggg                                                                                           >  2:505896/1‑148 (MQ=255)
                                                       gCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGt                                                                                          >  2:2736346/1‑148 (MQ=255)
                                                       gCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGt                                                                                          >  2:2780163/1‑148 (MQ=255)
                                                          aCTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTaaa                                                                                       <  2:1494930/148‑1 (MQ=255)
                                                           cTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAg                                                                                      <  2:551339/148‑1 (MQ=255)
                                                           cTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAg                                                                                      >  1:3323684/1‑148 (MQ=255)
                                                            tGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAgg                                                                                     <  2:1724057/148‑1 (MQ=255)
                                                            tGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAgg                                                                                     >  2:1790590/1‑148 (MQ=255)
                                                            tGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAgg                                                                                     <  2:959317/148‑1 (MQ=255)
                                                             gAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGa                                                                                    <  2:3025764/148‑1 (MQ=255)
                                                               aaTTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATa                                                                                  >  1:3027316/1‑148 (MQ=255)
                                                                attgattgCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATAc                                                                                 >  1:1622605/1‑148 (MQ=255)
                                                                attgattgCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATAc                                                                                 <  2:1670039/148‑1 (MQ=255)
                                                                attgattgCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATAc                                                                                 >  1:2654215/1‑148 (MQ=255)
                                                                 ttgattgCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACa                                                                                >  1:2135888/1‑148 (MQ=255)
                                                                     ttgCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAgcct                                                                            >  2:2106905/1‑148 (MQ=255)
                                                                      tgCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAgcctg                                                                           >  1:2457888/1‑148 (MQ=255)
                                                                      tgCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTAGAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAgcctg                                                                           >  1:692600/1‑148 (MQ=255)
                                                                        cccTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAgcctgcc                                                                         <  1:975068/148‑1 (MQ=255)
                                                                        cccTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAgcctgcc                                                                         <  2:1387119/148‑1 (MQ=255)
                                                                        cccTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAgcctgcc                                                                         <  1:110331/148‑1 (MQ=255)
                                                                          cTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTCTAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGcctgcctg                                                                       >  2:2917124/1‑148 (MQ=255)
                                                                          cTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGcctgcctg                                                                       >  1:1350615/1‑148 (MQ=255)
                                                                          cTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGcctgcctg                                                                       <  2:3259807/148‑1 (MQ=255)
                                                                          cTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGcctgcctg                                                                       >  2:1166709/1‑148 (MQ=255)
                                                                            gTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAg                                                                     >  2:2983577/1‑148 (MQ=255)
                                                                             ttAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGc                                                                    <  1:1035381/148‑1 (MQ=255)
                                                                               aGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCt                                                                  <  1:2197174/148‑1 (MQ=255)
                                                                               aGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTAAGCCt                                                                  <  2:3243871/148‑1 (MQ=255)
                                                                                gCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTa                                                                 <  2:2002197/148‑1 (MQ=255)
                                                                                  gAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATc                                                               <  1:3073494/148‑1 (MQ=255)
                                                                                     cGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAg                                                            <  2:354131/148‑1 (MQ=255)
                                                                                       gTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGtt                                                          >  1:777754/1‑148 (MQ=255)
                                                                                        tGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTa                                                         >  2:1895032/1‑148 (MQ=255)
                                                                                         gAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTAc                                                        <  2:2652254/148‑1 (MQ=255)
                                                                                          aaTTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACt                                                       <  1:495268/148‑1 (MQ=255)
                                                                                                aCTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGc                                                 >  1:3493615/1‑148 (MQ=255)
                                                                                                 ctctGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTCTCCAGTTACTTAATGCt                                                >  2:3109541/1‑148 (MQ=255)
                                                                                                   ctGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTaa                                              <  1:232030/148‑1 (MQ=255)
                                                                                                    tGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTaaa                                             <  1:1127553/148‑1 (MQ=255)
                                                                                                      gCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACa                                           <  2:2283298/148‑1 (MQ=255)
                                                                                                       cGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAg                                          <  1:820706/148‑1 (MQ=255)
                                                                                                        gAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGa                                         <  1:2537084/148‑1 (MQ=255)
                                                                                                         aGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGat                                        >  1:2330546/1‑148 (MQ=255)
                                                                                                         aGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTATTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGat                                        >  1:1229001/1‑148 (MQ=255)
                                                                                                          gCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGata                                       <  2:895219/148‑1 (MQ=255)
                                                                                                              ttGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTgg                                   >  1:503699/1‑148 (MQ=255)
                                                                                                                gggTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGt                                 <  2:27167/148‑1 (MQ=255)
                                                                                                                gggTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGt                                 <  2:728034/148‑1 (MQ=255)
                                                                                                                gggTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGt                                 >  1:409699/1‑148 (MQ=255)
                                                                                                                gggTGAAACGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGt                                 <  1:1486221/148‑1 (MQ=255)
                                                                                                                       aCGCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTgg                          <  2:2862280/148‑1 (MQ=255)
                                                                                                                         gCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGAt                        <  1:3161902/148‑1 (MQ=255)
                                                                                                                         gCTGGGAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGAt                        <  1:1115487/148‑1 (MQ=255)
                                                                                                                            gggAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGc                     <  2:2402881/148‑1 (MQ=255)
                                                                                                                            gggAATGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGc                     <  2:442507/148‑1 (MQ=255)
                                                                                                                                aTGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGTa                 <  1:3239707/148‑1 (MQ=255)
                                                                                                                                 tGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGTAg                >  2:3217248/1‑148 (MQ=255)
                                                                                                                                 tGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGTAg                >  2:2271228/1‑148 (MQ=255)
                                                                                                                                 tGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGTAg                >  1:581176/1‑148 (MQ=255)
                                                                                                                                 tGAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGTAg                >  2:2303745/1‑148 (MQ=255)
                                                                                                                                  gAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGTAgt               <  2:1536298/148‑1 (MQ=255)
                                                                                                                                  gAGCCGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGTAgt               >  2:2567059/1‑148 (MQ=255)
                                                                                                                                      cGGGCGGCTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGTAGTGTAg           <  1:1525067/148‑1 (MQ=255)
                                                                                                                                        ggcggcTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGTAGTGTAGcc         <  1:1844154/148‑1 (MQ=255)
                                                                                                                                            gcTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGTAGTGTAGCCGtgc     <  1:705767/148‑1 (MQ=255)
                                                                                                                                            gcTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGTAGTGTAGCCGtgc     <  1:181747/148‑1 (MQ=255)
                                                                                                                                             cTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGTAGTGTAGCCGtgct    >  2:3315210/1‑148 (MQ=255)
                                                                                                                                             cTATAAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGTAGTGTAGCCGtgct    >  2:697816/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GGTGCATCAGGTGAAGTGAGCGCAGTGGAGACAATTTCATGAAGGATAACACCGTGCCACTGAAATTGATTGCCCTGTTAGCGAACGGTGAATTTCACTCTGGCGAGCAGTTGGGTGAAACGCTGGGAATGAGCCGGGCGGCTATTAATAAACACATTCAGACACTGCGTGACTGGGGCGTTGATGTCTTTACCGTTCCGGGTAAAGGATACAGCCTGCCTGAGCCTATCCAGTTACTTAATGCTAAACAGATATTGGGTCAGCTGGATGGCGGTAGTGTAGCCGTGCTGC  >  CP009273/4162972‑4163262

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: