Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,566,100 G→T 25.5% R255S (CGC→AGC)  glgB ← 1,4‑alpha‑glucan branching enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,566,1000GT25.5% 730.5 / 190.9 320R255S (CGC→AGC) glgB1,4‑alpha‑glucan branching enzyme
Reads supporting (aligned to +/- strand):  ref base G (116/122);  new base T (46/36);  total (162/158)
Fisher's exact test for biased strand distribution p-value = 3.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 323 total aligned reads displayed.

CGGTTGGCTGATAACCCCAACTGCCATCGAAGGGATGCTCGTTAATGGGCAGTAGTTCGAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCATTTGCGCTTCAAAGGCATAAG  >  CP009273/3565953‑3566249
                                                                                                                                                   |                                                                                                                                                     
 ggTTGGCTGATAACCCCAACTGCCATCGAAGGGATGCTCGTTAATGGGCAGTAGTTCGAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACg                                                                                                                                                    >  2:368471/1‑150 (MQ=255)
  gTTGGCTGATAACCCCAACTGCCATCGAAGGGATGCTCGTTAATGGGCAGTAGTTCGAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGc                                                                                                                                                   >  1:3047762/1‑150 (MQ=255)
     ggCTGATAACCCCAACTGCCATCGAAGGGATGCTCGTTAATGGGCAGTAGTTCGAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAg                                                                                                                                                >  1:1011410/1‑150 (MQ=255)
      gCTGATAACCCCAACTGCCATCGAAGGGATGCTCGTTAATGGGCAGTAGTTCGAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAgg                                                                                                                                               <  1:4682494/150‑1 (MQ=255)
          aTAACCCCAACTGCCATCGAAGGGATGCTCGTTAATGGGCAGTAGTTCGAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCTACGCCAGGAAcc                                                                                                                                           >  1:1537471/1‑150 (MQ=255)
          aTAACCCCAACTGCCATCGAAGGGATGCTCGTTAATGGGCAGTAGTTCGAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCTACGCCAGGAAcc                                                                                                                                           >  1:1624900/1‑150 (MQ=255)
          aTAACCCCAACTGCCATCGAAGGGATGCTCGTTAATGGGCAGTAGTTCGAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAAcc                                                                                                                                           >  1:5132398/1‑150 (MQ=255)
            aaCCCCAACTGCCATCGAAGGGATGCTCGTTAATGGGCAGTAGTTCGAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCTACGCCAGGAACCCa                                                                                                                                         <  1:1726504/150‑1 (MQ=255)
            aaCCCCAACTGCCATCGAAGGGATGCTCGTTAATGGGCAGTAGTTCGAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCa                                                                                                                                         <  2:1038005/150‑1 (MQ=255)
            aaCCCCAACTGCCATCGAAGGGATGCTCGTTAATGGGCAGTAGTTCGAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCa                                                                                                                                         <  1:3717823/150‑1 (MQ=255)
              ccccAACTGCCATCGAAGGGATGCTCGTTAATGGGCAGTAGTTCGAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAgg                                                                                                                                       <  1:5302555/150‑1 (MQ=255)
               cccAACTGCCATCGAAGGGATGCTCGTTAATGGGCAGTAGTTCGAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGt                                                                                                                                      <  1:5038929/150‑1 (MQ=255)
                ccAACTGCCATCGAAGGGATGCTCGTTAATGGGCAGTAGTTCGAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTg                                                                                                                                     <  1:3028653/150‑1 (MQ=255)
                  aaCTGCCATCGAAGGGATGCTCGTTAATGGGCAGTAGTTCGAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGaa                                                                                                                                   >  2:4747713/1‑150 (MQ=255)
                  aaCTGCCATCGAAGGGATGCTCGTTAATGGGCAGTAGTTCGAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGGTGTCGGTGTGGCGACGCCAGGAACCCAGGTGaa                                                                                                                                   >  2:4283877/1‑150 (MQ=255)
                     tGCCATCGAAGGGATGCTCGTTAATGGGCAGTAGTTCGAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACtt                                                                                                                                >  2:3589274/1‑150 (MQ=255)
                        cATCGAAGGGATGCTCGTTAATGGGCAGTAGTTCGAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCTACGCCAGGAACCCAGGTGAACTTCAt                                                                                                                             <  1:4617710/150‑1 (MQ=255)
                        cATCGAAGGGATGCTCGTTAATGGGCAGTAGTTCGAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCAt                                                                                                                             <  1:5011940/150‑1 (MQ=255)
                         aTCGAAGGGATGCTCGTTAATGGGCAGTAGTTCGAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCTACGCCAGGAACCCAGGTGAACTTCATa                                                                                                                            <  2:1353813/150‑1 (MQ=255)
                             aaGGGATGCTCGTTAATGGGCAGTAGTTCGAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATa                                                                                                                        <  1:1405838/150‑1 (MQ=255)
                                ggATGCTCGTTAATGGGCAGTAGTTCGAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCTACGCCAGGAACCCAGGTGAACTTCATAAATagag                                                                                                                     <  1:1216349/150‑1 (MQ=255)
                                ggATGCTCGTTAATGGGCAGTAGTTCGAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATagag                                                                                                                     <  1:3297111/150‑1 (MQ=255)
                                   tGCTCGTTAATGGGCAGTAGTTCGAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGAtt                                                                                                                  <  1:5313682/150‑1 (MQ=255)
                                   tGCTCGTTAATGGGCAGTAGTTCGAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGAtt                                                                                                                  >  2:756324/1‑150 (MQ=255)
                                    gCTCGTTAATGGGCAGTAGTTCGAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCTACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTg                                                                                                                 <  2:212774/150‑1 (MQ=255)
                                    gCTCGTTAATGGGCAGTAGTTCGAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTg                                                                                                                 <  2:3702143/150‑1 (MQ=255)
                                     cTCGTTAATGGGCAGTAGTTCGAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCTACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTgg                                                                                                                >  2:3738113/1‑150 (MQ=255)
                                     cTCGTTAATGGGCAGTAGTTCGAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTgg                                                                                                                >  2:518721/1‑150 (MQ=255)
                                     cTCGTTAATGGGCAGTAGTTCGAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTgg                                                                                                                >  2:520300/1‑150 (MQ=255)
                                      tCGTTAATGGGCAGTAGTTCGAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGgc                                                                                                               <  2:5246868/150‑1 (MQ=255)
                                       cGTTAATGGGCAGTAGTTCGAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCTACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGgcg                                                                                                              <  1:4854215/150‑1 (MQ=255)
                                        gTTAATGGGCAGTAGTTCGAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGgcgc                                                                                                             >  2:487462/1‑150 (MQ=255)
                                          tAATGGGCAGTAGTTCGAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCAt                                                                                                           <  2:4633088/150‑1 (MQ=255)
                                             tGGGCAGTAGTTCGAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCaa                                                                                                        <  2:3132937/150‑1 (MQ=255)
                                                     agtTCGAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTc                                                                                                >  2:4086149/1‑150 (MQ=255)
                                                     agtTCGAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTc                                                                                                >  2:776297/1‑150 (MQ=255)
                                                     agtTCGAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATAGTTGTCGGTGTGGCTACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTc                                                                                                <  1:1939116/150‑1 (MQ=255)
                                                       ttCGAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCCGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGc                                                                                              >  2:1645942/1‑150 (MQ=255)
                                                        tCGAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCt                                                                                             <  2:2030525/150‑1 (MQ=255)
                                                          gAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCttt                                                                                           >  1:4239834/1‑150 (MQ=255)
                                                          gAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCttt                                                                                           <  1:404202/150‑1 (MQ=255)
                                                          gAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCttt                                                                                           <  1:258496/150‑1 (MQ=255)
                                                           aGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCtttt                                                                                          <  2:195032/150‑1 (MQ=255)
                                                             gTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCtttttt                                                                                        <  1:58495/150‑1 (MQ=255)
                                                              tGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTg                                                                                       >  2:3264049/1‑150 (MQ=255)
                                                              tGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTg                                                                                       >  2:4257963/1‑150 (MQ=255)
                                                              tGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGACAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTg                                                                                       >  1:3856781/1‑150 (MQ=255)
                                                               gggTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTgc                                                                                      >  1:40075/1‑150 (MQ=255)
                                                               gggTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTgc                                                                                      <  2:4704256/150‑1 (MQ=255)
                                                                   aaaGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCACTTTTTTGCGctc                                                                                  >  2:972451/1‑150 (MQ=255)
                                                                    aaGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCTACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGctct                                                                                 >  2:4788817/1‑150 (MQ=255)
                                                                    aaGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGctct                                                                                 <  1:3935917/150‑1 (MQ=255)
                                                                     aGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCTACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCtt                                                                                <  1:922286/150‑1 (MQ=255)
                                                                     aGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCTACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCtt                                                                                >  2:1468749/1‑150 (MQ=255)
                                                                         catccatTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGt                                                                            >  1:1272539/1‑150 (MQ=255)
                                                                         catccatTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGt                                                                            >  2:3760707/1‑150 (MQ=255)
                                                                            ccatTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCTACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTg                                                                         >  1:2254129/1‑150 (MQ=255)
                                                                             catTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCTACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGt                                                                        >  2:1536174/1‑150 (MQ=255)
                                                                             catTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCTACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGt                                                                        >  1:3122736/1‑150 (MQ=255)
                                                                              atTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTATCTCAAACATAAATTGTTGTCGGTGTGGCGACGACAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTa                                                                       <  2:3423510/150‑1 (MQ=255)
                                                                                ttAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACa                                                                     >  1:558709/1‑150 (MQ=255)
                                                                                   gCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAAcc                                                                  >  1:115338/1‑150 (MQ=255)
                                                                                    cATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCt                                                                 >  2:4013932/1‑150 (MQ=255)
                                                                                     aTAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCtt                                                                >  2:3102124/1‑150 (MQ=255)
                                                                                     aTAAGGCACCAGTTGATCGGCCAGCTCGAGGTAGCTCAACCAGAAATTGTTGTCGGTGTGTCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGAATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCtt                                                                <  2:121415/150‑1 (MQ=255)
                                                                                      tAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCTACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCttt                                                               <  1:4790063/150‑1 (MQ=255)
                                                                                      tAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACAttt                                                               >  1:2054092/1‑150 (MQ=255)
                                                                                        aGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCTACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTc                                                             <  2:1689174/150‑1 (MQ=255)
                                                                                          gCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCg                                                           <  2:2456801/150‑1 (MQ=255)
                                                                                           cACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCTACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCgg                                                          >  1:4697303/1‑150 (MQ=255)
                                                                                           cACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCgg                                                          >  1:4484264/1‑150 (MQ=255)
                                                                                                 ttGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCTACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGccc                                                    <  1:2307101/150‑1 (MQ=255)
                                                                                                 ttGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGccc                                                    <  2:871044/150‑1 (MQ=255)
                                                                                                 ttGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGccc                                                    <  1:756324/150‑1 (MQ=255)
                                                                                                  tGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCTACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCg                                                   >  1:3418043/1‑150 (MQ=255)
                                                                                                  tGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCg                                                   >  1:2813310/1‑150 (MQ=255)
                                                                                                       ggCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAAt                                              <  2:1790250/150‑1 (MQ=255)
                                                                                                       ggCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAAt                                              >  1:2953370/1‑150 (MQ=255)
                                                                                                        gCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAATTa                                             <  2:4384663/150‑1 (MQ=255)
                                                                                                              tCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGAc                                       <  1:2522870/150‑1 (MQ=255)
                                                                                                               cgcgGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGAcg                                      <  1:3677059/150‑1 (MQ=255)
                                                                                                                  ggTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAATATACGCgg                                   >  2:4108573/1‑150 (MQ=255)
                                                                                                                    tAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGtt                                 >  1:2070209/1‑150 (MQ=255)
                                                                                                                      gCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTc                               >  2:1747917/1‑150 (MQ=255)
                                                                                                                       cTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTcc                              <  2:2966475/150‑1 (MQ=255)
                                                                                                                                aaaTTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCAtt                     >  2:2808908/1‑150 (MQ=255)
                                                                                                                                 aaTTGTTGTCGGTGTGGCTACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCAttt                    <  1:4392731/150‑1 (MQ=255)
                                                                                                                                  aTTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCATTTg                   <  1:4557713/150‑1 (MQ=255)
                                                                                                                                     gttgtCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCATTTgcgc                >  2:3336949/1‑150 (MQ=255)
                                                                                                                                          cGGTGTGGCTACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCATTTGCGCTTCaa           >  1:4150436/1‑150 (MQ=255)
                                                                                                                                               tgGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCATTTGCGCTTCAAAGGCa      <  2:5325138/150‑1 (MQ=255)
                                                                                                                                                ggCTACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCATTTGCGCTTCAAAGGCAt     <  2:3529350/150‑1 (MQ=255)
                                                                                                                                                ggCTACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCATTTGCGCTTCAAAGGCAt     >  2:3464589/1‑150 (MQ=255)
                                                                                                                                                ggCTACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCATTTGCGCTTCAAAGGCAt     <  1:1468749/150‑1 (MQ=255)
                                                                                                                                                ggCTACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCATTTGCGCTTCAAAGGCAt     >  2:2220355/1‑150 (MQ=255)
                                                                                                                                                ggCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCATTTGCGCTTCAAAGGCAt     >  1:4940865/1‑150 (MQ=255)
                                                                                                                                                ggCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCATTTGCGCTTCAAAGGCAt     <  2:2370812/150‑1 (MQ=255)
                                                                                                                                                  cGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCATTTGCGCTTCAAAGGCATaa   <  1:518721/150‑1 (MQ=255)
                                                                                                                                                   gACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCATTTGCGCTTCAAAGGCATAAg  >  2:5035831/1‑150 (MQ=255)
                                                                                                                                                   gACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCATTTGCGCTTCAAAGGCATAAg  <  1:879455/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CGGTTGGCTGATAACCCCAACTGCCATCGAAGGGATGCTCGTTAATGGGCAGTAGTTCGAGGTGGGTAAAGCCCATCCATTTAGCATAAGGCACCAGTTGATCGGCCAGCTCGCGGTAGCTCAACCAGAAATTGTTGTCGGTGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCATTTGCGCTTCAAAGGCATAAG  >  CP009273/3565953‑3566249

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: