Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,648,870 A→T 12.7% W177R (TGG→AGG)  rspA ← bifunctional D‑altronate/D‑mannonate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,648,8700AT12.7% 993.7 / 119.0 331W177R (TGG→AGG) rspAbifunctional D‑altronate/D‑mannonate dehydratase
Reads supporting (aligned to +/- strand):  ref base A (149/139);  new base T (21/21);  total (171/160)
Fisher's exact test for biased strand distribution p-value = 8.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.19e-01

Only 100 of 332 total aligned reads displayed.

CAAAGCGCGCCGCTTCAATAGGCGTTAAGCGATGGTGCATGTCATGCAGCAAATGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATAC  >  CP009273/1648723‑1649019
                                                                                                                                                   |                                                                                                                                                     
 aaaGCGCGCCGCTTCAATAGGCGTTAAGCGATGGTGCATGTCATGCAGCAAATGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTg                                                                                                                                                    <  1:3910659/150‑1 (MQ=255)
   aGCGCGCCGCTTCAATAGGCGTTAAGCGATGGTGCATGTCATGCAGCAAATGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTc                                                                                                                                                  <  1:3694323/150‑1 (MQ=255)
      gcgcCGCTTCAATAGGCGTTAAGCGATGGTGCATGTCATGCAGCAAATGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCtt                                                                                                                                               >  2:1745664/1‑150 (MQ=255)
        gccgcTTCAATAGGCGTTAAGCGATGGTGCATGTCATGCAGCAAATGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTg                                                                                                                                             <  2:3099941/150‑1 (MQ=255)
           gcTTCAATAGGCGTTAAGCGATGGTGCATGTCATGCAGCAAATGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGtt                                                                                                                                          >  2:3417215/1‑150 (MQ=255)
            cgtCAATAGGCGTTAAGCGATGGTGCATGTCATGCAGCAAATGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTg                                                                                                                                         <  2:863465/148‑1 (MQ=255)
            cTTCAATAGGCGTTAAGCGATGGTGCATGTCATGCAGCAAATGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCTCTGTCCTTTGGTTg                                                                                                                                         >  1:1151277/1‑150 (MQ=255)
            cTTCAATAGGCGTTAAGCGATGGTGCATGTCATGCAGCAAATGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTg                                                                                                                                         <  1:4462281/150‑1 (MQ=255)
             ttCAATAGGCGTTAAGCGATGGTGCATGTCATGCAGCAAATGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGc                                                                                                                                        >  1:2517324/1‑150 (MQ=255)
               cAATAGGCGTTAAGCGATGGTGCATGTCATGCAGCAAATGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCgg                                                                                                                                      <  1:74447/150‑1 (MQ=255)
                aaTAGGCGTTAAGCGATGGTGCATGTCATGCAGCAAATGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCTCTGTCCTTTGGTTGCggg                                                                                                                                     >  2:5005755/1‑150 (MQ=255)
                  gaGGCGTTAAGCGATGGTGCATGTCATGCAGCAAATGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGtt                                                                                                                                   <  1:2806480/149‑1 (MQ=255)
                      cGTTAAGCGATGGTGCATGTCATGCAGCAAATGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATa                                                                                                                               >  1:2277165/1‑150 (MQ=255)
                       gTTAAGCGATGGTGCATGTCATGCAGCAAATGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATaa                                                                                                                              <  1:2111599/150‑1 (MQ=255)
                       gTTAAGCGATGGTGCATGTCATGCAGCAAATGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATaa                                                                                                                              >  2:3828675/1‑150 (MQ=255)
                        ttAAGCGATGGTGCATGTCATGCAGCAAATGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAg                                                                                                                             <  2:1101302/150‑1 (MQ=255)
                        ttAAGCGATGGTGCATGTCATGCAGCAAATGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAg                                                                                                                             <  2:3920001/150‑1 (MQ=255)
                         tAAGCGATGGTGCATGTCATGCAGCAAATGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGc                                                                                                                            >  1:119095/1‑150 (MQ=255)
                         tAAGCGATGGTGCATGTCATGCAGCAAATGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGc                                                                                                                            >  1:3317097/1‑150 (MQ=255)
                             cGATGGTGCATGTCATGCAGCAAATGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGa                                                                                                                        >  1:681893/1‑150 (MQ=255)
                             cGATGGTGCATGTCATGCAGCAAATGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGa                                                                                                                        >  1:2697380/1‑150 (MQ=255)
                              gATGGTGCATGTCATGCAGCAAATGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCTCTGTCCTTTGGTTGCGGGTTCATAAGCCAGac                                                                                                                       >  1:458424/1‑150 (MQ=255)
                              gATGGTGCATGTCATGCAGCAAATGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGac                                                                                                                       <  1:1554306/150‑1 (MQ=255)
                                tGGTGCATGTCATGCAGCAAATGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGacct                                                                                                                     <  2:591806/150‑1 (MQ=255)
                                tGGTGCATGTCATGCAGCAAATGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCAAGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGacct                                                                                                                     <  2:2511279/150‑1 (MQ=255)
                                 ggTGCATGTCATGCAGCAAATGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGacctt                                                                                                                    >  2:4285541/1‑150 (MQ=255)
                                   tGCATGTCATGCAGCAAATGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGaccttta                                                                                                                  <  2:3937962/150‑1 (MQ=255)
                                   tGCATGTCATGCAGCAAATGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGaccttta                                                                                                                  <  2:4214101/150‑1 (MQ=255)
                                    gCATGTCATGCAGCAAATGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCTCTGTCCTTTGGTTGCGGGTTCATAAGCCAGacctttac                                                                                                                 <  2:2840754/150‑1 (MQ=255)
                                    gCATGTCATGCAGCAAATGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGacctttac                                                                                                                 >  1:4415370/1‑150 (MQ=255)
                                     cATGTCATGCAGCAAATGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGacctttacc                                                                                                                <  1:2100779/150‑1 (MQ=255)
                                      aTGTCATGCAGCAAATGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGacctttacct                                                                                                               <  2:3162596/150‑1 (MQ=255)
                                      aTGTCATGCAGCAAATGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGacctttacct                                                                                                               >  1:4906103/1‑150 (MQ=255)
                                      aTGTCATGCAGCAAATGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGacctttacct                                                                                                               >  2:5037725/1‑150 (MQ=255)
                                       tGTCATGCAGCAAATGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGacctttacctt                                                                                                              <  2:2060852/150‑1 (MQ=255)
                                        gTCATGCAGCAAATGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGacctttaccttt                                                                                                             <  2:814219/150‑1 (MQ=255)
                                              cagcaAATGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACAt                                                                                                       <  2:4745406/150‑1 (MQ=255)
                                                gcaAATGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGc                                                                                                     <  1:835830/150‑1 (MQ=255)
                                                 caAATGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGcc                                                                                                    >  1:545780/1‑150 (MQ=255)
                                                  aaaTGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCg                                                                                                   >  1:4695476/1‑150 (MQ=255)
                                                  aaaTGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCg                                                                                                   >  1:1260749/1‑150 (MQ=255)
                                                    aTGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCTCTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTa                                                                                                 >  1:884191/1‑150 (MQ=255)
                                                     tGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAg                                                                                                <  1:4475674/150‑1 (MQ=255)
                                                       ttCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAggt                                                                                              >  2:2161872/1‑150 (MQ=255)
                                                        tCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAggtg                                                                                             <  2:4491680/150‑1 (MQ=255)
                                                         cATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCTCTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAggtgg                                                                                            >  2:932188/1‑150 (MQ=255)
                                                             aaaaCCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGtttt                                                                                        >  1:4879712/1‑150 (MQ=255)
                                                              aaaCCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTc                                                                                       <  2:2280356/150‑1 (MQ=255)
                                                                aCCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCAt                                                                                     >  2:525457/1‑150 (MQ=255)
                                                                aCCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGACTCTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCAt                                                                                     <  1:2389447/150‑1 (MQ=255)
                                                                   aaaCTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCTCTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATAcc                                                                                  <  1:5005755/150‑1 (MQ=255)
                                                                   aaaCTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATAcc                                                                                  >  2:3281460/1‑150 (MQ=255)
                                                                    aaCTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCa                                                                                 <  2:4316577/150‑1 (MQ=255)
                                                                       ttGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAggg                                                                              <  1:1458534/150‑1 (MQ=255)
                                                                        tGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGa                                                                             <  2:2885944/150‑1 (MQ=255)
                                                                        tGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGa                                                                             <  2:214559/150‑1 (MQ=255)
                                                                        tGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGa                                                                             >  2:2613840/1‑150 (MQ=255)
                                                                          ttACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTACGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGAtt                                                                           >  2:4417029/1‑150 (MQ=255)
                                                                           tACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTc                                                                          <  2:1949051/150‑1 (MQ=255)
                                                                           tACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTc                                                                          <  2:483985/150‑1 (MQ=255)
                                                                              gTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCTCTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGc                                                                       <  2:886618/150‑1 (MQ=255)
                                                                                aCCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCTCTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCAc                                                                     <  1:1415556/150‑1 (MQ=255)
                                                                                aCCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCAc                                                                     <  2:4460853/150‑1 (MQ=255)
                                                                                    cgTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCa                                                                 >  2:4564791/1‑150 (MQ=255)
                                                                                    cgTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCGTTCTACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCa                                                                 >  1:2401745/1‑150 (MQ=255)
                                                                                     gTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGTCCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCAc                                                                <  1:1813984/150‑1 (MQ=255)
                                                                                     gTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCAc                                                                >  1:4522234/1‑150 (MQ=255)
                                                                                     gTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCAc                                                                <  1:5080537/150‑1 (MQ=255)
                                                                                             aTTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATgg                                                        >  2:4298718/1‑150 (MQ=255)
                                                                                               ttCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCt                                                      >  2:4656290/1‑150 (MQ=255)
                                                                                               ttCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCt                                                      >  2:733381/1‑150 (MQ=255)
                                                                                                tCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCAATGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATTCCTTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGTCtt                                                     <  1:1043770/150‑1 (MQ=255)
                                                                                                  ggCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTg                                                   <  2:1643478/150‑1 (MQ=255)
                                                                                                   gCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGa                                                  >  2:2531683/1‑150 (MQ=255)
                                                                                                      tGACATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCCCTGCACGCGGATGGCTTTGTATc                                               <  2:131033/150‑1 (MQ=255)
                                                                                                      tGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATc                                               <  1:3491174/150‑1 (MQ=255)
                                                                                                       gAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATTCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATcc                                              >  2:2998649/1‑150 (MQ=255)
                                                                                                        aaaTCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCa                                             >  1:3567309/1‑150 (MQ=255)
                                                                                                          aTCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATAACAGGGATTCCGCACTGCACGCGGATGTCTTTGAATCCAAg                                           >  2:359119/1‑150 (MQ=255)
                                                                                                           tCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCTCTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGc                                          >  2:3950646/1‑150 (MQ=255)
                                                                                                            cGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGct                                         <  1:1498358/150‑1 (MQ=255)
                                                                                                             gAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCTCTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACAAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAtctc                                        >  1:2688945/1‑150 (MQ=255)
                                                                                                             gAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGctc                                        <  1:1944540/150‑1 (MQ=255)
                                                                                                              aGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGctct                                       >  1:4151081/1‑150 (MQ=255)
                                                                                                               ggTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCTCTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCtt                                      <  1:2605565/150‑1 (MQ=255)
                                                                                                                 tATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTtga                                    >  2:2805782/1‑150 (MQ=255)
                                                                                                                 tATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTtga                                    <  1:1288175/150‑1 (MQ=255)
                                                                                                                  aTTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTtgat                                   >  1:247308/1‑150 (MQ=255)
                                                                                                                   tttCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTtgatg                                  >  2:3333495/1‑150 (MQ=255)
                                                                                                                          gTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCa                           <  2:2697246/150‑1 (MQ=255)
                                                                                                                            cGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATa                         <  1:1721145/150‑1 (MQ=255)
                                                                                                                               ccACAGCTGCTCTTCCGGCCTCTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAatc                      <  2:2314612/150‑1 (MQ=255)
                                                                                                                                  caGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAatcatc                   >  2:1497074/1‑150 (MQ=255)
                                                                                                                                   aGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATcc                  <  1:1111559/150‑1 (MQ=255)
                                                                                                                                    gctgctCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCa                 <  1:29321/150‑1 (MQ=255)
                                                                                                                                      tgctCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCaga               <  1:1760084/150‑1 (MQ=255)
                                                                                                                                         tctTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCt            >  2:4584618/1‑150 (MQ=255)
                                                                                                                                           ttCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTtc          >  1:3108414/1‑150 (MQ=255)
                                                                                                                                               ggCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCAAGAGCTtcatca      >  1:3201430/1‑150 (MQ=255)
                                                                                                                                                 ccACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAAt    <  1:3798030/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CAAAGCGCGCCGCTTCAATAGGCGTTAAGCGATGGTGCATGTCATGCAGCAAATGTTCATTAAAACCAAACTTGTTACGTACCGCGTCAAACAATTTCGGCATGAAATCGAGGTATTTCTCCGTCGACCACAGCTGCTCTTCCGGCCACTGTCCTTTGGTTGCGGGTTCATAAGCCAGACCTTTACCTTTCGACATGCCGTAGGTGGTTTTCATACCAGGGATTCCGCACTGCACGCGGATGGCTTTGAATCCAAGCTCTTGATGACGGGCATAATCATCCAGAGCTTCATCAATAC  >  CP009273/1648723‑1649019

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: