Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,931,021 C→A 84.3% A123E (GCG→GAG)  rbsK → ribokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,931,0210CA84.3% 632.2 / 133.0 279A123E (GCG→GAG) rbsKribokinase
Reads supporting (aligned to +/- strand):  ref base C (26/18);  new base A (111/124);  total (137/142)
Fisher's exact test for biased strand distribution p-value = 1.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 280 total aligned reads displayed.

CCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGC  >  CP009273/3930875‑3931166
                                                                                                                                                  |                                                                                                                                                 
   gcagcTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCtgg                                                                                                                                              <  1:3175114/149‑1 (MQ=255)
        tCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGaa                                                                                                                                         >  1:1024821/1‑149 (MQ=255)
          gCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAgc                                                                                                                                       <  2:3892175/149‑1 (MQ=255)
            cACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAgcgc                                                                                                                                     >  2:598490/1‑149 (MQ=255)
              cTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCaa                                                                                                                                   >  1:2780401/1‑149 (MQ=255)
                  tAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTg                                                                                                                               >  2:1753418/1‑149 (MQ=255)
                   aaCATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGa                                                                                                                              >  2:3093100/1‑149 (MQ=255)
                   aaCATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGa                                                                                                                              <  2:4012930/149‑1 (MQ=255)
                    aCATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAg                                                                                                                             >  2:1735165/1‑149 (MQ=255)
                    aCATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAg                                                                                                                             <  1:4711045/149‑1 (MQ=255)
                     cATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGc                                                                                                                            >  2:2708719/1‑149 (MQ=255)
                      aTTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCg                                                                                                                           <  1:4618789/149‑1 (MQ=255)
                       ttGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGt                                                                                                                          >  2:2904610/1‑149 (MQ=255)
                       ttGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGt                                                                                                                          >  1:2035185/1‑149 (MQ=255)
                        tGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTa                                                                                                                         <  2:2960625/149‑1 (MQ=255)
                        tGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTa                                                                                                                         <  2:1318686/149‑1 (MQ=255)
                         gATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTAt                                                                                                                        <  1:972235/149‑1 (MQ=255)
                         gATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTAt                                                                                                                        >  2:788695/1‑149 (MQ=255)
                           tatTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTg                                                                                                                      >  1:1906618/1‑149 (MQ=255)
                            atTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGc                                                                                                                     <  2:3175261/149‑1 (MQ=255)
                            atTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGc                                                                                                                     >  1:1753338/1‑149 (MQ=255)
                            atTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGc                                                                                                                     <  2:484873/149‑1 (MQ=255)
                             ttACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGcc                                                                                                                    >  2:3953643/1‑149 (MQ=255)
                               aCTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCaa                                                                                                                  <  2:419673/149‑1 (MQ=255)
                               aCTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCaa                                                                                                                  <  2:450758/149‑1 (MQ=255)
                                    ggTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGt                                                                                                             >  1:3392116/1‑149 (MQ=255)
                                      tCAGCGTGATCAAAGGCGAATCAACAGGTGAGGCGCGGATTTTTCTTAAAGCCCAAGGTGATAAAGACAACTGTATTAATTCCTGCGCTAATTCAGCCCAATCCCCCGCGCAGCTTGAAGCCCAACGTGAGAGTATATCCAACTCGACa                                                                                                           >  1:3653693/1‑149 (MQ=255)
                                        aGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGc                                                                                                         >  1:3237999/1‑149 (MQ=255)
                                         gCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCa                                                                                                        >  1:3991831/1‑149 (MQ=255)
                                         gCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCa                                                                                                        <  1:2975586/149‑1 (MQ=255)
                                         gCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCa                                                                                                        >  1:1602470/1‑149 (MQ=255)
                                            tGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATTTCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCCGAGCTGGTGGAATCGCAACGTGAGCGTATTTCAAACTCGTCAGCatta                                                                                                     >  1:490060/1‑149 (MQ=255)
                                            tGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCatta                                                                                                     >  1:3947728/1‑149 (MQ=255)
                                            tGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCatta                                                                                                     >  1:2466060/1‑149 (MQ=255)
                                              aTCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCattatt                                                                                                   >  1:3072525/1‑149 (MQ=255)
                                              aTCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCattatt                                                                                                   <  1:2129037/149‑1 (MQ=255)
                                               tCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCAttatta                                                                                                  >  2:1547954/1‑149 (MQ=255)
                                               tCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCAttatta                                                                                                  >  2:617759/1‑149 (MQ=255)
                                                 aaaGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAAt                                                                                                >  1:1449075/1‑149 (MQ=255)
                                                        aaTCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTg                                                                                         <  1:3532705/149‑1 (MQ=255)
                                                        aaTCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTg                                                                                         <  1:3453020/149‑1 (MQ=255)
                                                         aTCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTgg                                                                                        <  2:106814/149‑1 (MQ=255)
                                                             aCAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATc                                                                                    >  2:2122123/1‑149 (MQ=255)
                                                             aCAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATc                                                                                    <  2:660869/149‑1 (MQ=255)
                                                             aCAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATc                                                                                    <  2:3703975/149‑1 (MQ=255)
                                                             aCAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATc                                                                                    >  1:229595/1‑149 (MQ=255)
                                                                ggTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATcacc                                                                                 >  2:3837298/1‑149 (MQ=255)
                                                                ggTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAAAGTCATCGGTATTCATTCCGGCGCTAATTCTTCCCATTACCCGGCGCAGGTTGAAGCGCAACGTGAGCGTATTTCCAACGCGTCAGCATTATTAATTCAGCAGGAATcacc                                                                                 >  1:1321457/1‑149 (MQ=255)
                                                                  tgtgGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATcaccac                                                                               >  1:4558529/1‑149 (MQ=255)
                                                                       cgcTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGaa                                                                          <  1:1682867/149‑1 (MQ=255)
                                                                        gcTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGaaa                                                                         <  1:1276973/149‑1 (MQ=255)
                                                                         cTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAg                                                                        >  1:3846416/1‑149 (MQ=255)
                                                                          tGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAgt                                                                       <  1:1776240/149‑1 (MQ=255)
                                                                          tGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAgt                                                                       <  2:2274278/149‑1 (MQ=255)
                                                                           gATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAgtg                                                                      <  1:4280529/149‑1 (MQ=255)
                                                                           gATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAgtg                                                                      >  2:248092/1‑149 (MQ=255)
                                                                           gATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAgtg                                                                      >  2:223402/1‑149 (MQ=255)
                                                                            aTTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAgtgt                                                                     >  1:2109582/1‑149 (MQ=255)
                                                                               tttGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGAt                                                                  >  2:891856/1‑149 (MQ=255)
                                                                                ttGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATg                                                                 >  1:1156935/1‑149 (MQ=255)
                                                                                  gTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGc                                                               <  1:3036256/149‑1 (MQ=255)
                                                                                   ttAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCa                                                              >  2:3555566/1‑149 (MQ=255)
                                                                                      aTGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAgcg                                                           <  2:727497/149‑1 (MQ=255)
                                                                                        ggCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAgcggc                                                         <  2:3203198/149‑1 (MQ=255)
                                                                                              ggTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAAt                                                   >  2:2538080/1‑149 (MQ=255)
                                                                                                tGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCg                                                 <  1:4392249/149‑1 (MQ=255)
                                                                                                tGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGACAATCg                                                 >  1:4762979/1‑149 (MQ=255)
                                                                                                tGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCg                                                 >  1:4271100/1‑149 (MQ=255)
                                                                                                tGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCg                                                 <  1:1071719/149‑1 (MQ=255)
                                                                                                    aaTGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCa                                             <  2:1968483/149‑1 (MQ=255)
                                                                                                      tGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATc                                           >  1:1040839/1‑149 (MQ=255)
                                                                                                        tCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCaa                                         >  2:4374673/1‑149 (MQ=255)
                                                                                                                 ttCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACt                                >  1:2851889/1‑149 (MQ=255)
                                                                                                                 ttCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACt                                >  1:3690230/1‑149 (MQ=255)
                                                                                                                   cATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTAt                              <  2:2299003/149‑1 (MQ=255)
                                                                                                                    aTGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATc                             >  2:4548535/1‑149 (MQ=255)
                                                                                                                    aTGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATc                             >  2:3738635/1‑149 (MQ=255)
                                                                                                                     tGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCg                            <  2:2819397/149‑1 (MQ=255)
                                                                                                                     tGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCg                            <  1:2991251/149‑1 (MQ=255)
                                                                                                                      gCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATAAAAATAAGACTATCGt                           <  1:2328296/149‑1 (MQ=255)
                                                                                                                         ggCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTgc                        <  1:4566695/149‑1 (MQ=255)
                                                                                                                          gcgcTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTgcg                       >  1:1238486/1‑149 (MQ=255)
                                                                                                                            gcTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCt                     <  1:3492673/149‑1 (MQ=255)
                                                                                                                             cTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCtt                    <  2:1172223/149‑1 (MQ=255)
                                                                                                                             cTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCtt                    >  2:3882715/1‑149 (MQ=255)
                                                                                                                               aaTGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCAGCAAAATAAGACTATCGTTGCGCTTaa                  >  2:4036502/1‑149 (MQ=255)
                                                                                                                                aTGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAAc                 >  2:4128903/1‑149 (MQ=255)
                                                                                                                                  gctgcCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAAccc               <  2:1906605/149‑1 (MQ=255)
                                                                                                                                   ctgcCCTTTCCCCGGCGCGGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACccg              <  1:96662/149‑1 (MQ=255)
                                                                                                                                   ctgcCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACccg              <  1:597876/149‑1 (MQ=255)
                                                                                                                                    tgcCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACccgg             >  2:3232065/1‑149 (MQ=255)
                                                                                                                                    tgcCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACccgg             <  1:3619494/149‑1 (MQ=255)
                                                                                                                                    tgcCCTGTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACccgg             <  2:4057554/149‑1 (MQ=255)
                                                                                                                                     gcCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACccggc            <  1:2607473/149‑1 (MQ=255)
                                                                                                                                     gcCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACccggc            <  1:1095048/149‑1 (MQ=255)
                                                                                                                                           tCCCCGGCGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACccggctccggc      >  1:2178933/1‑149 (MQ=255)
                                                                                                                                           tCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACccggctccggc      <  1:4116276/149‑1 (MQ=255)
                                                                                                                                             cccGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCcggctccggctc    <  2:1343020/149‑1 (MQ=255)
                                                                                                                                             cccGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCcggctccggctc    >  2:1420644/1‑149 (MQ=255)
                                                                                                                                              ccGGCGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTcg   >  1:935947/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGC  >  CP009273/3930875‑3931166

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: ‑