Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,898,210 G→T 7.4% L250I (CTA→ATA)  bglF ← fused beta‑glucoside‑specific PTS enzymes: IIA component/IIB component/IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,898,2100GT 7.4% 1246.7 / 53.2 326L250I (CTA→ATA) bglFfused beta‑glucoside‑specific PTS enzymes: IIA component/IIB component/IIC component
Reads supporting (aligned to +/- strand):  ref base G (143/159);  new base T (7/17);  total (150/176)
Fisher's exact test for biased strand distribution p-value = 9.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 328 total aligned reads displayed.

ATTTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAG  >  CP009273/3898062‑3898357
                                                                                                                                                    |                                                                                                                                                   
    gCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAAt                                                                                                                                                 <  2:2792778/149‑1 (MQ=255)
     ccAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATg                                                                                                                                                <  1:870497/149‑1 (MQ=255)
       aGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATATCAATGgt                                                                                                                                              >  1:4257083/1‑149 (MQ=255)
       aGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGgt                                                                                                                                              >  1:3226798/1‑149 (MQ=255)
          aGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGtgtg                                                                                                                                           >  1:643179/1‑149 (MQ=255)
           gccgccCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTga                                                                                                                                          >  1:1947961/1‑149 (MQ=255)
            ccgccCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATATCAATGGTGTgaa                                                                                                                                         <  2:1427569/149‑1 (MQ=255)
            ccgccCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTgaa                                                                                                                                         <  2:2461242/149‑1 (MQ=255)
            ccgccCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTgaa                                                                                                                                         <  2:5144858/149‑1 (MQ=255)
            ccgccCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTgaa                                                                                                                                         <  2:3361142/149‑1 (MQ=255)
             cgccCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTgaag                                                                                                                                        <  2:2145255/149‑1 (MQ=255)
                ccATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGaaa                                                                                                                                     >  1:2773663/1‑149 (MQ=255)
                 cATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAAt                                                                                                                                    <  2:15714/149‑1 (MQ=255)
                 cATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAAt                                                                                                                                    >  2:1688367/1‑149 (MQ=255)
                    tACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAAtttt                                                                                                                                 <  2:3403422/149‑1 (MQ=255)
                       cgcgCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGa                                                                                                                              <  1:1980944/149‑1 (MQ=255)
                          gcCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTg                                                                                                                           >  1:1367793/1‑149 (MQ=255)
                           cccGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGc                                                                                                                          <  2:2073252/149‑1 (MQ=255)
                           cccGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGc                                                                                                                          <  2:4666879/149‑1 (MQ=255)
                             cGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCg                                                                                                                        <  2:2635548/149‑1 (MQ=255)
                               cAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGAc                                                                                                                      <  1:2363576/149‑1 (MQ=255)
                                 aaTGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACgg                                                                                                                    <  2:3275050/149‑1 (MQ=255)
                                  aTGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGt                                                                                                                   <  2:3871296/149‑1 (MQ=255)
                                    gCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTaa                                                                                                                 >  1:3580531/1‑149 (MQ=255)
                                      aGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAAcc                                                                                                               >  1:5193411/1‑149 (MQ=255)
                                       ggAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCa                                                                                                              <  1:1032077/149‑1 (MQ=255)
                                            cGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATATCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCAt                                                                                                         <  2:4922708/149‑1 (MQ=255)
                                             gCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCAtt                                                                                                        <  2:2296878/149‑1 (MQ=255)
                                              ccTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTa                                                                                                       >  1:2516687/1‑149 (MQ=255)
                                               cTGATACAGCCAGAGATAACCCGCGGCAATCAGTTCGCTGATCCAGGTTGATAGCGGCCCAACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGAGAAGAAATTATTGATTGCCGACGGTAACCACGCATTaa                                                                                                      <  2:2917133/149‑1 (MQ=255)
                                               cTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTaa                                                                                                      >  1:1538568/1‑149 (MQ=255)
                                                tGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAg                                                                                                     <  2:5114285/149‑1 (MQ=255)
                                                 gATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATATCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGt                                                                                                    <  1:4870230/149‑1 (MQ=255)
                                                 gATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGt                                                                                                    <  2:2370754/149‑1 (MQ=255)
                                                       gCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGt                                                                                              <  2:3469145/149‑1 (MQ=255)
                                                        ccAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGtt                                                                                             <  1:2162367/149‑1 (MQ=255)
                                                         cAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTc                                                                                            >  2:381055/1‑149 (MQ=255)
                                                             gaTAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGa                                                                                        >  2:2844951/1‑149 (MQ=255)
                                                               tAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAAt                                                                                      >  2:2613907/1‑149 (MQ=255)
                                                               tAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAAt                                                                                      >  2:576087/1‑149 (MQ=255)
                                                                 cccGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAAtgc                                                                                    <  2:2117813/148‑1 (MQ=255)
                                                                 aCCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAAtgc                                                                                    <  2:662698/149‑1 (MQ=255)
                                                                  ccGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAAtgct                                                                                   <  1:3805703/149‑1 (MQ=255)
                                                                         cAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAAc                                                                            <  2:1213450/149‑1 (MQ=255)
                                                                          aaTCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAAcc                                                                           >  1:768686/1‑149 (MQ=255)
                                                                           aTCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCa                                                                          <  1:2472305/149‑1 (MQ=255)
                                                                           aTCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCa                                                                          >  1:942538/1‑149 (MQ=255)
                                                                            tCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAg                                                                         <  2:2018666/149‑1 (MQ=255)
                                                                            tCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAg                                                                         >  1:3208239/1‑149 (MQ=255)
                                                                            tCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAg                                                                         <  1:207705/149‑1 (MQ=255)
                                                                                ttCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAg                                                                     >  1:1085794/1‑149 (MQ=255)
                                                                                  cGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGaa                                                                   >  1:1555921/1‑149 (MQ=255)
                                                                                  cGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGaa                                                                   >  1:4151776/1‑149 (MQ=255)
                                                                                  cGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGaa                                                                   >  1:2954113/1‑149 (MQ=255)
                                                                                    cTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGaaaa                                                                 >  1:1146844/1‑149 (MQ=255)
                                                                                        tCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAaata                                                             <  2:2335674/149‑1 (MQ=255)
                                                                                            ggTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATgg                                                         <  2:768686/149‑1 (MQ=255)
                                                                                             gTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATggg                                                        >  1:2604484/1‑149 (MQ=255)
                                                                                             gTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATggg                                                        <  1:1837755/149‑1 (MQ=255)
                                                                                              ttGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGa                                                       <  1:1104504/149‑1 (MQ=255)
                                                                                                gATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAAt                                                     <  1:3902598/149‑1 (MQ=255)
                                                                                                  tAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATaa                                                   <  1:3484982/149‑1 (MQ=255)
                                                                                                  tAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATaa                                                   >  2:3978113/1‑149 (MQ=255)
                                                                                                    gCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAAcc                                                 >  1:3992914/1‑149 (MQ=255)
                                                                                                      ggCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGa                                               <  1:558528/149‑1 (MQ=255)
                                                                                                       gCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGAt                                              >  1:4650862/1‑149 (MQ=255)
                                                                                                       gCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGAt                                              <  1:3679269/149‑1 (MQ=255)
                                                                                                        ccccACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATg                                             <  1:1419027/149‑1 (MQ=255)
                                                                                                        ccccACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATg                                             <  1:46100/149‑1 (MQ=255)
                                                                                                        ccccACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATg                                             >  2:2648364/1‑149 (MQ=255)
                                                                                                         cccACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGa                                            <  1:2919429/149‑1 (MQ=255)
                                                                                                               agcagAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTa                                      >  1:201437/1‑149 (MQ=255)
                                                                                                                 cagAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAAt                                    >  2:2991507/1‑149 (MQ=255)
                                                                                                                 cagAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAAt                                    <  2:4433808/149‑1 (MQ=255)
                                                                                                                  agAAAGGTGACGGGTGTGATAACCATCAGACATATCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAAtt                                   <  2:924235/149‑1 (MQ=255)
                                                                                                                     aaGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTcaa                                >  1:3726402/1‑149 (MQ=255)
                                                                                                                     aaGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTcaa                                <  1:380034/149‑1 (MQ=255)
                                                                                                                       ggTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTcaaca                              >  1:1797141/1‑149 (MQ=255)
                                                                                                                         tGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAAt                            <  2:1119271/149‑1 (MQ=255)
                                                                                                                           aCGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAAtgt                          <  1:5028314/149‑1 (MQ=255)
                                                                                                                             gggTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGa                        >  1:4847020/1‑149 (MQ=255)
                                                                                                                               gtgtGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGAcc                      >  1:3013959/1‑149 (MQ=255)
                                                                                                                                 gtgATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCgg                    <  2:2156401/149‑1 (MQ=255)
                                                                                                                                  tgATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCCTTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGa                   >  2:4634591/1‑149 (MQ=255)
                                                                                                                                   gATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGaa                  >  2:714434/1‑149 (MQ=255)
                                                                                                                                    aTAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAAt                 >  2:852231/1‑149 (MQ=255)
                                                                                                                                    aTAAACATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAAt                 <  1:3912064/149‑1 (MQ=255)
                                                                                                                                     tAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATa                <  1:3171010/149‑1 (MQ=255)
                                                                                                                                      aaCCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATAc               >  2:4303753/1‑149 (MQ=255)
                                                                                                                                       aCCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATAcc              >  2:1976897/1‑149 (MQ=255)
                                                                                                                                        ccATCAGACATAGCAATGGTGTGTAGACATTTTTGATTGCCGACGGTAACCACGCATTTAGTCGGCTTTCCAGAATGCTGCACATCCAGGCAGTTAAAATTATGGTAATAACCGTTGACGAGTTATTCTACAATGTGAGCGGAATAccc             <  2:501827/149‑1 (MQ=255)
                                                                                                                                        ccATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATAccc             >  1:3967671/1‑149 (MQ=255)
                                                                                                                                         cATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCa            >  1:513941/1‑149 (MQ=255)
                                                                                                                                          aTCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAg           >  1:2550712/1‑149 (MQ=255)
                                                                                                                                            cAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGACTACCCAGGa         >  2:4797086/1‑149 (MQ=255)
                                                                                                                                             aGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGaa        <  2:4571050/149‑1 (MQ=255)
                                                                                                                                             aGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGaa        >  1:2301690/1‑149 (MQ=255)
                                                                                                                                                  tAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCa   >  2:1969993/1‑149 (MQ=255)
                                                                                                                                                  tAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCa   <  2:4159057/149‑1 (MQ=255)
                                                                                                                                                   aGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAg  >  1:3310370/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATTTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAG  >  CP009273/3898062‑3898357

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: