Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,564,984 C→A 7.0% E627* (GAA→TAA)  glgB ← 1,4‑alpha‑glucan branching enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,564,9840CA 7.0% 965.4 / 29.7 256E627* (GAA→TAA) glgB1,4‑alpha‑glucan branching enzyme
Reads supporting (aligned to +/- strand):  ref base C (133/105);  new base A (10/8);  total (143/113)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 261 total aligned reads displayed.

CGCCATTTGCCCGGCTGGTTTATGCCGAAGCGATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCAC  >  CP009273/3564836‑3565132
                                                                                                                                                    |                                                                                                                                                    
cGCCATTTGCCCGGCTGGTTTATGCCGAAGCGATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTc                                                                                                                                                      <  1:3460442/149‑1 (MQ=255)
 gCCATTTGCCCGGCTGGTTTATGCCGAAGCGATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCa                                                                                                                                                     >  1:4749613/1‑149 (MQ=255)
            ggCTGGTTTATGCCGAAGCGATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACgg                                                                                                                                          <  1:1712133/149‑1 (MQ=255)
             gCTGGTTTATGCCGAAGCGATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACggg                                                                                                                                         >  2:4027869/1‑149 (MQ=255)
             gCTGGTTTATGCCGAAGCGATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACggg                                                                                                                                         <  2:1015471/149‑1 (MQ=255)
               tGGTTTATGCCGAAGCGATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTc                                                                                                                                       >  2:4782958/1‑149 (MQ=255)
               tGGTTTATGCCGAAGCGATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTc                                                                                                                                       <  2:182976/149‑1 (MQ=255)
                  tttATGCCGAAGCGATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCaaa                                                                                                                                    <  2:164588/149‑1 (MQ=255)
                     aTGCCGAAGCGATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATc                                                                                                                                 >  2:2028568/1‑149 (MQ=255)
                      tGCCGAAGCGATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATcc                                                                                                                                >  2:33524/1‑149 (MQ=255)
                      tGCCGAAGCGATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATcc                                                                                                                                <  1:4469529/149‑1 (MQ=255)
                          gAAGCGATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGtt                                                                                                                            <  2:3031549/149‑1 (MQ=255)
                          gAAGCGATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGtt                                                                                                                            >  1:577625/1‑149 (MQ=255)
                           aaGCGATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGACGTACGGGTCAAAATCCAGTTc                                                                                                                           <  1:1537974/149‑1 (MQ=255)
                            aGCGATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCa                                                                                                                          <  1:1232539/149‑1 (MQ=255)
                            aGCGATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTAAAAGCCGTACGGGTCAAAATCCAGTTCa                                                                                                                          <  2:1696240/149‑1 (MQ=255)
                             gCGATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCAt                                                                                                                         >  1:2240440/1‑149 (MQ=255)
                                aTAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCa                                                                                                                      >  2:3526535/1‑149 (MQ=255)
                                 tAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCAt                                                                                                                     <  2:2887023/149‑1 (MQ=255)
                                   aTCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTg                                                                                                                   >  2:4631175/1‑149 (MQ=255)
                                       tGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCttt                                                                                                               >  2:3151935/1‑149 (MQ=255)
                                        gACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTa                                                                                                              <  1:393902/149‑1 (MQ=255)
                                        gACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTa                                                                                                              <  2:4421198/149‑1 (MQ=255)
                                              gTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTAtggtgg                                                                                                        <  2:545611/149‑1 (MQ=255)
                                                 ccGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCgg                                                                                                     >  2:1905230/1‑149 (MQ=255)
                                                    gCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAg                                                                                                  >  1:2119112/1‑149 (MQ=255)
                                                    gCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCAGTGCTTTATGGTGGCGGTAg                                                                                                  >  1:1737634/1‑149 (MQ=255)
                                                       tAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACCAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGGTCATGCATTGCTTTATGGTGGCGGTAGGTg                                                                                               <  1:2006406/149‑1 (MQ=255)
                                                       tAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTg                                                                                               >  2:2314949/1‑149 (MQ=255)
                                                       tAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTg                                                                                               >  1:367291/1‑149 (MQ=255)
                                                        aaaGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGa                                                                                              >  1:573645/1‑149 (MQ=255)
                                                         aaGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAg                                                                                             >  2:4285416/1‑149 (MQ=255)
                                                          aGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAgg                                                                                            <  2:4518538/149‑1 (MQ=255)
                                                           gTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGt                                                                                           >  2:790461/1‑149 (MQ=255)
                                                            ttACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGtt                                                                                          <  1:4350347/149‑1 (MQ=255)
                                                             tACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTTGGTGAGGTTc                                                                                         <  1:1219910/149‑1 (MQ=255)
                                                              aCTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCa                                                                                        <  1:2420191/149‑1 (MQ=255)
                                                               cTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAg                                                                                       >  2:2383914/1‑149 (MQ=255)
                                                                tGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGa                                                                                      >  2:1187251/1‑149 (MQ=255)
                                                                 ggCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGAt                                                                                     <  1:2314949/149‑1 (MQ=255)
                                                                  gCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATc                                                                                    <  1:4246661/149‑1 (MQ=255)
                                                                  gCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATc                                                                                    >  2:606291/1‑149 (MQ=255)
                                                                      cGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCgcgc                                                                                >  1:4019212/1‑149 (MQ=255)
                                                                      cGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTAAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCgcgc                                                                                <  1:3283370/149‑1 (MQ=255)
                                                                      cGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTAAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCgcgc                                                                                <  1:712060/149‑1 (MQ=255)
                                                                       gatgatTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCa                                                                               >  2:2230718/1‑149 (MQ=255)
                                                                               cGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGc                                                                       <  1:4215101/149‑1 (MQ=255)
                                                                               cGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGc                                                                       <  1:1085465/149‑1 (MQ=255)
                                                                                gTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCt                                                                      >  1:1754237/1‑149 (MQ=255)
                                                                                gTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCt                                                                      <  1:4782958/149‑1 (MQ=255)
                                                                                  tACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTgg                                                                    >  2:3003853/1‑149 (MQ=255)
                                                                                  tACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTgg                                                                    >  1:1809464/1‑149 (MQ=255)
                                                                                    cccTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGac                                                                  <  1:2973635/149‑1 (MQ=255)
                                                                                    cccTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGac                                                                  >  1:3802026/1‑149 (MQ=255)
                                                                                       tctTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACAcc                                                               >  2:368970/1‑149 (MQ=255)
                                                                                          ttATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTAAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCgtg                                                            <  1:361567/149‑1 (MQ=255)
                                                                                           tATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCCTTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCgtgg                                                           >  1:2578172/1‑149 (MQ=255)
                                                                                             tCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCgtggtg                                                         >  2:956673/1‑149 (MQ=255)
                                                                                             tCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCgtggtg                                                         >  2:1667784/1‑149 (MQ=255)
                                                                                               gcgACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGcc                                                       <  2:4730733/149‑1 (MQ=255)
                                                                                               gcgACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGcc                                                       <  1:3183026/149‑1 (MQ=255)
                                                                                                 gACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAg                                                     >  1:3134330/1‑149 (MQ=255)
                                                                                                  aCGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGt                                                    <  1:1989279/149‑1 (MQ=255)
                                                                                                  aCGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGt                                                    >  1:2391376/1‑149 (MQ=255)
                                                                                                  aCGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGt                                                    >  1:3016983/1‑149 (MQ=255)
                                                                                                   cGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTAAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGtt                                                   <  1:4592852/149‑1 (MQ=255)
                                                                                                    gCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTa                                                  <  1:1399695/149‑1 (MQ=255)
                                                                                                     cacaAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTAt                                                 >  1:4280474/1‑149 (MQ=255)
                                                                                                     cacaAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTAt                                                 <  2:2954685/149‑1 (MQ=255)
                                                                                                      acaAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATc                                                >  2:2140408/1‑149 (MQ=255)
                                                                                                        aaaGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATcgc                                              <  2:2269295/149‑1 (MQ=255)
                                                                                                        aaaGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATcgc                                              >  2:2575226/1‑149 (MQ=255)
                                                                                                         aaGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCTCGCACCAGACGCTGGACACCGTGGTGCCAGTTATcgcc                                             >  2:1265461/1‑149 (MQ=255)
                                                                                                          aGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATcgccg                                            >  2:2819952/1‑149 (MQ=255)
                                                                                                            aTCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCgccgcc                                          <  2:3545879/149‑1 (MQ=255)
                                                                                                             tCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCt                                         <  1:1198441/149‑1 (MQ=255)
                                                                                                             tCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTAAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCt                                         >  2:4232096/1‑149 (MQ=255)
                                                                                                                gCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTcc                                      >  2:2336/1‑149 (MQ=255)
                                                                                                                gCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTcc                                      >  2:2421692/1‑149 (MQ=255)
                                                                                                                     gAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAg                                 >  2:1046800/1‑149 (MQ=255)
                                                                                                                       gCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGAt                               >  2:2506018/1‑149 (MQ=255)
                                                                                                                         gTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGc                             >  2:4349704/1‑149 (MQ=255)
                                                                                                                          ttCTTTTACATCCACCAACAGCCATTAAAAGCCGTACGGGTCAAAATCAAGATCATGAATAGCTATATGGTGGCGGTAGGTGAGGTACACATCGCTCACCAGACTCTGGACACCGTGGTGCCAGTTATCTCCGCCTTCCAACAGATGcc                            >  2:1541107/1‑149 (MQ=255)
                                                                                                                          ttCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGcc                            >  1:2085912/1‑149 (MQ=255)
                                                                                                                           tCTTTGTCATCCACCACCAGCCATTAAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCa                           >  2:2831083/1‑149 (MQ=255)
                                                                                                                            cTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAg                          <  1:284500/149‑1 (MQ=255)
                                                                                                                               ggTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCg                       <  2:4110947/148‑1 (MQ=255)
                                                                                                                                gTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGa                      <  2:2020415/149‑1 (MQ=255)
                                                                                                                                gTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGa                      >  1:636787/1‑149 (MQ=255)
                                                                                                                                    tCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCAGTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCt                  >  1:1561929/1‑149 (MQ=255)
                                                                                                                                     ccaccaccaGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTg                 <  1:595130/149‑1 (MQ=255)
                                                                                                                                        ccaccaGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTGATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCg              >  2:3743083/1‑149 (MQ=255)
                                                                                                                                         caccaGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGt             >  1:3471529/1‑149 (MQ=255)
                                                                                                                                           ccaGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCa           <  1:1898144/149‑1 (MQ=255)
                                                                                                                                             aGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCCGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATg         <  2:3243984/149‑1 (MQ=255)
                                                                                                                                             aGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATg         >  2:2991000/1‑149 (MQ=255)
                                                                                                                                               ccATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGt       <  2:2705462/149‑1 (MQ=255)
                                                                                                                                                cATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGtt      >  1:2908041/1‑149 (MQ=255)
                                                                                                                                                  ttCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTcc    >  2:3062302/1‑149 (MQ=255)
                                                                                                                                                    cAAAGCCGGACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCAc  <  2:2284243/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGCCATTTGCCCGGCTGGTTTATGCCGAAGCGATAATCATGACGCGGTACCGGCGTAAAGTTACTGGCAACGATGATTTCGTTACCCTCTTTATCGCGACGCACAAAGATCAGCACCGAGCGTTCTTTGTCATCCACCACCAGCCATTCAAAGCCGTACGGGTCAAAATCCAGTTCATGCATTGCTTTATGGTGGCGGTAGGTGAGGTTCAGATCGCGCACCAGACGCTGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCAC  >  CP009273/3564836‑3565132

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: