Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,565,147 A→T 10.1% F572L (TTT→TTA glgB ← 1,4‑alpha‑glucan branching enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,565,1470AT10.1% 993.0 / 74.7 330F572L (TTT→TTAglgB1,4‑alpha‑glucan branching enzyme
Reads supporting (aligned to +/- strand):  ref base A (136/160);  new base T (10/24);  total (146/184)
Fisher's exact test for biased strand distribution p-value = 7.08e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 336 total aligned reads displayed.

TCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAA  >  CP009273/3564998‑3565296
                                                                                                                                                     |                                                                                                                                                     
    aaTCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGAGAGGTTCAGAACGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAAtt                                                                                                                                                   >  2:1741847/1‑150 (MQ=255)
     aTCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTc                                                                                                                                                  <  2:4280152/150‑1 (MQ=255)
     aTCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTc                                                                                                                                                  <  1:880868/150‑1 (MQ=255)
     aTCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTc                                                                                                                                                  <  2:1046833/150‑1 (MQ=255)
     aTCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTc                                                                                                                                                  <  2:2168000/150‑1 (MQ=255)
      tCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCTAATTCg                                                                                                                                                 <  1:4933024/150‑1 (MQ=255)
        cAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGtt                                                                                                                                               <  1:644924/150‑1 (MQ=255)
          gTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTAc                                                                                                                                             <  2:1551556/150‑1 (MQ=255)
             cATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCTAATTCGTTACCCa                                                                                                                                          >  1:1816962/1‑150 (MQ=255)
              aTGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCAt                                                                                                                                         <  1:29659/150‑1 (MQ=255)
               tGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATg                                                                                                                                        <  2:3969377/150‑1 (MQ=255)
                 cATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGaa                                                                                                                                      >  2:3327434/1‑150 (MQ=255)
                  aTTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAAc                                                                                                                                     >  1:175057/1‑150 (MQ=255)
                    tGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACag                                                                                                                                   <  2:2475932/150‑1 (MQ=255)
                       tttATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACagtat                                                                                                                                <  1:2743577/150‑2 (MQ=255)
                        ttATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACagtagt                                                                                                                               <  1:1081695/150‑1 (MQ=255)
                         tATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCTAATTCGTTACCCATGAACAGTAGtt                                                                                                                              >  1:480202/1‑150 (MQ=255)
                         tATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGtt                                                                                                                              >  1:4544020/1‑150 (MQ=255)
                           tggtggCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTc                                                                                                                            >  2:2223631/1‑150 (MQ=255)
                           tggtggCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTc                                                                                                                            <  1:267557/150‑1 (MQ=255)
                            ggtggCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCt                                                                                                                           <  2:594531/150‑1 (MQ=255)
                             gtggCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCtt                                                                                                                          >  2:2127210/1‑150 (MQ=255)
                              tggCGGTAGGTGAGTTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTg                                                                                                                         <  1:3408384/150‑1 (MQ=255)
                              tggCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCTAATTCGTTACCCATGAACAGTAGTTTCTTg                                                                                                                         <  1:1175875/150‑1 (MQ=255)
                                gCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGcc                                                                                                                       <  2:3203229/150‑1 (MQ=255)
                                  ggTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCg                                                                                                                     >  1:624775/1‑150 (MQ=255)
                                  ggTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCg                                                                                                                     >  1:2485139/1‑150 (MQ=255)
                                  ggTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCg                                                                                                                     >  1:1459477/1‑150 (MQ=255)
                                  ggTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCg                                                                                                                     >  1:3070433/1‑150 (MQ=255)
                                  ggTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCg                                                                                                                     >  1:4364101/1‑150 (MQ=255)
                                    tAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCTAATTCGTTACCCATGAACAGTAGTTTCTTGCCCggg                                                                                                                   >  1:1995155/1‑150 (MQ=255)
                                    tAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCggg                                                                                                                   >  2:4270454/1‑150 (MQ=255)
                                    tAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCggg                                                                                                                   >  2:4158317/1‑150 (MQ=255)
                                     aGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGa                                                                                                                  <  1:2609822/150‑1 (MQ=255)
                                       gTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAAt                                                                                                                <  2:2020255/150‑1 (MQ=255)
                                        tGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCTAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATg                                                                                                               <  1:4904751/150‑1 (MQ=255)
                                        tGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATg                                                                                                               <  1:3010883/150‑1 (MQ=255)
                                          aGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGcc                                                                                                             >  1:5216801/1‑150 (MQ=255)
                                          aGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGca                                                                                                             >  2:3947238/1‑149 (MQ=255)
                                           ggTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCTAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGccc                                                                                                            >  2:1897090/1‑150 (MQ=255)
                                           ggTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGccc                                                                                                            >  1:1426077/1‑150 (MQ=255)
                                           ggTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGccc                                                                                                            <  2:4684163/150‑1 (MQ=255)
                                              tCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCcaca                                                                                                         >  1:5013677/1‑150 (MQ=255)
                                              tCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCcaca                                                                                                         <  1:2844022/150‑1 (MQ=255)
                                                  aTCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCa                                                                                                     <  1:3366749/150‑1 (MQ=255)
                                                    cgcgCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGc                                                                                                   >  2:1277469/1‑150 (MQ=255)
                                                       gcACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCAt                                                                                                >  2:3170798/1‑150 (MQ=255)
                                                        cACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCAta                                                                                               <  1:5085489/150‑1 (MQ=255)
                                                          ccAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCAtagt                                                                                             >  1:3596909/1‑150 (MQ=255)
                                                             gACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAgg                                                                                          <  2:3407323/150‑1 (MQ=255)
                                                              aCGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGgc                                                                                         >  2:3738831/1‑150 (MQ=255)
                                                                   ggACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCa                                                                                    <  2:4885303/150‑1 (MQ=255)
                                                                     acacCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAgg                                                                                  >  1:2114921/1‑150 (MQ=255)
                                                                        ccGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTc                                                                               >  2:2496546/1‑150 (MQ=255)
                                                                         cGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTcg                                                                              >  2:457320/1‑150 (MQ=255)
                                                                          gtggtgCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTcgc                                                                             >  1:87241/1‑150 (MQ=255)
                                                                                ccAGTTATCGCCGCCTTCCAACAGATGCCAGTCGATGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACAGATGAACAGTAGTTTCTTGCCCGGGTATGCCCACATCCATCCATAGTAGGCGCGCTGGTTCGCGAAttt                                                                       <  2:1070894/150‑1 (MQ=255)
                                                                                ccAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAAttt                                                                       <  2:4188999/150‑1 (MQ=255)
                                                                                ccAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAAttt                                                                       <  2:3264103/150‑1 (MQ=255)
                                                                                ccAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAAttt                                                                       >  1:3213628/1‑150 (MQ=255)
                                                                                ccAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAAttt                                                                       <  1:3910668/150‑1 (MQ=255)
                                                                                ccAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAAttt                                                                       <  2:1205657/150‑1 (MQ=255)
                                                                                   gTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTg                                                                    >  1:1353990/1‑150 (MQ=255)
                                                                                          ccgccTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGACTGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCg                                                             >  1:3465921/1‑150 (MQ=255)
                                                                                          ccgccTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCg                                                             <  1:3741590/150‑1 (MQ=255)
                                                                                           cgccTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGt                                                            >  2:4604327/1‑150 (MQ=255)
                                                                                             ccTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCTAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCg                                                          <  2:3748587/150‑1 (MQ=255)
                                                                                             ccTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCg                                                          >  2:50236/1‑150 (MQ=255)
                                                                                              cTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGc                                                         >  1:2542254/1‑150 (MQ=255)
                                                                                               ttCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGcc                                                        >  1:3509461/1‑150 (MQ=255)
                                                                                                tCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGccc                                                       >  1:4249712/1‑150 (MQ=255)
                                                                                                 ccAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCg                                                      >  2:2445676/1‑150 (MQ=255)
                                                                                                   aaCAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGc                                                    <  2:487251/150‑1 (MQ=255)
                                                                                                    aCAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCa                                                   >  2:2389149/1‑150 (MQ=255)
                                                                                                         tGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCgg                                              >  2:2332267/1‑150 (MQ=255)
                                                                                                            cAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCTAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCg                                           >  2:349406/1‑150 (MQ=255)
                                                                                                                 gAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAAt                                      <  1:366375/150‑1 (MQ=255)
                                                                                                                 gAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAAt                                      <  1:1865471/150‑1 (MQ=255)
                                                                                                                   ggCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCg                                    <  2:1602819/150‑1 (MQ=255)
                                                                                                                      tGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGACTCGAtt                                 >  1:398661/1‑150 (MQ=255)
                                                                                                                      tGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGAtt                                 <  2:2425271/150‑1 (MQ=255)
                                                                                                                      tGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGAtt                                 <  1:4682721/150‑1 (MQ=255)
                                                                                                                       ggCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGAttt                                <  2:2363313/150‑1 (MQ=255)
                                                                                                                         cGTCATGGGTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGAttttt                              <  1:2362602/150‑1 (MQ=255)
                                                                                                                          gTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGAtttttt                             >  2:2819957/1‑150 (MQ=255)
                                                                                                                          gTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGAtttttt                             >  2:2669570/1‑150 (MQ=255)
                                                                                                                          gTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCAGAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGAtttttt                             >  1:2347816/1‑150 (MQ=255)
                                                                                                                           tCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCCGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTa                            <  2:3340770/150‑1 (MQ=255)
                                                                                                                           tCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTa                            >  2:2959667/1‑150 (MQ=255)
                                                                                                                            cATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTAc                           <  2:4247192/150‑1 (MQ=255)
                                                                                                                               ggTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGt                        >  2:3855011/1‑150 (MQ=255)
                                                                                                                                  tCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGa                     <  2:1426724/150‑1 (MQ=255)
                                                                                                                                  tCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGa                     >  2:2171432/1‑150 (MQ=255)
                                                                                                                                           ggCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGACTCGATTTTTTACCCTGGACCACTtcat            >  1:4057412/1‑150 (MQ=255)
                                                                                                                                             cccTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTtcatca          >  1:3875629/1‑150 (MQ=255)
                                                                                                                                               cTGGGCTAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCAGGCGGTCGAGCATCGATTTTTTACCGTGGACCACTTCATCATg        >  1:1560818/1‑150 (MQ=255)
                                                                                                                                               cTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATg        >  1:194058/1‑150 (MQ=255)
                                                                                                                                                 gggCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCg      <  1:780966/150‑1 (MQ=255)
                                                                                                                                                  ggCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGa     >  2:587401/1‑150 (MQ=255)
                                                                                                                                                   gCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGAc    <  2:2435136/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAA  >  CP009273/3564998‑3565296

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: